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Abstract 



VNI is a general-purpose Monte-Carlo event-generator, which includes the simulation of lepton-lepton, 
pi \ lepton-hadron, lepton- nucleus, hadron-hadron, hadron-nucleus, and nucleus-nucleus collisions. It uses the 
real-time evolution of parton cascades in conjunction with a self-consistent hadronization scheme, as well as 
the development of hadron cascades after hadronization. The causal evolution from a specific initial state 
(determined by the colliding beam particles) is followed by the time-development of the phase-space den- 
sities of partons, pre-hadronic parton clusters, and final-state hadrons, in position-space, momentum-space 
and color-space. The parton-evolution is described in terms of a space-time generalization of the familiar 
momentum-space description of multiple (semi)hard interactions in QCD, involving 2 — > 2 parton collisions, 
2 — > 1 parton fusion processes, and 1 — > 2 radiation processes. The formation of color-singlet pre-hadronic 
clusters and their decays into hadrons, on the other hand, is treated by using a spatial criterion motivated 
by confinement and a non-perturbative model for hadronization. Finally, the cascading of produced pre- 
hadronic clusters and of hadrons is includes a multitude of 2 — > n processes, and is modeled in paralell to the 
parton cascade description. This paper gives a brief review of the physics underlying VNI, as well as a de- 
tailed description of the program itself. The latter program description emphasizes easy-to-use pragmatism 
and explains how to use the program (including simple examples), annotates input and control parameters, 
and discusses output data provided by it. 
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Originally published in: Computer Physics Communications 104, 70 (1997) 
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1 PROGRAM SUMMARY 



Title of program: VNI 

Computer fo which the program has been designed and others on which it has been tested: IBM RS-6000, 
Sun Sparc, Hewlett Packard UX A-9000 

Operating systems: IBM-AIX, Sun-OS, and any other UNIX operating systems, as well as LINUX. 
Programming language used: Fortran 77 

Memory required to execute with typical data: 2000 kwords 
No. of bits in a word: 64 

No. of lines in distributed program: 25760 lines of main program, plus 49 and 244 lines for two example programs. 



Keywords: Monte Carlo simulation, event generator, QCD kinetic theory, parton cascades, parton coalescence, 
hadronic final states. 

Nature of physical problem: 

In high-energy particle collisions certain phase-space regions can be populated by a large number of quanta, such that 
statistical correlations among them (e.g., in space, momentum, or color) become of essential importance. Examples 
are deep-inelastic lepton-hadron scattering and hadron-hadron collisions in the region of very small Bjorken-a;, or, 
collisions involving heavy nuclei in the central rapidity region. In these cases the produced particles evolve in a 
complicated non-equilibrium environment created by the presence of neighboring ones. The 'deterministic' quantum 
evolution of particle states due to self-interactions (depending only on the particle itself), receives a new 'statistical' 
kinetic contribution due to mutual interactions (depending crucially on the local density). The theoretical basis 
for addressing the solution for the dynamics of such particle systems is a quantum-kinetic formulation of the QCD 
equations of motion, an approximation that combines field-theoretical aspects associated with the renormalization 
group (including well-known resummation techniques) with aspects of transport theory associated with non-equilibrium 
multi-particle dynamics (including important quantum effects beyond the classical level) . 

Method of solution: 

The solution of the underlying quantum-kinetic equations of motion for non-equilibrium multi-particle QCD by Monte- 
Carlo simulation of collisions allowing for a variety of combinations of beam and target particles. To simulate the 
real-time evolution of the collision system in position space and momentum space on the basis of the equations of 
motions, the procedure is four-fold: i) the construction of the initial state including the decomposition of the beam and 
target particles into their partonic substructure, (ii) the evolution of parton cascades including multiple scatterings, 
emission- and fusion-processes, (iii) the self-generated conversion of partons into hadrons using a phenomenological 
model for parton-coalcscence into pre-hadronic clusters and subsequent decay into final-state particles, and (iv) the 
final-state interactions among produced clusters and hadrons including scattering, absorption, resonance production. 

Restriction on the complexity of the problem: 

For very high collision energy {^/s ^ 10 TeV in hadronic collisions, and ^ 5 TeV/nucleon in nuclear collisions) 
numerical inaccuracies due to repeated Lorentz boosts, etc., may accumulate to cause problems, although the program 
uses double precision throughout. 

Typical running time: 

The CPU time for a typical simulation is strongly dependent on the type of beam and target, the magnitude of 
collision energy, as well as on the time interval At chosen to follow an event in its real-time evolution. Examples are 
(for At = 35 fm): a) e+ + at = 100 GcV: 10000 events/hour; b) p p at = 200 GeV: 5000 events/hour; c) 
p + Au at ^/s = 200 GeV/nucleon: 100 events/hour; d) Au + Au at ^/s = 200 GeV/nucleon: 1 event/hour; All of the 
above quotes are approximate, and refer to a typical 133 Mhz or 166 Mhz processor on a modern Power-Workstation 
or Power-PC. 
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2 LONG WRITE-UP 



2.1 Introduction 

VNI Q is the Monte Carlo implementation of a relativistic quantum-kinetic approach |2| to the dynamics 
of high-energy particle collisions, inspired by the QCD parton picture of hadronic interactions fs], |5|. It 
is a product of several years of development in both the improving physics understanding of high-energy 
multiparticle dynamics in QCD, as well as the technical implementation in the form of a computer simulation 
program. The most relevant references for the following are Refs. g 0, |, I, 0, 0], where details of the 
main issues, discussed below, can be found. The puropse of VNI is to provide a comprehensive description 
of particle collisions involving beams of leptons, hadrons, or nuclei, in terms of the space-time evolution 
of parton cascades and parton-hadron conversion. The program VNI is concepted as a useful tool (and 
nothing more) to study the causal development of the collision dynamics in real time from a specified initial 
state of beam and target particles, all the way to the final yield of hadrons and other observable particles. 
The collision dynamics is traced in detail on the microscopic level of quark and gluon interactions in the 
framework of perturbative QCD, supplemented by a phenomenological treatment of the non-perturbative 
dynamics, including soft parton-collisions and parton-hadronization. The generic structure of the simulation 
concept is illustrated in Fig. 1. 

The main strength of VNI lies in addressing the physics of high-density QCD, which becomes an in- 
creasingly popular object of research, both from the experimental, phenomenological interest, and from the 
theoretical, fundamental point of view. Presently, and in the near future, the collider facilities HERA (ep, 
eA7), Tevatron {pp, pA), RHIC and LHC {pp, AA) are able to probe new regimes of dense quark-gluon matter 
at very small Bjorken-x or/and at large A, with rather different dynamical properties. The common feature 
of high-density QCD matter that can be produced in these experiments, is an expected novel exhibition 
of the interplay between the high-momentum (short-distance) perturbative regime and the low-momentum 
(long-wavelength) non-perturbative physics. For example, with HERA and Tevatron experiments, one hopes 
to gain insight into problems concerning the saturation of the strong rise of the proton structure functions at 
small Bjorken-x, possibly due to color-screening effects that are associated with the overlappping of a large 
number of small-x partons. Another example is the anticipated formation of a quark-gluon plasma in RHIC 
and LHC heavy ion collisions, where multiple parton rescattering and cascading generates a high-density 
environment, in which the collective motion of the quanta imust be taken into consideration. In this context 
the most advantageous and novel feature of VNI is the space-time cascade picture that provides a potentially 
powerful tool to study high-density particle systems in QCD, where accounting for the dynamically evolving 
space-time structure of the collisions is most important. The necessity of including space-time variables in 
addition to energy-momentum variables for high-density systems may be common knowledge in the field 
of heavy-ion physics, where the space-time aspect of many-body transport theory is essential ingredient to 
describe nucleus-nucleus collisions, but it is new to many high-energy particle physicists which only recently 
began to acknowleged the necessity to include time and space on top of the commonly used, pure momentum 
space description of, e.g., lepton and hadron collisions. 

Where does VNI fit into the diverse family of modern event generators for physics simulation of particle 
collisions (where 'particle' stands for any beam/target particle from an electron to a heavy nucleus)? 

The Monte Carlo models for high-energy particle collisions that are on the market, may be crudely divided 
into two classes: 

(i) The first class embodies event generators that restrict to "clean" reactions involving lepton or proton 
beams only, and which aim at a high-precision description of experimental tests of QCD's first princi- 



ples. Popular examples are PYTHIA ||, HERWIG [|5|, ARIADNE |4|], LEPTO |3, and ISAJET 
|48|. The common feature of these event generators is a combination of well understood perturba- 
tive QCD parton-shower description and a non-perturbative hadronization prescription to convert the 
partonic final state into hadrons. 



^ The three letters VNI do not mean anything profound. VNI is pronounced "Vinnie", short for "Vincent Le CuCurulIo Con 
GiGinello", the Uttle guy who hkes to hang out with his pals, the quarks ( citcitruUos) and gluons (gigmellos). That is QCD in 
Wonderland, and that is the whole, true story. 
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VNI 




Figure 1: Generic flow-chart of the simulation concept of VNI: It starts with the initial beam target 
particles A and B, decomposing them (except for leptonic A and/or B) in their hadronic constituents 
with partonic substructure, then proceeds through the parallely evolving stages of parton-cascade, parton- 
hadron conversion, fragmentation of bcani/targct remnants, final-state hadron-cascade, and finishes up with 
the final particle yield that reflects observables in a detector. 
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(ii) The second class are event generators that aim to describe "dirty" reactions involving nuclei, much less 
based on first-principle knowledge, but instead rely on phenomenological models to mimic the unknown 
details of the underlying physics. Here the widely used concept is to visualize a nuclear collision in terms 
of nucleon-nucleon collisions on the basis of a constituent valence-quark picture plus string-excitation 



and -fragmentation. Examples for these models are FRITIOF Q, DPM [|g], VENUS [|5T|, RQMD 
1 52 1, and HIJET [^^. Distinct from these is HIJING |54], which also incorporates a perturbative QCD 



approach to multiple minijet production, however it does not incorporate a space-time description. 

With the exception of HERWIG, all of the above Monte Carlo generators utilize some form of string frag- 
mentation phenomenology to model the non-perturbative hadronization and final-state particle production. 
Most commonly used is the Lund string model [^^. HERWIG on the other hand is built on a very dif- 
ferent parton- cluster formation/fragmentation approach, which forms the basis of the hadronization scheme 
developed in VNI. 

What are the shortcomings of the above-mentioned Monte Carlo models with respect to the particle dynamics 
in finite- density regions created by high-energy collisions? 

The high-energy particle physics generators of the first class lack the inclusion space-time variables in 
the dynamics description in both the perturbative QCD parton evolution and the non-perturbative process 
of parton-hadron transition. These models therefore cannot account for statistical interactions due to the 
presence of a finite density of particles close-by in space, such as rescattering, absorption of recombination 
processes. Hence, although these models to large extent use QCD's fundamental quark-gluon degrees of 
freedom, important aspects of parton dynamics and interactions at finite density are left out, because the 
particles are assumed to propagate unscathed in free space. 

The event generators of the second class, on the other hand, mostly do utilize a space-time description, 
however on the level of hadronic degrees of freedom (strings, baryons, mesons and resonances) rather than 
partonic degrees of freedom. For ultra-relativistic nucleus-nucleus collisions the parton approach appears 
to be more realistic than the hadronic or string picture, as it has been realized that short-range parton 
interactions play a major role for heavy ion collisions at collider energies of ^/s ^ 100 GeV, at least during 
the early and most dissipative stage of the first few fm. Here copiously produced quark-gluon mini-jets 
cannot be considered as isolated rare events, but must be embedded in complicated multiple cascade-type 
processes. Thus, the short range character of these interactions implies that perturbative QCD can and must 
be used, and that the picture of comparably large distance excitations of strings or hadronic resonances does 
not apply in this kinematic regime. 

In view of this discussion, the program VNI can be viewed in between the above two classes of Monte 
Carlo models: It provides a kinetic space-time description of parton evolution by utilizing well-developed 
techniques for perturbative QCD simulations at zero density or free space, as the event generators of the 
first class. On the other hand, it applies this concept also to the physics of finite-density particle systems, 
e.g., in collisions involving nuclei, as the event generators of the second class. Comparing VNI with the 
above-mentioned Monte Carlo models, the essential differences and partly new aspects are the following: 

a) the aspects of the space-time evolution of the particle distributions in addition to the evolution in 

momentum space Q, and the concepts of quantum kinetic theory and statistical physics [||, 

b) the self-consistent interplay of coherent (angular ordered) perturbative parton evolution according to the 



DGLAP 1 31] equations, with the fully dynamical cluster-hadronization according to the phenomeno- 



logical model of Ref . . 

c) the microscopic tracing of color degrees of freedom and the effect of color-correlations by using explicit 

color-labels for each parton, which allows to investigate final state interactions in the process of hadron 
formation [p^ . 

d) the diverse advantages of a stochastic simulation technique with which the various particle interaction 

processes are determined by the dynamics itself, through the local density of particles as they evolve 
causally in time. 

e) the statistical many-particle description for general non-equilibrium systems, which allows to study 

thermodynamic behaviour of the bulk matter |55|, such as the evolution of macroscopic energy density. 
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pressure, etc., or the dynamical development of the sytem to thermal/chemical equilibrium in heavy 
ion reactions. 

In summary, the improvement to be expected from VNI for the physics simulation of high-energy par- 
ticle collisions lies clearly in the 'dirty' high-density parton regime, where the space-time aspects are most 
important, and which currently and in the future is of central interest in experiments at, e.g., HERA, RHIC 
and LHC. On the other hand, VNI may also be perceived as a valuable alternative to high-energy event 
generators for the study of 'clean', zero-density collisions as e~^e~ annihilation or pp collisions, where the 
space-time aspects can provide useful additional insight in the collision dynamics for experiments at, e.g., 
LEP or the Tevatron. Recent applications to e~^e~ and ep collisions may be found in Refs. §, |lO|, 0], and 
to heavy-ion colhsions in Refs. [0, 111, 



2.2 General concept 

The central element in the physics description implemented in VNI is the use of QCD transport theory Q 
and quantum field kinetics Q to follow the evolution of a generally mixed multi-particle system of partons 
and hadrons in 7-dimensional phase-space d^rSkdE. Included are both the parton-cascade development 
^, 15, 0] which embodies the renormalization-group improved evolution of multiple parton collisions 



including inelastic (radiative) processes, and the phenomenological parton-hadron conversion model of Refs. 
[P, 10, 11], in which the hadronization mechanism is described in terms of dynamical parton-cluster formation 



as a local, statistical process that depends on the spatial separation and color of nearest-neighbor partons, 
followed by the decay of clusters into hadrons. 

In contrast to the commonly- used momentum-space description, the microscopic history of the 
dynamically-evolving particle system is traced in space-time and momentum space, so that the correla- 
tions of partons in space, time, and color can be taken into account for both the cascade evolution and 
the hadronization mechanism. It is to be emphasized, that the interplay between perturbative and non- 
perturbative regimes is controlled locally by the space-time evolution of the mixed parton-hadron system 
itself (i.e., the time-dependent local parton density), rather than by an arbitrary global division between 
parton and hadron degrees of freedom (i.e., a parametric energy /momentum cut-off). In particular the 
parallel evolution of the mixed system of partons, pre-hadronic clusters, and hadrons, with the relative 
proportions determined by the dynamics itself, is a novel feature that is only possible by keeping track of 
both space-time and energy-momentum variables. 

Probably the greatest strength of this approach lies in its application to the collision dynamics of com- 
plicated multi-particle systems, as for example in collisions involving nuclei {eA, pA and AB)^ for which a 
causal time evolution in position space and momentum space is essential: Here statistical, non-deterministic 
particle interactions are most important, which can only be accounted for by following the time-evolution 
of the particle densities in space and energy-momentum. This approach allows to study the time evolution 
of an initially prepared beam/target collision system in complete phase-space from the instant of collisional 
contact, through the QCD-evolution of parton distributions, up to the formation of final hadronic states. 
It provides a self-consistent scheme to solve the underlying equations of motion for the particle densities as 
determined by the microscopic dynamics. 

The model as a whole consists of four major building-blocks, which are illustrated schematically in Fig. 

2: 

1. The initial state associated with the incoming nuclei involves their decomposition into nucleons and of 

the nucleons into partons on the basis of the experimentally measured nucleon structure functions 
and elastic form- factors. This procedure translates the initial nucleus-nucleus system into two collid- 
ing clouds of virtual partons according to the well-established parton decpomposition of the nuclear 
wavefunctions at high energy Q. 

2. The parton cascade development starts from the initial interpenetrating parton clouds, and involving 

the space-time development with mutual- and self- inter actions of the system of quarks and gluons. 
Included are multiple elastic and inelastic interaction processes, described as sequences of elementary 
2 — > 2 scatterings, 1^2 emissions and 2^1 fusions. Moreover, correlations are accounted for be- 
tween primary virtual partons, emerging as unscathed remainders from the initial state, and secondary 
real partons, materialized or produced through the partonic interactions. 
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3. The hadronization dynamics of the evolving system in terms of parton-coalescence to color-neutral clusters 

is decribed as a local, statistical process that depends on the spatial separation and color of nearest- 
neighbor partons. Each pre-hadronic parton-cluster fragments through isotropic two-body decay into 
primary hadrons, according to the density of states, followed by the decay of the latter into final stable 
hadrons. 

4. The 'afterburner' hadron cascade describes the evolution of produced pre-hadronic clusters and hadrons, 

emerging both from the hadronization of cascading partons, as well as from primary remnant partons 
which represent the fraction of unscathed initial state nucleons. Pre-hadronic clusters can mutually 
rescatter, or scatter off close-by hadrons, before they decay into stable final-state hadrons. Similarly, 
already formed hadrons may deflected by elestic collisions with other hadrons and clusters, or may be 
excited by inelastic collisions and resonance formation/decay. 

Such a pragmatical division, which assumes complex interference between the different physics regimes to 
be negligible, is possible if the respective dynamical scales are such that the short-range (semi)hard parton 
interactions (scattering, radiation, fusion) of perturbative nature, and the non-perturbative mechanism 
of hadron formation (parton-coalescence and cluster-decay), occur on well-separated space-time scales (or 
momentum scales, by virtue of the uncertainty principle). Loosely speaking, the typical momentum scale 
associated with parton collisions, radiative emissions, or parton fusion, has to be larger than the inverse 
'confinement length scale' ~ 1 fm which separates perturbative and non-perturbative domains. Further 
discussion of this condition of validity can be found in, e.g., M, 17|. 



2.3 Equations of motion from quantum kinetics for multi-particle dynamics 



A firm theoretical basis for the above space-time cascade description of a multiparticle system in high-energy 
collisions can be derived systematically from quantum-kinetic theory on the basis of QCD's first principles 
in a stepwise approximation scheme (see e.g., Refs. |2| and references therein). This framework allows to 
cast the time evolution of the mixed system of individual partons, composite parton-clusters, and physical 
hadrons in terms of a closed set of integro-differential equations for the phase-space densities of the different 
particle excitations. The definition of these phase-space densities, denoted by F^, where a = p,c,h labels 
the species of partons, pre-hadronic clusters, or hadrons, respectively, is: 



Fa{r,k) = F^{t,f;E,k) 



dN^ 



d^rd^kdE 



(1) 



where r = r^ = {t,f), k = k^^ = {E,k), and k^ = E'^ — k"^ can be off-shell (space-like k"^ < m?, time-like 
k'^ > vn?) or on-shell (/c^ = m?). The densities (|l|) measure the number of particles of type a at time t with 
position in r + dr, momentum in /c + dk, and energy \n E + dE (or equivalently invariant mass in k"^ + dk"^). 
The Fa are the quantum analogues of the classical phase-space distributions, including both off-shell and on- 
shell particles, and hence contain the essential microscopic information required for a statistical description 
of the time evolution of a many-particle system in complete 7-dimensional phase-space d^rd^kdE, thereby 
providing the basis for calculating macroscopic observables. 

The phase-space densities (||) are determined by the self-consistent solutions of a set of transport equa- 
tions (in space-time) coupled with renormalization-group type evolution equations (in momentum space). 
Referring to Refs. [^, ||] for details, these equations can be generically expressed as convolutions of the den- 
sities Fa of particle species a, interacting with specific cross sections Ij for the processes j. The resulting 
coupled equations for the space-time development of the densities of partons Fp, clusters Fc and hadrons 
Fh is a self-consistent set in which the change of the densities Fa is governed by the balance of the various 
possible interaction processes among the particles. Fig. 3 represents these equations pictorially. For the 
densities of partons, the transport equation (governing the space-time change with r^) and the evolution 
equation (controlling the change with momentum scale A;^), read, respectively. 



d 



FpiiFpiii o 



FpFp, o 



lyp p 



II III i\ 

>pp) 



+ I{p"p"' P) 



FpFp, o 



I{pp p"p") + I{pp 



I{p'p" p) — I{pp p") 



Fp Fp, o 



i{pp 



(2) 
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a) initial state 



beam A 




beam B 



b) parton cascade development 




primary 

. partons 

c) parton-cluster formation, hadronization, hadron cascade development 

central 



beam 
hadrons 




beam 
ladrons 



beam cluster B 

parton clusters 



Figure 2: The four components of the model: a) the initial state constructed in terms of the parton 
distribution of the incoming nuclei; b) the time-evolution of parton cascades in 7-dimensional phase-space 
c) the formation of color neutral clusters from secondary partons emerging from cascading, as well of the 
remnant primary partons from the initial state, followed by the fragmentation of the clusters into final 
hadrons. d) the final-state hadron-cascade evolution 8f produced clusters and hadrons. 



d ^ ^ 

Fp{r, k) = Fp> o I{p' pp") - Fpo I{p p'p") . 



(3) 



For the densities of pre-hadronic clusters and hadrons, the evolution equations are homogeneous to good 
approximation, so that one is left with non-trivial transport equations only, Q ^ Hence, for the evolution of 
the cluster densities Fc, one has 

k^-^ Fc{r, k) = Fp Fp. o i{pp' c) - F^o i{c h) + F^^F^,. o i{c"c"' ^ cc') - F,F,, o i{cc' ^ c"(('i) 



+ F^^Fh' o i{c'h' ch) - F^Fh o i{ch ch!) + 
0, 



(5) 



and similarly, for the evolution of the hadron densities Fh-, the equations read (Fig. 3) 



k,-^^ Ft,{r,k) = F,oi{c^h) + 



Fh, o i{h' -^h) - Fho i{h h') 



+ Fh'>Fy„ o i{h"h"' hh') (6) 



- FhFy o i{hh' h"h"') + F^>Fy o I{c'h' ch) - F^Fh o i{ch ch') 



. 



(7) 



In each convolution F o I oi the density of particles F entering a particular vertex I includes a 

sum over contributing subprocesses, and a phase-space integration weighted with the associated subprocess 
probability distribution of the squared amplitude. Explicit expressions are given in Refs. [S, ^. In (pl)- 



(^) , each convolution F o I of the density of particles F entering a particular vertex / includes a sum over 
contributing subprocesses, and a phase-space integration weighted with the associated subprocess probability 
distribution of the squared amplitude. Explicit expressions are given in Refs. HI, |9|. 

As mentioned before, the equations (||)-® reflect a probabilistic interpretation of the multi-particle 
evolution in space-time and momentum space in terms of sequentially-ordered interaction processes j, in 
which the rate of change of the particle distributions Fa (a = p, c, h) in a phase-space element d^rd'^k is 
governed by the balance of gain (+) and loss (— ) terms. The left-hand side describes free propagation 
of a quantum of species a, whereas on the right-hand side the interaction kernels / are integral operators 
that incorporate the effects of the particles' self- and mutual interactions. This probabilistic character is 
essentially an effect of time dilation, because in any frame where the particles move close to the speed of 
light, the associated wave-packets are highly localized to short space-time extent, so that comparatively 
long-distance quantum interference effects are generally very small. 

The terms on the right-hand side of the transport- and evolution-equations @-(|6|) correspond to one of 
the following categories (c.f. Fig. 3): 

• parton scattering and parton fusion through 2-body collisions, 

• parton multiplication through radiative emission processes on the perturbative level, 

• colorless cluster formation through parton coalescence depending on the local color and spatial con- 
figuration, 

• hadron formation by decay of the cluster excitations into final-state hadrons. 

^ It is worth noting that eq. (^) embodies the momentum space (fc^) evolution of partons through the renormalization of 
the phase-space densities Fp, determined by their change k^dFp{r,k)/dk^ with respect to a variation of the mass (virtuality) 
scale in the usual QCD evolution framework ^ ^ . On the other hand, for pre-hadronic clusters and hadrons, renormal- 
ization effects are comparatively small, so that their mass fluctuations Afc'^/fc^ can be ignored to first approximation, implying 
k^dFc{r,k)/dk^ = k^dFh{r,k)/dk^ = 0. 

It is worth noting that eq. (^) embodies the momentum space (fe^) evolution of partons through the renormalization of 
the phase-space densities Fp, determined by their chaiige k^dFp{r, k)/dk^ with respect to a variation of the mass (virtuality) 
scale k^ in the usual QCD evolution framework 0, nffl. approximation, implying k^dFc{r, k) /dk^ = k^ dFh(r,k) / dk^ = 0. 
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parton-cluster 
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parton-branching 



CLUSTERS: 
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FIG. 3. Graphical representation of the equations (2)-(5) for the particle phase-space densities Fp of partons, Fc of 
pre-hadronic clusters, and Fh of hadrons. The 'dots' represent further processes not explicitely shown, such as cluster-hadron 
absorption, as well as inelastic collisions among clusters and hadrons. 
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• scattering of pre-hadronic clusters or already formed hadrons with other clusters or hadrons. Note: 
this includes also (as indicated by the 'dots') cluster-cluster fusion, cluster-hadron absorption, as well 
as inelastic collisions among clusters and hadrons, in which energy-momentum is transferred into 
excitation. 



The equations ©-(^ reflect a probabilistic interpretation of the multi-particle evolution in space-time 
and momentum space in terms of sequentially-ordered interaction processes j, in which the rate of change of 
the particle distributions Fa {a = p, c, h) in a phase-space element d^rd^k is governed by the balance of gain 
(+) and loss (— ) terms. The left-hand side describes free propagation of a quantum of species a, whereas 
on the right-hand side the interaction kernels I are integral operators that incorporate the effects of the 
particles' self- and mutual interactions. This probabilistic character is essentially an effect of time dilation, 
because in any frame where the particles move close to the speed of light, the associated wave-packets are 
highly localized to short space-time extent, so that comparatively long-distance quantum interference effects 
are generally small for these short-distance interactions (however it does not mean they are unimportant for 
global particle dynamics such as collective behavior). 

It is worth emphasizing that the model goes beyond a purely classical cascadel It incorporates a number 
of quantum mechanical features in the particle dynamics and production, for example: 

• the inclusion of off-shell particles, which treats space- like or time-like particles as energy fluctuations 
with finite life-time given by the dedgee of virtuality oc l/-\/|A?]"; 

• the feature of coherence in parton scatterings which relates momentum transfer and impact param- 
eter of 2-body scatterings and restricts the possible momentum exchange to be consistent with the 
uncertainty relation; 

• the concept of (angular ordered) coherent parton evolution including important gluon interference 
effects is used, which is well known to be an important non-classical property of soft (small-x) partons, 

• particle scatterings and decays (both partonic and hadronic) are treated consistent with the uncertainty 
principle, by requiring a scattering/formation time before particles appear as incoherent quanta in the 
system; 

• particles are produced not at a sharply localized point in space and momentum space (as it would 
be classically, but impossible quantum mechanically); rather their position and momentum at the 
production vertex is smeared with exp[— r^/cr^ — u^p^] (with a = 0.5 fm as inferred from Bose-Einstein 
data from e^e~ collisions), a prescription which imposes again the uncertainty principle and guarantees 
ArAp > h. 

Therefore the particle's densities Fa{r,k) are not just the classical phase-space densities (although they 
are just the number of countable particles per unit phase-space), but incorporate all the above quantum 
mechanical features implicitely. The space-time and energy integrated Fa are in fact observable quantities 
because after all what is seen in a detector are a particle distributions. 



2.4 Scheme of solution in global Lorentz-frame of reference 

In the above kinetic approach to the multi-particle dynamics, the probabilistic character of the transport- 
and evolution equations (§)-(^) allows one to solve for the phase-space densities Fa{r,k) by simulating the 
dynamical development as a Markovian process causally in time. Because it is an 'initial- value problem', 
one must specify some physically appropriate initial condition Fa{to,f, k) at starting time to; such that all 
the dynamics prior to this point is effectively embodied in this initial form of Fa. The set of equations 
(0)-(|6|) can then be solved in terms of the evolution of the phase-space densities F^ for t > to using Monte 
Carlo methods to simulate the time development of the mixed system of partons, clusters, and hadrons in 
position and momentum space |^, ^. 
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With the initial state specified as discussed below, the phase-space distribution of particles at t = Iq = 
can be constructed and then evolved in small time steps At = 0{W~^ fm) forward throughout the par- 
allely evolving stages of parton cascade, parton-cluster formation, and cluster-hadron decays, until stable 
final-state hadrons and other particles (photons, leptons, etc.) are left as freely-streaming particles. The 
partons propagate along classical trajectories until they interact, i.e., collide (scattering or fusion process), 
decay (emission process) or recombine to pre-hadronic composite states (cluster formation). Similarly, the 
so-formed pre-hadronic parton-clusters travel along classical paths until they convert into primary hadrons 
(cluster decay), followed by the hadronic decays into stable final state particles. The corresponding proba- 
bilities and time scales of interactions are sampled stochastically from the relevant probability distributions 
in the kernels / of eq. (0)-(|6|). 

It is important to realize, that the spatial density and the momentum distribution of the particles are 
intimately connected: The momentum distribution continously changes through the interactions and deter- 
mines how the quanta propagate in coordinate space. In turn, the probability for subsequent interactions 
depends on the resulting local particle density. Consequently, the development of the phase-space densities 
is a complex interplay, which - at a given point of time - contains implicitely the complete preceeding history 
of the system. 

It is clear that the description of particle evolution is Lorentz- frame dependent, and a suitable refer- 
ence frame (henceforth called global frame) must be chosen (not necessarily the laboratory frame of an 
experiment). When computing Lorentz-invariant quantities, such as cross sections or final-state hadron 
spectra, the particular choice is irrelevant, whereas for non-invariant observables, such as energy distibu- 
tions or space-time-dependent quantities, one must at the end transform from the arbitrarily-chosen frame 
of theoretical description to the actual frame of measurement. 

For calculational convenience, it is most suitable to choose the global center- of-mass (cm) frame of the 
colliding beam particles, with the collision axis in the z-direction. In this global cm-frame, the incoming 
particles A, B (= lepton, hadron, or nucleus) have four-momenta, 

P^S = (SA,B,0,0,±Pem) 

Ea,b = ^ (s + Ml ± Ml) s = El„, = {P^ + P^f (8) 

Pcm = ^ ^[s - {Ma + Mb)2] [s - {Ma - Mb)^] , 

where Ma^b are the masses of A and B, and so incoming particles then carry well-defined fractions of Pcm, 
having only a non-vanishing longitudinal momentum along the z-axis. Particularly in the case of nuclei A, 
the daughter nucleons Ni = 1, . . . A, have momenta, Pn^ = (0, 0, ±Pcm/A). 

In the following, the global cm-frame of A and B is assumed to be the reference frame, with the initial 
energy-momentum of the collision system given by (^). Furtermore the terms 'hadron', respectively 'nucleon', 
are used to distinguish initial states A + B with A and/or B being a single hadron, respectively a nucleus 
with A {B) nucleons. 



2.5 Initial state 

If one or both beam particles A, B are leptons, they are considered as point-like objects which carry 
the full beam energy, meaning that any QED or QCD substructure of the leptons, as well as initial-state 
photon radiation by the leptons, is neglected. For lepton-lepton annihilation, it is assumed that the colliding 
leptons produce a time-like 7 or boson of invariant mass = -\-q'^ at time t = —Q~^, so that t = 
characterizes the point when the 7 {Z^) decays into a quark-antiquark pair. Similarly, for lepton-hadron 
(nucleus) collisions, the lepton is emitting a space- like virtual 7 of invariant mass = —q^ at time 
t = —Q^^, and hence t = is the point when the 7 hits the hadron (nucleus). 

In the general case of collisions involving hadrons and/or nuclei, the incoming particles A and B are 
decomposed into their parton substructure by phenomenological construction of the momentum and spatial 
distributions of their daughter partons on the basis of the known hadron (nucleon) structure functions and 
elastic hadron (nucleon) form-factors. In the cm-frame, where the two incoming particles A, B (= hadron. 
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nucleus), are moving close to the speed of light, the parton picture is applicable and the parton substructure 
of the hadrons or nucleons can be resolved with reference to some initial resolution scale Qq. This resolution 
scale generally varies with beam-energy and mass number A, B, in that it depends on the typical momentum 
and spatial density of partons as well as on their interaction probability. To be specific, it may be identified 
with the statistically estimated expectation value for the interaction scale of all primary parton-parton 
collisions (that are those in which at least one initial state parton is involved) |]l^]: 

Ql ^ Ql{x,P,A) = A- (-^ + Mk^] ' , (9) 

where {p\)prim is the average relative transverse momentum squared generated by primary parton-parton 
collisions, and Rq = IGeV^^. The pre-factor A" (0 < a < 4/3 a parameter) accounts for the nuclear 
dependence for A,B > 1, and x is the parton's momentum fraction of the mother hadron or nucleon which 
carries momentum P (= Pcm for single hadron, or = Pcm/A for a nucleon in a nucleus A). The scale Qq 
defines the initial point in momentum space above which the system of partons is treated as an ensemble 
of incoherent quanta. The dynamics prior to this point is contained in the initial parton-structure function 
of the mother hadron (nucleon). Clearly, the convention (|9|) yields only an average value dominated by the 
most probable (semi)hard parton collisions with relatively small momentum transfers of a few GeV. However, 
primary parton collisions with a momentum scale ^ Qq, which correspond to relatively rare fluctuations, 
are taken into account individually by the Q^-evolution of the hadron (nucleon) structure functions through 
space-like and time-like radiation processes, discussed later. 

The actual form of the phase-space distribution of partons, eq. (^ is initially to be specified at the 
reference point Qq. It is constructed as the following superposition of the parton distributions in the 
individual hadrons at Qq and at time t = to = 0, the point of collision of A and B: 

(10) 

The right hand side is a sum over all 1 < A'^j < hadrons or nucleons contained in the collision system 
A + B. Each term in the sum is a parton phase-space density given by a convolution of an initial momentum 
distribution P^^ and a spatial distribution Ra^- The subscript a = g, qj, qj (j = 1, . . . , nj) labels the species 
of partons, and Ni refers to the type of the i-th hadron (for nuclei this is just proton or neutron). The 
four-vectors k = = {E, k) and r = r^ = (t, r) refer to the partons. The 3-momenta and initial positions 
of the hadrons or nucleons are denoted P, respectively R. All vectors are understood with respect to the 

global cm-frame, at time t = = 0. The partons' energies E = Ea{k'^,q^) = \Jk'^ + m"^ + take into 
account their initial space-like virtualities q^ < which are distributed around = —Qq/A, under the 

constraint that for each hadron in the initial state, the total invariant mass of the daughter partons must 
equal the mother hadron mass. This mimics the fact that the initial partons are confined inside their parent 
hadrons or nucleons and cannot be treated as free particles, i.e. they have not enough energy to be on mass 
shell, but are off mass shell by a space- like virtuality q^ (eq. ( |l6| ) below). 

As suggestively illustrated in Fig. 2 and discussed in more detail now, the initial state A-\- B involving 
hadrons and/or nuclei, appears in the global cm-frame as two approaching 'tidal waves' of large- a: partons 
(mainly valence quarks), where each 'tidal wave' has an extended tail of low— x partons (mostly gluons and 
sea-quarks) . 

2.5.1 Initial momentum distribution 

For each hadron or nucleon the number of partons, the distribution of the flavors, their momenta and 
associated initial space-like virtualities, are obtained from the function P^"- in ([T0|). Denoting P = Pcm/Nh-, 
with Nh the total number of hadrons (nucleons) in A -|- i?, it is represented in the form 

P^^ik,P-Ql) = (^1) Fi^n^,Qo) Paix) g{k±) 5{P, - P) 5'{P^) (11) 



Fa{r,k)lo^,^ = Y.P^^{k,P-Ql).R^^{k,r,R) 
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with the momentum and energy fractions x 
normalizations 



k,/P, X = E/P = Jx"^ + + m2 + g2)/p2) ^nd the 



E 

a 

E 



dx xF^'{x,Qf^) = 1 







dk"^ d^k 



The physics behind the ansatz (111]) is the following: 



EP^^{k\k,P-Ql) 



d^k± g{k_i 
n''^{QlP) . 



[12) 



(13) 



(i) The functions F^'-{x, Qq) are the scale-dependent measured hadron (nucleon) structure functions with x 
being the fraction of the parent hadron's or nucleon's longitudinal momentum P carried by the parton. 
The transverse momentum distribution g{k±) is specified below. The factor x/x in (O) is included to 
form the invariant momentum integral J d^k/[{2iT)^2E] out of the distribution Pj^' |19|| . In (|l^ ) the 



quantity n^' is the total number of partons in a given nucleon with momentum P = \P\ at Qq. 



(ii) The function Pa{x) in ( pT|) takes into account nuclear shadowing effects affecting mainly soft (small x) 
partons in a nucleus These shadowing effects are evident in the observations of the European Muon 
Collaboration ||2^ as a depletion of the nuclear structure functions at small x relative to those of a free 
nucleon. Several mechanisms have been proposed to explain this nuclear shadowing effect on the basis 
of the parton model |21, |2^. However, here instead the phenomenological approach of Wang and 
Gyulassy |24] is adopted which is based on the following parametrization for the A dependence of 



the shadowing for both quarks and gluons: 

. F^'^^{x,Ql) 



1 + 1.19 Inyl 



Pa{Ra 



1/6 

— 1.5(xo + xl)x'^ + SxqXlX 
1.08(^1/3-1^ 

exp(— x /Xq) 



(14) 



where F^'^ and F, 
nucleon, and xq 



ln{A + 1) 

^ represent the structure function of a nucleus A, respectively those of a free 
0.1, XL = 0.7. The function /^^(r) = 0.1 ( A^^^ 



1 



1 



takes into account 



the impact parameter dependence, with R±^ labeling the transverse distance of a nucleon from its 
nucleus center and Ra the radius of the nucleus. 

(iii) The distribution g{k±) specifies the primordial transverse momenta of partons according to a normal- 
ized Gaussian 



\k±\' 



(15) 



independent of the type of parton or hadron (nucleon). It takes into account the uncertainty of 
momentum (Fermi motion) due to the fact that the initial partons are confined within the nucleons. 
This intrinsic k± can be observed in Drell-Yan experiments where it is found that the distribution is 
roughly independent of s and [^5|. As inferred from these analyses, the width in (|l^) is ko = 0.42 
GeV, corresponding to a mean {k±) 0.38 GeV. 

The scheme (i)-(iii) to sample the fiavor and momentum distribution is carried out independently for 
each nucleon subject to the requirement 



E^^ 




E ^J/i E kz 



Mt 



(16) 



where the summation runs over all partons belonging to the same nucleon as determined by ( [1^ ) and (p^), 
and M/j is the mother hadron (nucleon) mass. With the partons' 4-momenta distributed as outlined above, 

^This feature is optional in the simulation, by default it is switched off. 
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the constramt (16) determines the distribution in the variable = k"^ — = E'^ — k'^ — rn^ < 0, the 
partons' initial space-like virtualities. The resulting distribution in q'^ is a strongly peaked Gaussian with a 
mean value of Qq/^- 



2.5.2 Initial spatial distribution 



The initial spatial distribution of the partons, Ra\ appearing in eq. (10), depends on the magnitude of 
their momenta, the positions of their parent hadrons (nucleons), 
the latter. It is represented as 



as well as on the spatial substructure of 



R^^{p,r,R) = 5-'[R- Rab 



boosted 



(17) 



where the momentum dependence is here purely due to boosting the distributions to the global cm-frame. 
The components of the ansatz (^) have the following meaning: 

(i) The incoming beam particles (a single hadron, or a nucleus) are assumed to have initial cm-positions 
Rab = (iAZyiB/2, -Byifi), corresponding to a chosen impact parameter Bab = iBx,By) and a min- 
imum longitudinal separation AZab- In the case where A and/or B is not a single hadron, but a 
nucleus, the individual nucleons are assigned positions around the centers of the parent nucleus in its 
rest-frame, according to a Fermi distribution for nuclei with mass number A > 12 and a Gaussian 
distribution for nuclei A < 12, 



Hn,{R) = I 

The parameters are c = roA^^^, tq 
mean square radius of the nucleus. 



^ (l + exp[(i?-c)/a])- 

\ 2 _ r .^2/^2] 



47r 0Ffe 



exp 



{A > 12) 
{A < 12) 



1.14 /m, a = 0.545 fm and b 



''A ! 



(18) 



where i?™* is the 



(ii) Next, the partons are distributed around the centers of their mother hadrons or nucleons, still in the 
rest frame of the A, respectively B, with an exponential distribution 



1 u-^ 
An Svr 



exp I 



-vr\ 



(19) 



where v = 0.84 GeV corresponds to the measured elastic formfactor of the mother hadron or nucleon, 

0.81 fm. 



with a mean square radius of i?^* 



(iii) Finally, as indicated by the subscript boosted in eq. (|T^, the positions of the hadrons or nucleons and 
their associated valence quarks are boosted into the global cm-frame of the colliding beam particles A 
and B. The valence quarks then occupy the Lorentz contracted region (Az)„ ~ 2Ra Mh/P, whereas 
the sea quarks and gluons are smeared out in the longitudinal direction by an amount {Az)g^s ~ 
1/kz < 2Ra {Ra = ^0^^/^) behind the valence quarks. This is an important feature of the partons 

As a consequence, the parton 



26, 27, 25]. 



when boosting a hadron or nucleus to high rapidities 
positions are correlated in longitudinal direction with their momenta, as required by the uncertainty 
principle. This leads to an enhancement of the densities of gluons and sea quarks with x < l/(2i2/jM/j) 
proportional to A^^^, because such partons from different nucleons overlap spatially when the nucleons 
are at the same impact parameter. 



2.6 Parton cascade development 

With the above construction of the initial state, the incoming beam particles A and B are decomposed in 
their associated parton content at the initial resolution scale Qq at time t = to = 0, and the subsequent 
dynamical development of the system for t > tQ can now be traced according to the kinetic equations (|2|)- 
(^). In the kinetic approach, the space-time evolution of parton densities may be described in terms of 
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> time > 



Figure 3: Schematical illustration of a typical parton cascade development initiated by a collision of two 

partons a and b. The incoming primary parton a evolves through a space-like cascade from the initial 
resolution scale Qq up to the scale Q^j^ at which it collides with parton b. The outgoing partons c and d 
both initiate time-like cascades which are described as a combination of multiple branchings (or emissions) , 
fusions (or absorptions), or secondary scatterings (or rcscattcrings). A particular branch in a cascade tree 
terminates (locally in space-time), when the partons in that branch recombine with neighboring ones to 
pre-hadronic clusters and their subsequent decay into hadrons. 
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parton cascades. Fig. 4 illustrates a typical parton cascade sequence. It is important to realize that, in 
general, there can be many such cascade sequences, being internetted and simultanously evolving (typical 
for hadron- nucleus or nucleus- nucleus collisions). 

Each cascade can be subdivided into elementary 2 — > 2 scatterings, 1 — > 2 branchings (emissions), and 
2^1 fusions (absorptions). In Fig. 4, a primary parton a that originates from one of the incident nuclei, 
collides with another parton b with some momentum transfer characterized by the scale at the collision 
vertex. The parton a has evolved from the initial scale Qqi at which it was resolved in its parent hadron or 
nucleus, up to by successive space-like branchings. From the scattering of a and b the partons c and 
d emerge, both of which can initiate sequences of time-like branchings. These newly produced partons can 
themselves branch, rescatter, or undergo fusion with other partons. 

In practice the dynamical interplay of 2 — > 2, 2 ^ 1 and 1^2 processes can be simulated as follows: 
Since the time evolution of the parton system is described in small discrete time steps Ai = O(10~'^/m), 
the time scale of an interaction process is compared with Ai to decide whether the interaction occurs within 
this time interval. A virtual parton e* with momentum k'^, k"^ ^ 0, that is produced via a + h — > e*, has 
a short life-time Tg* oc 1/ if its invariant mass | /c^ | is large and it is likely to decay within jt^,* < At 
(the full formulae are derived in [^9|). Therefore, the 2 — > 2 process a + 6— >e* — >c + d preferably occurs. 
On the other hand, if \k'^\ is small, corresponding to a long Te*, parton e* is likely not to decay within the 
time step At and the fusion process a + 6 — > e* rather happens. In this case, e* propagates as a quasi-stable 
particle until in the following time step it has an increased decay probability ||2^. This may then result in 
the 1 — > 2 decay e* c + d. Alternatively, e* might collide with another parton before its decay or may 
propagate freely until the following time step, and so on. In this manner the elementary 2 — > 2, 2 ^ 1, and 
1 — > 2 processes are treated on equal footing and double counting is prevented: either a2— >2, ora2^1, 
and possibly a subsequent 1^2 process can occur, but not both. The relative probability is determined by 
the uncertainty principle, i.e. by relating the momentum scale of the process to the time scale as explained. 

For the collisional processes 2 — > 2 and 2 — > 1, the statistical occurence is determined by the 2-body 
cross-section in Born apprpoximation, with higher-order inelastic corrections effectively included in the 
parton evolution before and after each collision. This parton-evolution is calculated using the well-known 
jet calculus [15, iHl ba sed on the 'modified leading logarithmic approximation' (MLLA) to the QCD evolution 
of hard processes |17|, 3C]. Each individual parton-collision factorizes then in an elementary 2-body collision, 
in which the in- and out-going partons undergo an ordered sequence of elementary branchings a ^ b + c, 
accounting for higher-order radiative corrections to the lowest-order Born approximation. These branchingss 
can be described stochastically as a Markov cascade in position and momentum space. One distinguishes 
initial-state, space-like branchings of the two partons entering a collision vertex, and final-state, time-like 
radiation off the collided partons after the collision. 

The specific feature of the present approach is that, in addition to the definite virtuality and momentum, 
each elementary vertex has a certain space and time position which is obtained by assuming that the partons 
in the cascade propagate on straight-line trajectories in between their interactions. In the MLLA framework, 
the basic properties of both space-like and time-like showers are determined by the DGLAP equations 
[^ ], but with essential differences in time ordering, kinematics and the treatment of infrared singularities 
associated with soft gluon emission. 



2.6.1 Elementary parton-parton collisions 

For the elementary parton scatterings a + 6 ^ c + d, and fusion processes a + 6 — > c*, two distinct classes 
of processes are considered: 

(i) truly perturbative QCD hard parton collisions with a sufficiently large momentum transfer or invari- 

ant mass s; 

(ii) phenomenological treatment of soft parton collisions ^ with low momentum transfer p'^ or invariant 

mass s. 

The motivation for this differentiation of parton-parton collisions is to regulate the singular behavior of 
the collision integrals in (^) which results from the divergence of the associated Born-amplitudes squared 

^The inclusion of soft parton collisions is optional in the program, and by default is switched off. 
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for small momentum transfers (except for s channel processes such as qq annihilation). To render the 
parton-parton cross-sections finite, an invariant hard-soft division scale Pq is introduced such that collisions 
occurring at a momentum scale > Pq are treated as perturbative (semi)hard collisions, whereas for those 
with Q"^ < Pq a, soft, non-perturbative interaction is assumed to occur. That is, the total parton-parton 
cross-section for a collision between two partons a and b is represented as 

^at{s) = V ^ I dQ'l I + ( ^^^^ 1 } , (20) 




Po 



dQ^ 



where 



q2 ^ q2/^ £ f P_L - -^(7^) for t(n) channel .^i) 

1 m? = s for s channel 



S = {Pa+ Pbf , i = {pa - Pcf , U= {pa- Pdf (22) 

That is, the scale of the collision is set by p\ for scatterings a + b ^ a + b, or hy s for annihilation 
processes and a + b — > c+d and fusion processes a + b ^ c. The sum over c, {d) corresponds to summing over 
all possible reaction channels (processes). The specific value of pq is generally dependent on beam energy 
Ecra = {Pa + Ps f of colliding hadrons (nuclei) A + B, as well as on the mass numbers A and B. It is 
treated as a parameter that is determined by existing experimental data. A suitable form is 



Po 



a f2Ecni/GeV\'' 



Po{E,m,A,B) = - ( j;^^ ] (a = 2.0GeV, 6 = 0.27). (23) 



The complementation of both hard and soft contributions renders the parton cross-section finite and 
well defined for all Q^. It must be emphasized that the effect of soft parton collisions on the global dy- 
namics is not essential for collisions involving only leptons and/or hadrons, but plays an important role in 
hadron- nucleus, or nucleus- nucleus collions. The soft parton collisions naturally involve comparably small 
momentum transfer, so that their contribution to transverse energy production is small, but the effect on 
soft particle production is significant. 

In accord with the two-component structure of the cross-section (|20|), these parton collisions are distin- 
guished depending on the momentum transfer: 

(i) Hard collisions above po: 

For the perturbative (semi)hard collisions above pq, the momentum scale is determined by the 
corresponding differential cross-sections 

j~hard -i 



where |Mp is the process-dependent spin- and color-averaged squared matrix element in Born approx- 



imation, published in the literature (see e.g. Refs. |32| and for massive quarks Refs. ||3 



(ii) Soft collisions below po: 

In the case of a non-perturbative soft collision between two partons it is assumed that a very low 
energy double gluon exchange occurs. This provides a natural continuation to the harder collisions 
above pq where the dominant one gluon exchange processes gg — > gg, gq — > gq and qq qq have the 
same overall structure A simple, and physically plausible, form for the soft cross-section that 

continues the hard cross-section for below po down to = 0, may be modelled by introducing a 
screening mass-term /x^ in the nominator of 



^ oc (25) 
dQ^ Q^ + fi^ 

where /x is a phenomenological parameter that governs the overall magnitude of the integrated cj'^"-^* oc 
ln[(pg + /x2)//i2]. The value of /X is not known precisely, but it can be estimated to be in the range 0.3 
- 1 GeV. 
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Notice that both hard and soft scatterings are treated in this approach on completely equal footing. 
That is, with the four momenta of the incoming partons known, the momentum transfer and scattering 
angle are sampled from the respective differential cross-sections da/di (24) and (|25|), assuming azimuthal 
symmetry of the scattering geometry. 



2.6.2 Space-like parton evolution 

Space-like parton cascades are associated with QCD radiative corrections due to brems-strahlung emitted by 
primary space- like partons, which are contained in the initial state hadrons (nuclei), and which encounter 
their very first collision. Thereby such a primary, virtual parton is struck out of the coherent hadron 
(nucleus) wavefunction to materialize to a real excitation. A typical space-like cascade is illustrated in the 
left part of Fig. 4, where the parton oq, originating from the initial hadron (nucleus) state, undergoes 
successive space-like branchings oq aia'i, ai — >■ 0202, a„_i — > a^a'^ to become the parton a = an which 
then actually collides with another parton b. The branching chain proceeds by increasing the virtualities \qf\ 
of the partons aj in the cascade, starting from gq with {q^l ~ Qq/4 up to ~ the space-like virtuality 
of the scattering parton a, where = sets the scale of the collision between partons a and b, and hence 
for the evolution from Qq to Q^. The emitted partons on the side branches with momenta ki on the other 
hand have, due to energy-momentum conservation at the branching vertices, time-like virtualities and each 
of them can initiate a time-like cascade, discussed afterwards. 

The proper inclusion of both collinear, hard and coherent, soft branchings is achieved by describing the 
space-like cascade in both space-time and momentum space by using a so-called angular- ordered evolution 



variable (rather than the virtualities \qj\) [35] 



q]=E]Cj+i, 0+1 = ^^^^ =^ l-cos0o,i+i (0<j<n), (26) 

where qj = {ujj,qj) and kj+i = (E'j+i, £,+1) refer to the j^^ branching aj — ^ aj^laj_^_l with momentum 
assignment qj — > qj+ikj+i (see Fig. 4). The space-like cascade is then strictly ordered in the variable 
Qj^i > Qj, which is equivalent to the ordering of emission angles, Wj^o, j+i' < "^j+i^o, i+2'- 

The space-like cascade terminates with parton an = a entering the vertex of collision with parton b, that 
is, in Fig. 4. The history of parton a is however not known until after it has collided with parton b, 
because it is this very collision that causes the cascade evolution of parton a. Therefore one must reconstruct 
the cascade backwards in time starting from the time of the collisions at the vertex and trace the history 
of the struck parton a back to the initial state at time to = ^-t which it was originally resolved with Qq in its 
hadron (nucleon) mother. The method used here is a space-time generalization of the 'backward evolution 
scheme' [^]. To sketch the procedure, consider the space-like branching qn-i Qnkn which is closest to the 
collision vertex Q^^ in Fig. 4. The virtualities satisfy [16| > IfZ^.il, and qn,Qn-i < (space-like) but 
fe^ > (time-like) . The relative probability for a branching to occur between q'^ and q^ + dq'^ is given by 

q Z ZTT y r [TniXniq ) j 

where Xj = {qj)z/Pz (j = n, n — 1) are the fractions of longitudinal momentum Pz of the initial mother 
hadron (nucleon), with F{rj;xj,q^) = F{rj,qj) the corresponding parton distributions defined by (^, and 
the variables 

En Xn Xn—1 Xn /noN 

^ = p - ' 1 - ^ = -B - — (28) 

specify the fractional energy or longitudinal momentum of parton n and n', respectively, taken away from 
n — 1. The function as/{2TT) 'y(z) is the usual DGLAP branching probability [Q, |l^], with 7(2;) giving 
the energy distribution in the variable z. The last factor T^^\At) in (27) determines the time interval 



in the global cm- frame. At = tn — t„_i, that is associated with the branching process a^-i — > ana]^. It 
accounts for the formation time of a„ from a„_i on the basis of the uncertainty principle: At = Aw/|(jf^|, 
Auj ~ (xn — Xn~i) Pz- A Very simple form is taken here, 

T(^)(At) = 5(^^^^4^Pz - At) . (29) 

V Wnl ) 
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The backwards evolution of the space-hke branching Qn-i Qn + kn is expressed in terms of the probabihty 
that parton a„_i did not branch between the lower bound Qq, given by the initial resolution scale Qq, and 
q^ = fc^ Ql^. In that case, parton n can not originate from this branching, but must have been produced 
otherwise or already been present in the initial parton distributions. This non-branching probability is given 
by the Sudakov form-factor for space-like branchings: 

Sn{xn,f,qh At) = expj-^ r r^'\vi%_^{xn,z,q'^; At)] , (30) 

I a J^-i'l') ) 

where the sum runs over the possible species a = g,q,q oi parton On-i- The upper limit of the p^-integration 
is set by (f ^Q^ab^ associated with the collision vertex with parton b. The limits z± are determined by 
kinematics [Q: z^{q) = Qo/q and z+{q) = 1 — Qo/q. 

The knowledge of the space-like formfactor Sn{xn,q'^,qo', At) is enough to trace the evolution of the 
branching closest to the hard vertex backwards from q^ at t = t^, the time of collision in the global cm- 
frame, qn-i at tn-i = tn — Xn/\qfi\Pz- The next preceding branchings qn-2 — > qn-ik-n-i, etc., are then 
reconstructed in exactly the same manner with the replacements tn tn-i, Xn — > Xn-i, qn Qn-i^ ^^'^ ^o 
forth, until the initial point qQ at to = is reached. 



2.6.3 Time-like parton evolution 

Time-like parton cascades are initiated by secondary partons that emerge either from the side-branches of 
a preceding space-like or time-like cascade, or directly from a scattering or fusion process. Consider the 
time-like cascade initiated by the parton c in the right part of Fig. 4, with momentum k = kn- This parton 
has been produced in the collision a -\- b ^ c -\- d with a time-like off-shellness k'^ ~ Qab- 

Again an angular- ordered (rather than virtuality-ordered) evolution in space-time and momentum-space 
of the cascade is employed to incorporate interference effects of soft gluons emitted along the time-like 
cascade tree of Fig. A, c = Cm ^ Cm-ic^_i, ci — > cqc/q. In contrast to (p6|) , the time- like version of the 



angular evolution variable is |37] 



k] = Ci-i , = ' ' ^ 1 - cos%.i),(,-.i), (m > j > 1) . (31) 

so that the time-like cascade can be described by a fc^-ordered (rather than A;^-ordered) evolution, which 
corresponds to an angular ordering with decreasing emission angles 9jj' > (j_iy . 

Proceeding analogously to the space- like case (c.f. (p^)), the probability dV^^^i-i fo^' the first branching 
after the 'yq vertex, km k^-ik'^-i with /c^_i,/c^_i, is given by the space-time extension [|l|, ^ of the 
usual DGLAP probability distribution |||, |^], 



dV. 



(T) 
m, m- 



i(z,P; At) = ^dz ^^7^_(^.,)_(^.,),(z) T(^)(At) , (32) 



where T(^^(At) is the probabilty that parton m with virtuality k^ and corresponding proper lifetime 
Tm oc l/v^fc^ decays within a time interval At, 

T(^)(At) = 1 - expf-^) . (33) 



tm{k) 

The actual lifetime of the decaying parton m in the global cm- frame is then tm{k) = ^/Tra{k), where 
tq{k) ~ 2>E / {lask"^) for quarks and tg{k) ~ E / {2ask'^) for gluons As before, Fj denotes the local density 
of parton species j = m, m — 1, and Us/ {2t:)^{z) is the DGLAP branching kernel with energy distribution 
7(z). The probability (|3^ is formulated in terms of the energy fractions carried by the daughter partons, 

z = , 1-z = , (34) 
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with the virtuahty km of the parton m related to z and ^ through k^ = + ^m-i + '^^m^i'^ ~ ^.nd 
the argument in the running couphng Os in (^) is = 2z^(\ — z)'^Emi — k\. 

The branching probabihty ( |3^ ) determines the distribution of emitted partons in both coordinate and 
momentum space, because the knowledge of four-momentum and lifetime (or At between successive branch- 
ings) give the spatial positions of the partons, if they are assumed to propagate on straight paths between 
the vertices. The probability that parton m does not branch between and a minimum value Kq = /Iq is 
given by the exponentiation of (|32|), yielding the Sudakov form-factor for time-like branchings: 

- E / - dvZL-iiz,k'^;At)\ , (35) 

which is summed over the species a = g, q, q of parton m — 1. The integration limits /cg and z± are determined 
by the requirement that the branching must terminate when the partons enter the non-perturbative regime 
and begin to hadronize. This condition can be parametrized by the confinement length scale Lc = 0(1 fm) 

with ko^ ^L'"^ = fj,Q, and z^{km) = 1 — Z-{krn) = ^o/\/ 4^m, so that for 2;+(Aio) = •Z-(^o) = 1/^ phase 
space for the branching vanishes. 

The time-like form factor Tmik"^, ^oi At) determines the four- momenta and positions of the partons of a 
particular emission vertex as sketched above for the first branching from k^ at t = tm, the time of production 
of parton c in the global cm-frame, to /^^-i im-i = tm + Em/lk^l- Subsequent branchings are described 
completely analogously by replacing tm t^-i, k^ ^m-i; ^tc. Hence T{k'^, Kq] At) generates a time-like 
cascade as sequential branchings starting from t = at the hard vertex forward in time, until the partons 
eventually hadronize as discussed below. 



2.7 Cluster formation and hadronization 

In view of lack of knowledge about the details of confinement dynamics and the non-perturbative hadroniza- 
tion mechanism, one must rely at present on model-building. In the present approach, the cluster- 
hadronization scheme of Ref. |^ |ll| is employed. This phenomenological scheme is inspired by the 
Mar chesini- Webber model ||3^, however it works in space-time plus color-space. On the other hand it is 
very different from commonly used string-fragmentation models such as the Lund-model [^]. 

In VNI, both the cluster-formation from the collection of quarks and gluons at the end of the perturbative 
phase and the subsequent cluster-decay into final hadrons consist of two components: 

1. The recombination of the secondary time-like partons, their conversion into colorless parton clusters and 

the subsequent decay into secondary hadrons. 

2. The recombination of the primary space-like partons that remained spectators throughout the collision 

development into beam clusters and the fragmentation of these clusters. 

The important assumption here is that the process of hadron formation depends only on the local space-, 
time-, and color-structure of the parton system, so that the hadronization mechanism can be modelled as 
the formation of color-singlet clusters of partons as independent entities (pre-hadrons), which subsequently 
decay into hadrons. This concept is reminiscent of the 'pre-confinement' property [^] of parton evolution, 
which is the tendency of the produced partons to arrange themselves in color-singlet clusters with limited 
extension in both position and momentum space, so that it is suggestive to suppose that these clusters are 
the basic units out of which hadrons form. 

2.7.1 Cluster formation 
(i) Parton clusters: 

Parton clusters are formed from secondary partons, i.e. those that have been produced by the hard 
interaction and the parton shower development. The coalescence of these secondary partons to color- neutral 
clusters has been discussed in detail in Refs. |lO|. Throughout the dynamically-evolving parton cascade 
development, every parton and its nearest spatial neighbour are considered as as potential candidates for 
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a 2-parton cluster, which, if color neutral, plays the role of a 'pre-confined' excitation in the process of 
hadronization. Within each single time step, the probability for parton-cluster conversion is determined for 
each nearest-neighbor pair by the requirement that the total color charge of the two partons must give a 
composite color-singlet state (if necessary by accompanying gluon emission), and the condition that their 
relative spatial distance L exceeds the critical confinement length scale L^- The scale L is defined as the 
Lorentz- invariant distance Lij between parton i and its nearest neighbor j: 

L = Lij = min(Aji,...,Aij,...,Ai„) , (36) 



where Ajj = ^^ {fij)}) + {fij)"^ + {rij)^ + with r^- = r^ — r^ , and the probability for the coalescence of 

the two partons i, j to form a cluster is modelled by a distribution of the form 

Tlij-,c oc (^1 - exp(-AFLij)^ ~ 1 - exp ( ^^ Z^l^^ < - ' 

and Hij^c = (1) if Lij < Lq {Lij > L^)- Here l\F is the local change in the free energy of the system that 
is associated with the conversion of the partons to clusters, and the second expression on the right side is 
a parametrization in terms of Lq = 0.6 fm and L^ = 0.8 fm that define the transition regime. As studied 



in Ref. |1C], the aforementioned color constraint, that only colorless 2-parton configurations may produce 
a cluster, can be incorporated by allowing coalescence for any pair of color charges, as determined by the 
space-time separation Lij and the probability (|3^), however, accompanied by the additional emission of a 
gluon or quark that carries away any unbalanced net color in the case that the two coalescing partons are 
not in a colorless configuration. 

(ii) Beam clusters: 

If one or both of the colliding beam/target particles A and B were a hadron or a nucleus, the one, re- 
spectively two beam clusters are formed from the spectator partons that represent the receding beam/target 
remnants of the original particles A and B. More precisely, the remaining fraction of the longitudinal mo- 
mentum and energy that has not materialized or been redirected and harnessed during the course of the 
collision, is carried by those primary partons of the initial-state hadrons or nuclei, which remained spectators 
throughout. In the present approach these partons maintain their originally assigned momenta and their 
space-like virtualities. Representing the beam remnants of A and/or B, they may be pictured as the co- 
herent relics of the original hadron (or nucleus) wavefunctions. Therefore the primary virtual partons must 
be treated differently than the secondary partons which are real excitations that contribute incoherently 
to the hadron yield. In the global cm-frame, the primary partons are grouped together to form a massive 
beam cluster with its four-momentum given by the sum of the parton momenta and its position given by 
the 3-vector mean of the partons' positions. 



2.7.2 Hadronization of clusters 
(i) Parton clusters: 

For the decay of each parton cluster into final-state hadrons, the scheme presented in detail in Refs. Q 
is employed: If a cluster is too light to decay into a pair of hadrons, it is taken to represent the lightest 
single meson that corresponds to its partonic constituents. Otherwise, the cluster decays isotropically in its 
rest frame into a pair of hadrons, either mesons or baryons, whose combined quantum numbers correspond 
to its partonic constituents. The corresponding decay probability is chosen to be 

Ii,^h = T,{E,,ml) N ^ exp , 38 

Jmu rn-^ V "2-0/ 

where is a normalization factor, and the integrand is a Hagedorn spectrum that parametrizes quite 
well the density of accessible hadronic states below which are listed in the particle data tables, and 
mo = rn-jT- In analogy to (|33|), 7^ is a life-time factor giving the probability that a cluster of mass mf. decays 
within a time interval At in the global cm-frame, 

%{E„ml) = 1 - exp(- ) , (39) 
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with the Lorentz-boosted hfe time tc = JcTc ~ Ec/m^. In this scheme, a particular cluster decay mode is 
obtained from ( p8|) by summing over all possible decay channels, weighted with the appropriate spin, flavour, 
and phase-space factors, and then choosing the actual decay mode acording to the relative probabilities of 
the channels. 

(ii) Beam clusters: 

The fragmentation of the beam clusters containing the spectator partons mimics in the present model 
what is commonly termed the 'soft underlying event', namely, the emergence of those final-state hadrons 
that are associated with the non-perturbative physics which underlies the perturbatively-accessible dynamics 
of the hard interaction with parton shower fragmentation. 

In the spirit of the Marchesini- Webber model |^5|, a version (suitably modified for the present purposes) 
of the soft hadron production model of the UA5 collaboration |4^, which is based on a parametrization of 
the CERN pp collider data for minimum-bias hadronic collisions. The parameters involved in this model 
are set to give a good agreement with those data. 

Soft hadron production is known to be a universal mechanism that is common to all high-energy 
collisions that involve beam hadrons in the initial state, and that depends essentially on the total energy- 
momentum of the fragmenting final-state beam remnant. Accordingly, one may assume that the fragmenta- 
tion of the final-state beam cluster depends solely on its invariant mass M, and that it produces a charged- 
particle multiplicity with a binomial distribution p2| , 

(i7T I U\n 



T{n + h) {n/ky 
nT{k) {1 + n/ky 



P{n) = ' ' ' ^. , (40) 



where the mean charged multiplicity n = n{M'^) and the parameter k = k{M'^) depends on the invariant 
cluster mass R according to the following particle data parametrization 



n(M^) = 10.68 {my-^^^ - 9.5 k{M^) = 0.029 In(M^) - 0.064 . (41) 



Adopting the scheme of Marchesini and Webber the fragmentation of a beam cluster of mass M 
proceeds then as follows: First, a particle multiplicity n is chosen from (^0|), and the actual charged particle 
multiplicity is taken to be n plus the modulus of the beam cluster charge. Next, the beam cluster is split 
into sub-clusters (91^2), {Q2Q3), ■ ■ ■ (Qi = u,d), which are subsequently hadronized in the beam cluster rest 
frame, in the same way as the parton clusters described above. To determine the sub-cluster momenta, the 
following mass distribution is assumed, 

P{M) = c (M - 1) exp [-a{M - 1)] , (42) 

with c a normalization constant and a = 2 GeV~^, resulting in an average value of (M) ~ 1.5 GeV. The 
transverse momenta are taken from the distibution 



P{p±) = c p± exp 



bJpl + M2 



(43) 



with normalization c' and slope parameter 6 = 3 GeV~^, and the rapidities y are drawn from a simple flat 
distribution P{y) oc const, with an extent of 0.6 units and Gaussian tails with 1 unit standard deviation 
at the ends. Finally, all hadronization products of the sub-clusters are boosted from the rest frame of the 
original beam cluster back into the global cm-frame. 



Notice that in this model M fluctuates statistically, as a result of fluctuations of the initial-state parton configuration in 
the incoming hadrons (or nuclei), as well as due to the fluctuating number of remnant partons during the space-time evolution. 
Hence the distribution (E^) and the mean multiplicity vary from event to event. This is in contrast to the original UA5 
model, in which the fixedbeam energy \/s/2 controls the energy dependence of soft hadron production. 
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2.8 Cluster/hadron cascade development 



The 'afterburner' hadron cascade development is concepted in in close analogy with the parton cascade 
evolution. For a description of the physics ingredients and the technical details, the interested reader 
is referred to Ref. [^]. There are two sources for producing primary hadrons, namely, first, the pre- 
hadronic parton-clusters (those emerging from coalescence of materialized, interacted partons), and second, 
the beam/target fragmentation clusters (those resulting from reassembling of the spectator partons that 
have not interacted at all). Both these classes of clusters and hadrons participate in the hadron cascading. 
The clusters or hadrons propagate along classical paths until they either scatter off or absorb other clusters 
or hadrons, or form resonant states with subsequent decays. At each step the newly produced particles 
can reinteract themselves. The collisions among clusters and hadrons are treated on equal footing, and 
the corresponding interaction probabilities are calculated from the 'additive quark model' ||5^ within which 
one can associate a specific cross-section for any pair of colliding hadrons according to their valence quark 
content. Cluster-cluster and cluster-hadron collisions are straightforward to include in this approach, since 
each cluster that emerges from parton coalescence has a definite quark content, depending on the fiavor of 
the mother partons. Again, each newly produced hadron becomes a 'real' particle only after a characteristic 
formation time At^ = j/M^ depending on their invariant mass M/j and their energy through 7 = E/M^. 
Before that time has passed, a hadron must be considered as a still virtual object that cannot interact 
incoherently until it has formed according to the uncertainty principle. 



3 PROGRAM DESCRIPTION 



3.1 The package VNI 

The program package VNI is a completely new written code, using only fragments of its precedessors VNI- 
1.0 (Ref. ||6|) and VNI-2.0 (Ref. 0). The current program is a first release of a long term project that aims 
at unifying the simulation of a wide class of high energy QCD processes, ranging from elementary particle 
physics processes, e.g., e"'"e~ -fragmentation, deep inelastic ep-scattering, via hadronic collisions such as pp 
{pp), to reactions involving (heavy) nuclei, e.g. deep inelastic eA scattering, pA, or AA (AB) collisions. The 
types of collision processes that are available and are discussed in more detail in Section 3.4 below, include 
collisions involving e^, /i^ for leptons, p, n, vr. A, S . . . for hadrons, and De, He, O, . . . , Au, Pb, U for nuclei. 
Some of these collision processes have of course been investigated before; here they are revisited within the 
space-time picture of parton cascades and cluster dynamics in order to provide an alternative method for 
Monte Carlo simulation and also to check the consistency with the wide literature on both experimental 
measurements and theoretical predictions. 



3.2 Overall structure of the program 

The program package VNI is written entirely in Fortran 77, and should run on any machine with such a 
compiler. The program VNI adopts the common block structure, particle identification codes, etc., from the 



Monte Carlos of the Lund family, e.g., the PYTHIA program |44]. In addition, various program elements 
of the latter are used in the sense of a library, however modified to incorporate the additional aspects 
of the space-time and color label information. VNI also employs certain (significantly modified) parts of 
the HERWIG Monte Carlo |^5|, used for the final decays of clusters into hadronic resonances and stable 
particles. It is important to stress that the correspondence of VNI with PYTHIA and HERWIG is rather 
loose and should not lead to naive identification, in particular since VNI is taylored for a larger set of 
degrees of freedom that is necessary to describe the space-time dynamics and the information on color 
correlations. Nevertheless, by truncating the additional information of the particle record of VNI, a direct 
interface to the particle records of the Lund Monte Carlo PYTHIA, and to the program HERWIG, is actually 
straightforward. Instructions are given below in Section 3.10. 

The package VNI consists of the following files: 



(1) the documentation vni.ps (this one) 
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(2) the program VNI: vni.f 

(3) three include files vnil.inc, vni2.inc, vniS.inc 



Further useful, but supplementary software (see http: / / rhic.phys.columbia.edu / rhic/vni ) : 

(4) various example programs for applications of VNI. 

(5) a simple, portable histogram package vnibook.f . 

(5) a program to insert 2-particle probes into VNI, vniprob.f . 

(5) a program to study Bose-Einstein physics and HBT, vnihbt.f . 

(6) the PYTHIA program pythia.f, to serve as interface with VNI if desired 

The program vni.f is the Fortran source code that contains all the subroutines, functions, and block 
data, as categorized above. All default settings and particle data are automatically loaded by including the 
three include files vnil.inc - vni3.inc. The example programs, as well as the histogram package vnibook.f, 
the programs vniprob.f vnihbt.f and pythia.f are not integral part of the program, i.e. are distributed 
supplementary as useful guiding tools. 

The program source code vni.f is organized in 9 parts, in which subroutines and functions that are related 
in their performance duties, are grouped together (a summary list of all components with a brief description 
of their purpose is given in Appendix E): 

part 1: Main steering routines VNIXRIN (initialization), VNIXRUN (event generation), and VNIXFIN 
(finalization) , plus other general-duty subroutines. VNIXRIN and VNIXFIN are only to be called 
once, while VNIXRUN is to be called anew for each event, in which the particle system is evolved in 
discrete time steps and 7-dimensional phase-space. 

part 2: Initialization routines. They include, e.g. the process dependent event-by-event initialization of 
the chosen type of collision process, and various other slave routines associated with the initial set up 
of an event. 

part 3: Package for parton distributions (structure functions). A portable and self-contained collection of 
routines with a number of structure function parametrizations, plus a possible interface to the CERN 
PDFLIB, plus various slave routines. 

part 4: Evolution routines for the parton cascade / cluster-hadronization simulation, including both 
(process-dependent) particle evolution routines and (process-independent) universal routines for space- 
like and time-like shower evolution, for parton scatterings, parton recombination to clusters, and cluster 
decays to hadrons and other final state particles. 

part 5: Diverse utility routines and functions performing the slave-work for the part of perturbative parton 
evolution. 

part 6: Diverse utility routines and functions performing the slave- work for the part of cluster fragmentation 
to hadrons. 

part 7: Evolution routines for the 'afterburner' hadron cascade as a self-contained integrated package. This 
part was developed by R. Longacre, the Wizard of Poquott. 

part 8: Event study and analysis routines for accumulating, analysing and printing out statistics on ob- 
servables. 

part 9: Dummy routines and system dependent utilities. 
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3.3 Special features and machine dependence 

Important to remark are the following special features, being not strictly common in standard Fortran 77: 

1. The program uses of a machine- dependent timing routine, called VNICLOC, which measures the elapsed 

CPU time []. Two alternatives are provided: a) the TIM EX routine of the CERN library, in which 
case the latter must be linked with the program, or, b) the MCLOCK routine which is specific 
to IBM-AIX compilers. If neither of these options are desired, or does not work in the users local 
environment, the lines involving calls for TIMEX {MCLOCK) calls can easily commented out. ^. 

2. All common block variables and dimensions are defined globally within the include files vnil.inc, vni2.inc, 

vni3.inc. This greatly simplifies for instance possible modifications (e.g., enlarging or decreasing) of 
global common block arrays. Consequently, all subroutines contain INCLUDE statements for these 
include files, and therefore when compiling and linking the program, the files vnil.inc, vni2.inc, and 
vni3.inc must be properly accessible. 

3. Double precision is assumed throughout since version 4.0. For applications at extremely high energies, 

single precision for any real variable starts to become a problem. The main source of numerical 
precision-loss arises from multiple Lorentz boosts of particles throughout the simulation. Therefore 
the code has been converted now to full double precision. 

4. Implicit variable-type assignment is assumed, such that variables or functions beginning with letters I - 

M are of integer*^ type, and variables or functions with first letters A - H and O - Z are of type reahS 
(double precision). 

5. SAVE statements have been included in accordance with the Fortran standard. Since most ordinary 

machines take SAVE for granted, this part is not particularly well tried out, however. Users on 
machines without automatic SAVE are therefore warned to be on the lookout for any variables which 
may have been missed. 



3.4 The main subroutines 

There is a minimum of three routines that a user must know: 

1. VNIXRIN for the overall initialization of a sample of collision events, 

2. VNIXRUN for the subsequent actual generation of each event, 

3. VNIXFIN for finishing up simulation. 

These three routines, which are briefly descibed below, are to be called in that order by a user's program as 
is exemplified by the example program in Section 3.10. 

SUBROUTINE VNIXRIN (NEVT , TFIN , FRAME , BEAM , TARGET , WIN) 
Purpose: to initialize the overall simulation procedure. 

NEVT : integer specifying the number of collision events to be simulated for a selected physics process. 

TFIN : final time (in fm) in the global Lorentz frame up to which each collision event is followed. 

^ The usage of the time measurement is not necessary for the performance of the program as-a-whole, ahhough it is a nice 
convenience. 

* The subroutine VNICLOC is hsted at the very end in vm.f. 
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FRAME : a character variable used to specify the global Lorentz frame of the experiment. Uppercase 

and lowercase letters may be freely mixed. 

= 'CMS' : colliding beam experiment in CM frame, with beam momentum in +z direction and target 
momentum in —z direction. 

= 'FIXT': fixed target experiment, with beam particle momentum pointing in +z direction. 
= 'USER' : full freedom to specify frame by giving beam momentum in PSYS(1,1), PSYS(1,2) 
and PSYS(1,3) and target momentum in PSYS(2,1), PSYS(2,2) and PSYS(2,3) in 
commonblocks VNIRECl and VNIREC2. Particles are assumed on the mass shell, and 
energies are calculated accordingly. 

= 'FOUR' : as 'USER', except also energies should be specified, in PSYS(1,4) and 
PSYS(2,4), respectively. The particles need not be on the mass shell; effective 
masses are calculated from energy and momentum. (But note that numerical precision 
may suffer; if you know the masses the option 'FIVE' below is pref errable . ) 
= 'FIVE' : as 'USER', except also energies and masses should be specified, i.e the 
full momentum information in PSYS(1,1) - PSYS(1,5) and PSYS(2,1) - PSYS(2,5) should 
be given for beam and target, respectively. Particles need not be on the mass 
shell. Space-like virtualities should be stored as —y/—m? . Four-momentum and mass 
information must match (see Appendix C for details) . 

BEAM, TARGET : character variables to specify beam and target particles. 

Uppercase and lowercase letters may be freely mixed. An antiparticle may be denoted 
either by (tilde) or (bar) at the end of the name. It is also possible to 
leave out the charge for neutron and proton. 
= 'e^' : electron e~ . 

— 'e+' : positron e"*". 
= 'mu— ' : muon //~ . 

= 'mu+' : antimuon . 

= 'gamma': real photon 7 (not yet implemented) . 

= 'pi+' : positive pion 7r+ . 

= 'pi—' : negative pion 11" . 

= 'piO' : neutral pion tt" . 

= 'nO' : neutron n. 

= 'riD' : antineutron n. 

= 'p+' : proton p. 

— 'p— ' : antiproton p. 
= 'LambdaO': A baryon. 
= 'Sigma+' : S"*" baryon. 

— 'Sigma—' : baryon. 
= 'SigmaO' : E° baryon. 
= 'Xi— ' : H~ baryon. 

— 'XiO' : 2° baryon. 

= 'Omega—': Q~ baryon. 

= 'De' : ?De (Deuterium) nucleus . 

'He' : fHe (Helium) nucleus . 
= 'Ox' : (Oxygen) nucleus. 

= 'Al' : ?1A1 (Aluminium) nucleus. 
= 'Su' : f|s (Sulfur) nucleus. 
= 'Ca' : 20^^ (Calcium) nucleus. 
= 'Cu' : 29CU (Copper) nucleus. 
= 'Kr' : ifiKr (Krypton) nucleus . 
= 'Ag' : Ij^Ag (Silver) nucleus. 

'Ba' : 56 Ba (Barium) nucleus . 
= 'Wt' : ifW (Tungsten) nucleus. 
= 'Pt' : \fFt (Platinium) nucleus. 
= 'Au' : ra^Au (Gold) nucleus. 
= 'Pb' : if^Pb (Lead) nucleus. 
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= 'Ur' : g^^^U (Uranium) nucleus. 

WIN : related to energy of system (in GeV) . Exact meaning depends on FRAME, i.e. for 
FRAME='CMS' it is the total 

energy ^/s of system, and for FRAME='FIXT' it is the absolute 3-momentum P = \fp^ of 
beam particle. 



SUBROUTINE VNIXRUN(IEVT,TRUN) 

Purpose: to generate one event of the type specified by the VNIXRIN. This is the main routine, which 
administers the overall run of the event generation and calls a number of other routines for specific tasks. 

lEVT : integer labeling the current event number. 

TRUN : time (in fm) TRUN < TFIN in the global Lorcntz frame up to which the space time evolution 
of the current event lEVT is carried out. If TRUN < TFIN, then the particle record contains the 
history up to this time, and repeated calls resume from TRUN of the previous call (c.f. Section 3.10). 



SUBROUTINE VNIXFINO 

Purpose: to finish up overall simulation procedure, to write data analysis results to files, evaluate CPU time, 
close files, etc.. This routine must be called after the last event has been finished (c.f. Section 3.10). 

For a deeper understanding of the physics routines and their connecting structure, I refer to Appendix 
A for a brief description. Generally, for each of the included collision processes, there is an initialization 
routine that sets up the initial state, and an evolution routine that carries out the time evolution of the 
particle distributions in phase-space, starting from the initial state. Within the evolution routines, a number 
of universal slave routines (i.e. independent of the process under consideration) perform the perturbative 
parton evolution in terms of space-like branchings, time-like branchings, and parton collisions, as well as the 
cluster hadronization. 

Finally, the program provides some useful pre-programmed event-analysis routines which collect, analyze 
and print out information on the result of a simulation, and give the user immediate access to calculated 
particle data, spectra, and other observables. A more detailed description of their duties, and how to call 
them can be found in Appendix B. 
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Table 1: Available collision processes. 



IPRO 


type of collision process 


1 


^ 7/Z0 ^ 2 — jets hadrons 


2 


^ Zq^ W~^W~ ^ 4 - jets hadrons 


3 


+ h —> jets — > hadrons 


4 


^ + A ^ jets hadrons 


5 


I + h ^ I + jets + X ^ hadrons 


6 


t + A ^ I + jets + A — > hadrons 


7 


h + I ^ I + jets + X ^ hadrons 


8 


h + h' ^ jets + X ^ hadrons 


9 


h + A^ jets + X ^ hadrons 


10 


A + I ^ I + jets + X ^ hadrons 


11 


A + h ^ jets + X ^ hadrons 


12 


A + A' ^ jets + X — >■ hadrons 



3.5 The physics processes 

The program is structured to incorporate different classes of particle collision processes in a modular manner. 
The general classes that are part of the program are summarized in Table 1. All of the generic collision 
processes in Table 1 allow a further specification of beam and/or target particle, .e.g., I = e, fi for leptons, 
h = p,n,7r, . . . for hadrons, or ^ = De, He, O, . . . for nuclei. They are summarized in Table 2. 

As explained before, the actual evolution of any beam/target-particle collision is simulated on microscopic 
level of the system of partons, clusters and hadrons, and their interactions. For the parton interaction 
processes, a selection of elementary hard/soft 2^2 scatterings and 2^1 fusions are included in VNI, 
and future extensions are planned. In addition there are the standard space-like and time-like 1 — > 2 



Table 2: Available beam and target particles. 



leptons I 


hadrons h 


nuclei A 


e± 




fDe 




n,n 


2HC 




P,P 


16q 
8 ^ 




AO 


32 c 
16^ 






40 p„ 






64 p 




d- 


84 
36 






108 A o- 
47 -^g 






Bel 






^fwt 






195 Pi 

78 






}fAn 
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Table 3: Elementary partonic subprocesses. 



Class 


TOT TT~> 


type of subprocess 


a) 2 ^> 2 : 


1 


Qi.qj - 


qiqj 




2 


mi - 


qkqk 




Q 
O 


qm 


99 




4 


qm - 


91 




5 


qm - 


77 




6 


qi9 - 


qi9 




7 


qi9 - 


qa 




8 


99 — 


qkqk 




9 


99 - 


99 




10 


soft scattering 


h) 9 —>■ ^ ■ 

U 1 ? J- . 


± ± 


qm 






12 


qi9 - 


^qt 




13 


99 - 


-9* 


c) 1 ^ 2 (space-like) : 


21 


9*- 


> qm 




22 


q*- 


^ qi9 




23 


9*- 


99 




24 


qt- 


^ qa 


d) 1 ^ 2 (time-like) : 


31 


9*- 


> qm 




32 


qt- 


^ qi9 




33 


9*- 


99 




34 


qt- 


^ qa 



branching processes, which, when combined with the elementary tree processes, provide the usual parton 
cascade method of including higher order, real and virtual, corrections to the elementary tree processes. The 
parton-cluster formation and hadronization scheme includes a number of 2-parton recombination processes 
and the decay of the formed clusters into final state hadron. 

It is possible to select a combination of partonic subprocesses to simulate. For this purpose, all subpro- 
cesses are numbered according to an I SUB code. The list of allowed codes is given below. In the following g 
denotes a gluon, qi represents a quark of flavour i = 1, . . . ,nf, i.e. for n/ = 6 this translates to d, u, s, c, b, t. 
A corresponding antiquark is denoted iji. The notation 7 is for a real photon, i.e. on shell. An asterix * 
denotes an off-shell parton. 



3.6 The particle record 

Each new event generated is in its entity stored in the commonblocks VNIRECl and VNIREC2, which thus 
forms the event record. Here each particle that appears at some stage of the time evolution of the system, 
will occupy one line in the matrices. The different components of this line will tell which particle it is, from 
where it originates, its present status (fragmented/decayed or not), its momentum, energy and mass, and 
the space-time position of its production vertex. The structure of the particle record VNIRECl /VNIREC2 
follows closely the one of the Lund Monte carlos [0, |4^, ^ , employing the same overall particle classification 
of particle identification, status code, etc., however with important differences concerning color- and space- 
time degrees of freedom.. 

Note: When in the following reference is made to certain switches and parameters MSTV or PARV, 
these are described in Sec. 3.7 below. 
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PARAMETER (NV=100000) 

C0MM0N/VNIREC1/N,K(NV,5) ,L(NV,2) 
C0MM0N/VNIREC2/P(NV,5) ,R(NV,5) ,V(NV,5) 

Purpose: Contains the event record, i.e. the complete list of all partons and particles in the current event. 

N : number of lines in the K, L, P, R, and V arrays occupied by the current event. N is continuously 
updated as the definition of the original configuration and the treatment of parton cascading, cluster 
fragmentation and hadron production proceed. In the following, the individual parton/particle num- 
ber, running between 1 and A'^, is called I. The maximum A'^ is limited by the dimension NV of the 
arrays. 

K{I,1) : status code (KS), which gives the current status of the parton/particle stored in the line. The 
ground rule is that codes 1-10 correspond to currently existing partons/particles, while larger codes 
contain partons/particles which no longer exist. 
= : empty line. 

= 1 : an undecayed particle or an unfragmented parton. 

= 2 : an unfragmented parton, which is followed by more partons in the same color-singlet parton 
subsystem. 

= 3 : an unfragmented parton with special color-flow information stored in if (7, 4) and if (7, 5), such 
that color connected partons need not follow after each other in the event record. 
= 4 : an unfragmented parton 1, associated with an externally injected 2-particle probe, i.e. qq, gg, 
qj (only in conjunction with the package VNIEVPR). 

= 5 : an unfragmented parton 2, associated with an externally injected 2-particle probe, i.e. qq, gg, 
q'j (only in conjunction with the package VNIEVPR). 

= 6 : a final state hadron resulting from decay of clusters that have formed from materialized (inter- 
acted) partons. 

= 7 : a final state hadron resulting from soft cluster decay of beam/target remnants that have formed 
from left-over (non-interacted) initial state partons. 

= 8 : a final state hadron resulting from bound-state formation of qq pairs as externally inserted 
probes (only in conjunction with the package VNIEVPR). 

= 9 : a final state hadron resulting from jet fragmentation of q, q, or g as externally inserted probes 
(only in conjunction with the package VNIEVPR). 

= 10 : a prehadronic cluster that has been reconstructed from its original hadronic decay products if 
those have not passed through their formation time yet. 
= 11 : a decayed particle or a fragmented parton, c.f. = 1. 

= 12 : a fragmented jet, which is followed by more partons in the same color-singlet parton subsystem, 
c.f. = 2. 

= 13 : a parton which has been removed when special color-flow has been used to rearrange a parton 

subsystem, cf. = 3. 

= 14 : as = 4, but decayed or fragmented. 
= 15 : as = 5, but decayed or fragmented. 

= 16 : an unstable and removed hadron resulting from decay of clusters that formed from materialized 
(interacted) partons, c.f. = 6. 

= 17 : an unstable and removed hadron resulting from soft cluster decay of beam/target remnants 
formed from left-over (non-interacted) initial state partons, c.f. = 7. 
= 18 : as = 8, but decayed or fragmented. 
= 19 : as = 9, but decayed or fragmented. 

= 20 : as = 10, but decayed or fragmented. 

< : Particle entries with negative status codes (e.g. -1, -6, -7) parallel in their meaning those from 
above, but refer to particles which are only virtually present and become active only after a certain 
formation time, upon which their status code is set to the appropriate positive value. 

K{I,2) : Flavor code (KF) for partons, hadrons and electromagnetic particles included in the current 
version of VNI. The particle code is summarized in Tables 3-9. It is based on the flavor and spin 
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classification of particles, following the 1988 Particle Data Group numbering conventions. Nuclei are 
attributed a non-standard code that is internal to VNI, given by KF = 1000000 + 1000 • Nneutrons + 
Nprotons- A negative KF code, where existing, always corresponds to the antiparticle of the one listed 
in the Tables 3-9. 

K(I,3) : line number of parent particle or jet, where known, else 0. Note that the assignment of a particle 
to a given jet of a jet system is unphysical, and what is given there is only related to the way the event 
was generated. 

K(I,4) : Special color-flow information (for internal use only) of the form K{I,4) = MST\{2)-ICFR + 
ICTO, where ICFR and ICTO give the line numbers of the partons from which the color comes and 
to where it goes, respectively. 

-K'(/,5) : Special color-flow information (for internal use only) of the form K[I,5) = M.STY(2)-JCFR + 
JCTO, where JCFR and JCTO give the line numbers of the partons from which the anticolor comes 
and to where it goes, respectively. 

L{I,1) : color label (L = 1, . . . , NC) of a parton, where NC is number of colors specified by the parameter 
MSTV(5). Is = for antiquarks and all non-colored particles. 

L{I,2) : anticolor label (L = 1, . . . , NC) of a parton. Is = for quarks and all non-colored particles. 

P(/,l) : Px, momentum in the x direction in GcV. 
P{I,2) : py, momentum in the y direction in GeV. 
P{I,3) : Pz, momentum in the z direction in GeV. 
P(J,4) : E, energy, in GeV. 

P(J,5) : m = p^, mass in GeV. For partons with space-like virtualities, i.e. where = —m? > 0, its value 
is P(/,5) = -v^^=-Q. 



R{I,1) : current x position in frame of reference in 1/GeV. 
R{I,2) : current y position in frame of reference in 1/GeV. 
R{I,3) : current z position in frame of reference in 1/GeV. 

R{I,4) : time (in 1/GeV) of active presence in the system, i.e. tpres = t — tprod- Is set equal to PARW(2) 

at time of production tprod of the particle and then increased by PARW(2) in each timcstcp until the 
particle fragments or decays. If ( 0, giving the remaining formation time for virtually produced, but 
not yet real particles. 

P(/,5) : number and type of interactions of a particle, and is equal to MSTV(3)^-A^C0-|- MSTV(3)-A^S'i?-|- 
NTB, where - for partons - NCO is the number of 2-body collisions, NSB of space-like branchings, 
and NTB of time-like branchings, undergone up to current time. For clusters, hadrons and unstable 
particles, NSB is zero and NTB counts the number of 2-body decays. 

V{I,1) : X position of production vertex, in 1/GeV. 
V{I,2) : y position of production vertex, in 1/GeV. 
V(/,3) : z position of production vertex, in 1/GeV. 
V(J,4) : time t of production, in 1/GeV. 
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V(J,5) : encodes origin of production as 10-MSTV(4)-7MO + 10 • ITM + OR, where IMO is the direct 

mother, ITM the time of production, and OR the generation of particle IP in a cascade tree, i.e. the 
genetical origin with respect to the production vertex of its original 'Ur'-mother. (This keeps track of 
the production history even if the decayed mother has been removed and allows to follow the genealogy 
of the system.) 

Special cases: 

For the simulation of collisions involving hadrons or nuclei in the initial state, additional initial state 
information on the status and specific variables of " primary partons" of the incoming hadron (nucleus) are 
stored as follows. Note however, that after the first interaction of a "primary parton", this initial state 
information is lost and the assignments for the components of the arrays R, V are the same as above. 

K{I,4) : KF flavour code of hadron (nucleon), which is mother of the primary parton. 

K{I,^) : KF flavour code of nucleus from which a primary parton originates. (For hadrons K{I,A) and 
K{I,5) are equal). 

R{I,5) : Is equal to —\x\, where x is the longitudinal momentum fraction of primary parton that it carries 
off its mother hadron. 

V(/,4) : encodes location IB in the particle record of a primary sea-quark's brother, i.e. the sea-antiquark 
belonging to the same vacuum polarization loop, as MSTV(4)^ ■ IB+MSTY{A)-ITM + OR, analogous 
to V{1, 5) above. For primary valence quarks and gluons it is 0. 

Additional remark The arrays K{1, 3) — K{I, 5) and the arrays P, R, and V may temporarily take special 
meaning other than the above, for some specific internal use. 

In Section 3.10 below, a typical event listing of the particle record is printed, which serves as an example 
for the organization of the particle record, exhibiting the information for momentum space: It lists the 
particles contained in VNIREC1/VNIREC2 at a certain point during the evolution of an event, where KS 
and KF give status and flavor codes, C and A the color and anticolor index, P the energy-momentum-mass 
in GeV for each particle that appeared at some point in the event history. An analogous listing for the 
position space particle record (not printed here) contains the spatial coordinates and time of active presence 
for each particle, as well as the particle's rapidity and transverse momentum with respect to the 2;-axis 
(jet-axis, or beam-axis). Note: both listings can be obtained at any point by calling the routine VNILIST, 
described in Appendix B. 



Table 4: Quark, lepton and gauge boson codes. 



KF 


Name 


PriiiUxl 


KF 


Name 


Priiil (;d 


Kl' 


Name 


Priiilcd 


1 


d 


d 


11 


e~ 


e- 


21 


g 


g 


2 


u 


u 


12 




nu_e 


22 


7 


gamma 


3 


s 


s 


13 


^" 


mu- 


23 


ZO 


ZO 


4 


c 


c 


14 


I^M 


nu_mu 


24 


W+ 


W+ 


5 


b 


b 


15 


T 


tau- 


25 




HO 


6 


t 


t 


16 


Ur 


nu_tau 


26 






7 






17 






27 






8 






18 






28 






9 






19 






29 






10 






20 






30 







33 



3.7 The general input and control pEtrameters 

The VNIDAO common block contains the status code and parameters that regulate the performance of the 
program. All of them are provided with sensible default values, so that a novice user can neglect them, and 
only gradually explore the full range of possibilities. 

C0MM0N/VNIDA0/MSTV(200) ,PARV(200) ,MSTW(200) ,PARW(200) 

Purpose: to give access to status code and parameters that regulate the performance of the program. If the 
default values, denoted below by (D=. . . ), are not satisfactory, they must in general be changed before the 
VNIXRIN call. Exceptions, i.e. variables that can be changed for each new event, are denoted by (C). 



MSTV(200), PARV(200): control switches and physics parameters 



1. General: 



MSTV(l) : (D=100000) number of lines available in the commonblock VNIREC1/VNIREC2. Should 
always be changed if the dimensions of the K, P, R, V arrays are changed. 

MSTV(2) : (D=25000) is used to trace internal color flow information in the array K. 

MSTV(3) : (D=10000) is used to store number and type of interactions of a particle in the array R of 
VNIREC1/VNIREC2. 

MSTV(4) : (D=1000) is used to trace origin of a particle's production in the array V of 
VNIREC1/VNIREC2. 

MSTV(5) : (D=3) number of colors Nc included in the simulation. 

MSTV(6) : (D=6) number of quark flavors Nf included. 



Table 5: Special codes. 



KF 



Printed 



Meaning 



90 
91 
92 
93 
94 
95 
96 
97 
98 
99 
100 



CMsoft 
cluster 
Beam-REM 
Targ-REM 
CMshower 
SPHEaxis 
THRUaxis 
CLUS j et 
CELL j et 



Center-of-mass of beam/target remnant system before soft fragmentation 
Cluster from coalescence of secondary partons 
Cluster from remnant primary partons of beam particle 
Cluster from remnant primary partons of target particle 
Four-momentum of time-like showering system 
Event axis found with VNISPHE 
Event axis found with VNITHRU 
Jet (cluster) found with VNICLUS 
Jet (cluster) found with VNICELL 
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Table 6: Meson codes, part 1. 



KF 


Name 


Printed 


KF 


Name 


Printed 


211 


7r+ 


pi+ 


213 


P+ 


rho+ 


311 


K" 


KO 


313 


K*" 


K*0 


321 


K+ 


K+ 


323 


K*+ 


K*+ 


411 


D+ 


D+ 


413 


D*+ 


D*+ 


421 


D" 


DO 


423 




D*0 


431 


D,+ 


D_s+ 


433 


d:+ 


D*_s+ 


511 


B» 


BO 


513 




B*0 


521 


B+ 


B+ 


523 




B*+ 


531 


B° 


B_sO 


533 




B*_sO 


541 


Bc+ 


B_c+ 


543 




B*_c+ 


111 


tt" 


piO 


113 


p" 


rhoO 


221 


n 


eta 


223 




omega 


331 


1 

V 


eta' 


333 




phi 


441 




eta_c 


443 




J/psi 


551 


Vh 


eta_b 


553 


T 


Upsilon 


661 


Vt 


eta_t 


663 


e 


Tlieta 


130 


Kl 


K_LO 








310 


Ki 


K_SO 









MSTV(7) : (D=l) energy and momentum conservation options. Since both the initial generation of 
initial hadronic or nuclear parton system and the final cluster-hadron decay scheme require occasional 
reshuffling of 4-momentum, an energy and/or momentum imbalance can occur. This can be corrected, 
if desired, by renormalizing particle energies and momenta. 
= : no explicit conservation of any kind. 

= 1 : energy conservation is enforced at the end of each event (i.e. once only per event): particles 
share energy imbalance compensation according to their energy (roughly equivalent to boosting event 
to CM frame). 

= 2 : as =1, but now energy conservation is enforced at the end of each time step in each event (i.e. 

continously throughout an event). 

= 3 : as =1, but in addition also 3-momentum conservation the end of each event (i.e. once only per 
event), however at the cost of some particles shifted off-mass shell. 

= 4 : as =3, but in addition also 3-momentum conservation the end of each timestep in each event 
(i.e. continously throughout an event). 

MSTV(8) : (D=0) type of particles to be included in energy and momentum conservation options under 
MSTV(7). 

= : all living particles are included. 
= 1 : only partons. 
= 2 : only hadrons. 

MSTV(9) : - not used - 

MSTV(IO) : (D=llllll) initial seed for random number generation. 

MSTV(ll) : (D=5000000) maximum number of collision events A + B, after which program is forcibly 
terminated. 

MSTV(12) : (D=2500) maximum number of time steps per collision event, after which a new collision 
event is generated. When increased, the dimension of the arrays in commonblock VNITIM need to 
be altered accordingly. 
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Table 7: Meson codes, part 2. 



KF 


Name 


Printed 


KF 


Name 


Printed 


10213 


bi 


b_l+ 


10211 


+ 




a_0+ 


10313 


K? 


K_10 


10311 


K(*° 


K*_00 


10323 


1 


K_l+ 


10321 


K*+ 




K*_0+ 


10413 




D_l+ 


10411 


D*+ 


D*_0+ 


10423 


D'> 


D_10 


10421 




D*_00 


10433 




D_ls+ 


10431 


Dr,+ 


D*_Os+ 


10113 


bi' 


b_10 


10111 




a_00 


10223 


h? 


li_10 


10221 


fo° 


f_00 


10333 


hf 


h'_10 


10331 


r/0 
i-Q 


f '_00 


10443 


h?c 


h_lcO 


10441 


x!lc 


chi.OcO 


20213 


a+ 


a_l+ 


215 


a+ 

2 


a_2+ 


20313 


Kf 


K*_10 


315 




K*_20 


20323 


k:+ 

1 


K*_l+ 


325 


k;+ 


K*_2+ 


20413 




D*_l+ 


415 




D*_2+ 


20423 




D*_10 


425 




D*_20 


20433 




D*_ls+ 


435 




D*_2s+ 


20113 


a? 


a_10 


115 




a_20 


20223 


f? 


f_10 


225 


fO 


f_20 


20333 




f '_10 


335 


r/O 
l2 


f '_20 


20443 




chi.lcO 


445 


X'2'c 


chi_2c0 


30443 


i,' 


psi' 








30553 


^/ 


Upsilon' 









MSTV(13) : (D=0) choice of time grid TIME{I) and time increase TINC{I). Gcnerically, t{i) = 
t{i - 1) + [/(i) - f{i - 1)], with i = 0,..., MSTV{12), and t{0) = to =PARV(12), tf =PARV(13). 
= : power law increase f(i) = a ■ (if, where a =PARV(14) and b =PARV(15). 
= 1 : logarithmic increase f{i) = a ■ In(z)'', where a =PARV(14), 6 =PARV(15). 

MSTV(14) : (D=0) Option to select or switch on/off certain physics processes I SUB in the simulation, 
by specifying the array MSUB (c.f. Sec. 3.8 below), where MSUB{ISUB) = (1) switches a 
subprocess off (on). 

= : the relevant processes arc initialized automatically. 
= 1 : user selection of included processes is required. 

MSTV(15) : (D=0) Option to smear out classical (sharply localized) positions and three-momenta of 
particles at their production vertex, according to Gaussian distribution exp[—r'^/a^ —p'^/b'^]-, where 
the widths are a = PARV(16) in fm, and b = PARV(17) in GeV. 
= : switched off, i.e. no smearing. 

= 1 : smearing of particle position according to Gaussian with width PARV(16). 

= 2 : as = 1, plus smearing of particle momentum according to Gaussian with width PARV(17). 

MSTV(16) : (D=0) Option to use PYTHIA for generating the initial state particle system rather than 
VNI. Presently this applies only to Z+Z" coUisions, i.e., IPRO =1 or 2 in Table 1. NOTE: The 
subroutine package vnipythia.f must be linked in order to use this option. Also, in this case, one must 
comment out the dummy routines PYINIT and PYEVNT at the very end of vni.f. 
— : switched off. 
= 1 : switched on. 
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Table 8: Baryon codes. 



KF 


Name 


Printed 


KF 


Name 


Printed 








1114 


A" 


Delta- 


2112 


n 


nO 


2114 


AO 


DeltaO 


2212 


p 


P+ 


2214 


A+ 


Delta+ 








2224 


A++ 


Delta++ 


3112 


s- 


Sigma- 


3114 


s*- 


Sigma*- 


3122 




LambdaO 








3212 




SigmaO 


3214 




Signia*0 


3222 




Sigma+ 


3224 




Signia*+ 


3312 




Xi- 


3314 




Xi*- 


3322 


s° 


XiO 


3324 




Xi*0 








3334 




Omega- 


4112 




Signia_cO 


4114 




Signia*_cO 


4122 


A+ 

c 


Lainbda_c+ 








4212 


c 


Siffma_c+ 

o 


4214 




Sigma*_c+ 


4222 


S++ 

c 


Signia_c++ 


4224 


c 


Siffma* C++ 


4132 


c 


Xi_cO 








4312 


■J/0 
' — 


Xi'_cO 


4314 


— *0 

' — 'c 


Xi*_cO 


4232 




Xi_c+ 








4322 




Xi'_c+ 


4324 




Xi*_c+ 


4332 




Omega.cO 


4334 




Oiiiega*_cO 


5112 




Sigina_b- 


5114 




Sigma*_b- 


5122 




Lambda_bO 








5212 




Sigma_bO 


5214 


^b 


Sigma*_bO 


5222 




Sigma_b+ 


5224 




Sigiiia*_b+ 



Table 9: Nucleus codes. 







Priiil cd 


KF 




Printed 


1001001 
1002002 
1008008 
1014013 
1016016 
1020020 
1035029 


fDe 
^He 

fo 

f^Al 

32 q 
16'-' 

40 Pa 
2o*-^a 

iCu 


de 
he 
ox 
al 
su 
ca 
cu 


1061047 
1110074 
1117078 
118079 
1126082 
1146092 


i?^Ag 

238 TT 
92 


ag 
wt 
pt 
au 
pb 
ur 
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PARV(l) - PARV(ll) : - not used - 



PARV(12) : (D=0. fm) choice of initial point of time ti = PARV(12) when evolution of collision event 
begins. 

PARV(13) : (D=0. fm) choice of final point of time tj = PARV(13). If = 0., then it is set automatically 
to r/M£;(MSTW(2)), where MSTW(2) is the selected number of time steps. Otherwise, if > 0., it 
specifies tj such that the minimum min(tj, t^ax) is taken as the point when the evolution of a collision 
event stops, where tmax = TIME{MSTY{12)). 

PARV(14) : (D=0.05) prefactor a of the time-increment function /, see MSTV(13). 

PARV(15) : (D=l.l) exponent b of the time-increment function /, see MSTV(13). For longer collision 
time scales as in hadron-nucleus or nucleus nucleus collisions, it is automatically scaled by PARV(15) 
to 1.5 X PARV(15) to extend the time range. 

PARV(16) : (D=0.1 fm) For MSTV(15) = 1, width parameter a=PARV(16) of Gaussian exp[-r'^ / a'^] for 
smearing out classical (sharply localized) positions of particles at their production vertex. (Reasonable 
values are a = 0.1 — 1 fm, c.f. Zhang et. al, iiuc-th/9704041 ). 

PARV(17) : (D=0.1 GeV) For MSTV(15) = 2, width parameter 6=PARV(17) of Gaussian exp[-p'^ /b'^] 
for smearing out classical (sharply localized) momenta of particles at their production vertex. (Note, 
unertainty relation requires PARV(17) ~ 1/PARV(16) GeV/fm). 



2. Initial state of collision system: 



MSTV(20) : (D=0) choice of model for initial parton distributions. 
= : standard structure functions are used. 

= 1 : other parametrization, e.g. a la McLerran et al. (not implemented yet). 

MSTV(21) : (D=l) choice of nucleon structure-function set (c.f. MSTV(22)). 
= 1 : GRV94LO (Lowest order fit). 
= 2 : GRV94HO (Higher order: MSbar fit). 
= 3 : GRV94DI (Higher order: DIS fit). 
= 4 : - not used - 

= 5 : CTEQ2M (best higher order MSbar fit). 

= 6 : CTEQ2MS (singular at small x). 

= 7 : CTEQ2MF (flat at small x). 

= 8 : CTEQ2ML (large A). 

= 9 : CTEQ2L (best lowest order fit). 

= 10 : CTEQ2D (best higher DIS order fit). 

MSTV(22) : (D=l) choice of nucleon structure-function library. 

= 1 : the internal one with parton distributions chosen according to MSTV(21). 

= 2 : the PDFLIB one with MSTV(21)=1000 • NGROUP + NSET (requires PDFLIB to be linked). 

MSTV(23) : (D=l) choice of pion structure-function set (c.f. MSTV(24)). 
= 1 : Owens set 1. 
= 2 : Owens set 2. 



MSTV(24) : (D=l) choice of pion structure-function library. 

= 1 : the internal one with parton distributions chosen according to MSTV(23). 

= 2 : the PDFLIB one with MSTV{2'i) = WOO- NGROUP + NSET (requires PDFLIB to be linked). 
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MSTV(25) : (D=2) choice of minimum Q-value Qo used in parton distributions of initial hadronic or 

nuclear state. 

= : fixed at Qo =PARV(23). 

= 1 : determined by the average momentum tranfer of primary parton collisions Qq = (Qprim), with 
initial value Qo =PARV(23) before the first event. 

= 2 : as =1, but in addition varying with total energy y/s of the overall collision system, with Qq = 

max[Qo, f • {Eh/GeV)''], where Qo =PARV(23), a =PARV(24), b =PARV(25), and E,, = ^/i/hadron 
the average energy per hadron (nucleon). 

MSTV(26) : (D=l) option to choose x-dependent Qin in parton densities. 
= : switched off, Qo according to selection of MSTV(25) is used. 

= 1 : varying with Bjorken-x of struck initial state partons, as well as with total energy ^/s of collision 
system, with Qo = Q{x) where 1/Q^(a;) = I/Qq + Qo/{xEhadf , and Eh = -y/s/hadron the average 
energy per hadron (nucleon). 

MSTV(30) : (D=0) sub-selection of initial state for the following collision processes A + B, where specifi- 
cation of beam/target particles B is not unambigous, because collision may proceed via alternative 
channels (c.f. Tabic 1): 

a) X ^ hadrons: 

= : e^e~ '^/Z^ qq 

- 1 : e+e- ^ W+W' -^qq'lu 

- 2 : e+e" ^ W+W' qq' q'q 

MSTV(31) : (D=0) initial separation along the collision (z-)axis of beam and target particle in the overall 
CM frame. 

= : fixed at Az =PARV(31). 

> : shifted towards minimal separation, such that beam and target particles almost touch with 
Az =MSTV(31)-PARV(31). 

MSTV(32) : (D=2) selection of impact parameter b transverse to the collision {z-) axis of beam and 
target particle. 

= : fixed at 6 = Max[PARV(32),PARV(33)]. 

= 1 : randomly sampled in between bmin 

< b < 

bmaxi with bmin 

=PARV(32) and b„iax =PARV(33). 

= 2 : as 1, but with brain ~ 

PARV(32)-y^;;^, and b^ax = PARV(33)- V^^;^, where and is the 
non-diffractive crossection for given beam/target and beam energy. 

MSTV(33) : (D=2) Choice of nuclear matter distribution for collisions involving one or two nuclei, i.e. 
of nucleons' positions within nucleus. 
= : uniform distribution with 'sharp edge'. 

= 1 : Gaussian distribution for light nuclei A < 12, and 'smeared edge' distribution for intermediate 
and heavy nuclei ^ > 12 (c.f. PARV(27)). 

= 2 : Gaussian distribution for A < 12, and Fermi distribution for A> 12 (c.f. PARV(28)). 

MSTV(34) : (D=l) spatial distribution of initial partons within their parent hadrons of beam and target 
particles, in the hadron restframe within a sphere of the hadron radius. 
= : uniform distribution. 

= 1 : exponential distribution with width given by PARV(34). 
— 2 : Gaussian distribution with width given by PARV(35). 

MSTV(35) : (D=2) Lorentz-boosted spatial distribution of initial partons. 

= : spatial positions of all partons (valence, sea, glue) inside parent hadron are Lorentz-boosted 
with the parent hadron's (nucleon's) 7-factor. 

= 1 : only valence quarks are boosted with the parent hadron's (nucleon's) 7-factor; scaquarks and 
gluons are smeared out symmetrically(!) in hadron (nucleon) direction of motion around the valence 
quark disc. 

= 2 : as = 1, but now seaquarks and gluons are distributed all behind(!) the valence quark disc. 
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MSTV(36) : (D=l) choice of primordial A;_L-distribution of initial partons. 
= : no primordial k±. 

= 1 : exponential fc^-distribution with width specified by PARV(37). 
= 2 : Gaussian A;^ -distribution with width specified by PARV(38). 

MSTV(37) : (D=l) choice of masses of initial partons. 
= : on mass shell with quark current masses. 

= 1 : off mass shell with space-like virtuality as determined by the excess of the partons' momentum 

over their energies. 

MSTV(38) : (D=0) switch for including effects of parton shadowing in nuclear parton structure functions 
by multiplication with effective shadowing factor RaIx) = fA{x)/fN{x) as a simple parametrization of 
EMC data, where /a and /at are the measured nuclear and nucleon structure functions, respectively. 
= : parton shadowing switched off. 
= 1 : switched on. 

MSTV(39) : (D=0) option for asymmetric (in mass) collision systems involving nuclei, h + A with 
h = ir,p,n, . . ., OT A + B with Aj^B,to choose instead of global cm. frame or fixed target frame, the 
frame of equal (hadron)nuclcon-nuclcon velocity, that is the frame in which hadron (nucleons) from 
beam side have equal but opposite directed momentum than those from target side. 

— : switched off. 

— 1 : switched on. 



PARV(21) : (D=l. GeV"-*^) the minimum resolvable partonic Bjorken-cc, denoted Xc, is set by the 
'confinement radius' Rc =PARV(21), which defines the maximum possible longitudinal spread of a 
confined parton within Az = 1/xP < l/{xcP) = Rc, where P is the longitudinal momentum of the 
mother hadron. 

PARV(22) : (D=0.95) the minimum summed x-value {J^iXi)^^^, when 'filling' a hadron with partons of 
energy (momentum) fractions Xi of the hadron. The value depends, slightly on the value of Qo (c.f. 

PARV(23)-PARV(25)), or if cuts on x are made. It should be adjusted such that when adding one 
more parton, the average {J^i Xi) should be = 1. If the sum is smaller than PARV(22), another parton 
is added to the hadron, i.e. PARV(22) = 1 - (Ei xi). 

PARV(23) : (D=l. GeV) minimum value Qo that is used in initial parton distributions (structure 
functions or alternative distributions). I.e., \i Q < Qq, then the value PARV(23) is used. 

PARV(24) , PARV(25) : (D=2.50,0.25) possibility of using effective, energy-dependent Qin value in 
structure functions if MSTV(25)=2, according to the function Qo = max[Qo) | • {Eh/GeV)^\, where 
Qo =PARV(23), a =PARV(24), h = PARV(25), and Eh = y^/hadron the average energy per hadron 
(nucleon) . 

PARV(26) : (D=1.18 fm) nuclear radius parameter ro, with R{A) = roA^/3. 

PARV(27) : (D=3.0 fm) parameter R^ in 'smeared-edge' distribution of nucleons' positions in nucleus 

with A>12 (for MSTV(33) = 1). 

PARV(28) : (D=0.545 fm) parameter Ra in Fermi distribution of nucleons' positions in nucleus with 
^ > 12 (for MSTV(33) = 2). 

PARV(31) : (D=0.25 fm) initial separation along the collision (2;-)axis of beam and target particle in the 
overall CM frame for MSTV(31)> 1. 

PARV(32) : (D=0. fm) determines value of minimum impact parameter bmim transverse to the colli- 
sion (2;-)axis, such that bmin =PARV(32) for MSTV(32) = or 1, or bmin =PARV(32) Vdnd/Tr for 
MSTV(32)=2. 
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PARV(33) : (D=l. fm) determines value of maximum impact parameter bmax, transverse to the colli- 
sion (z-)axis, such that bmax =PARV(33) for MSTV(32) = or 1, or bmax =PARV{3S)■^/and/^^ for 
MSTV(32)=2. 

PARV(34) : (D=1.19 GeV^^) width w in exponential distribution exp[—r/w] of partons' positions around 
its parent-hadron center-of-mass (c.f. MSTV(33)). (Note: w = l/v, where v = 0.84 GeV is the 
measured nucleon form factor.) 

PARV(35) : (D=2.38 GeV~^) sigma s in Gaussian distribution exp[— r^/(2 s^)] of partons' positions 
around its parent-hadron centcr-of-mass (c.f. MSTV(33)). (Note: s = 2/u, where v = 0.84 GeV is 
measured nucleon form factor.) 

PARV(36) : (D=2.) cut-off for sampling initial partons' spatial distributions from exponential or Gaus- 
sian above, such that r < aRh, where a =PARV(36) and Rh is the radius of the mother hadron. 

PARV(37) : (D=0.44 GeV) width w in exponential distribution exp[—k±/w] of primordial k± distribution 
of initial partons (c.f. MSTV(36)). 

PARV(38) : (D=0.44 GeV) sigma s in Gaussian distribution exp{—kj_/{2s'^)] of primordial k± distribu- 
tion of initial partons (c.f. MSTV(36)). 

PARV(39) : (D=2. GeV) upper cut-off for sampling primordial k± from above exponential or Gaussian 
distribution. 



3. Parton scatterings: 



MSTV(40) : (D=3) master switch for parton collisions. 
= : generally switched off. 

= 1 : only primary-primary collisions, i.e. those involving 2 partons that both have not interacted 
before. 

= 2 : as = 1 plus primary-secondary collisions, i.e. those involving 2 partons, of which one has had 
already at least one previous interaction. 

= 3 : as = 2 plus secondary-secondary interactions, i.e. all multiple collisions included. 

MSTV(41) : (D=l) choice of (semi)hard QCD collision cross-sections. 

= 1 : standard 0{a'j.) perturbative cross-sections for 2^2 and 2—^1 processes. 
= 2 : alternative of 2 ^ n cross-sections (not implemented yet!). 

MSTV(42) : (D=l) choice of po to separate hard and soft parton collisions, and to act as a regulator for 
the parton cross-sections da/dQ"^, that are divergent for vanishing momentum transfer (or invariant 
mass) Q. 

^ : fixed at po = ina^[PARV{4:2),ina^{VKIN{3),VKIN{5)]. 

— 1 : energy and mass dependent, with po = ^0(^)^)5 according to the parametrization given by 
PARV(42), PARV(43). 

= 2 : initially, before the first event, set to po{s,A), c.f. =1, but in subsequent events set equal 
to the average momentum transfer of primary parton-parton collisions po = {Qprim} accumulated by 
statistics. 

MSTV(43) : (D=2) Q^-definition in 2 — >^ 2 parton collision processes; for s-channel or fusion processes, 
is always chosen to be m? = s. Note that the parton-parton collisions invariants are denoted as 
s = (pi +p2f, i = {pi -p'lf, u = {pi -pl^f. 
= 1 : = 2siu/{s^ + t^ + u^). 
= 2 : Q2 = o.5(m^i +m,^2)■ 
= 3 : = min(— £, —u). 
= 4: Q^ = s. 
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MSTV(44) : (D=2) choice of scheme for parton colhsion-time estimate. 
= : zero cohision-time, i.e. instantaneous process. 

= 1 : finite coUision-time sampled from an exponential distribution exp(— x), where x = t/r, t the 
time in the lab, and r the mean life-time of the parton in the lab, see PARV(44). 
= 2 : finite collision-time sampled from Gyulassy-Wang distribution (2/r) x/{l + x'^)'^ , where x = t/r, 
t the time in the lab, and r the mean life-time of the parton in the lab, see PARV(44). 

MSTV(45) : (D=l) selection of maximum allowed distance between two partons in their cm-frame in 
order to be able to collide (used for ruling out widely separated pairs to speed up simulation). 
= : fixed at a value given by ri2 max = PARV(45). 
= 1 : variable at a value ri2 max = 

MSTV(46) : (D=0) choice of probability distribution W(x) from which collision probability of 2 partons 

is sampled, where x = /r^\ with r[^^^ the transverse 2-parton separation at closest approach 

and the transverse radius of interaction as given by the parton-parton cross-section a. 

= : flat distribution, W{x) = 9{a- x), with a = PARV(46). 

= 1 : exponential distribution, W{x) = exp(— x/o), a = PARV(46). 

MSTV(47) : (D=2) special treatment of 'flavor excitation' by parton collisions involving one or two 
primary hadronic or nuclear (anti)quarks by requiring a minimum momentum transfer of the resolving 
parton with the heavy quark. 

= : no special treatment of heavy quarks over light quarks. 

= 1 : requirement of minimum momentum transfer for liberation of heavy quark out of initial parton 
distribution by scattering, where the minimum required momentum trasnfer is given by PARV(47) 
times the mass of the heavy quark. 

= 2 : as = 1, but now both the struck heavy quark as well as its initial antiquark sibling are liberated 
(i.e. get on mass shell). 

MSTV(48) : (D=0) switch for including soft parton collisions according to a phenomenological cross- 
section da/dQ"^ (X a^(po)/(Q^ -|- /x^) for soft parton collisions (c.f. eq. (|5[)), where po = PARV(42), 
/J, = PARV(48), and = for t, u channel and = s for s channel. 
= : switched off. 

= 1 : soft collisions are treated complementary to hard collisions, i.e. occur only if Q < po- 

= 2 : soft collisions are treated supplementary to hard collisions, i.e. compete with the latter, 

whichever cross-section dominates. 

MSTV(49) : (D=l) switch to allow for 2^1 parton fusion processes a + b ^ c* in competition with 
2 — > 2 scattering processes a + b ^ c* d + e, such that for a given pair a,b either of those processes 
is selected according to the relative probability given by the respective cross-sections and depending 
on the 'life-time' of the particle c* . 
= : fusion processes switched off. 

= 1 : switched on, with relative weight determined by PARV(49). 



PARV(41) : (D=0.25 GeV) nominal A-value used in running Og for parton-parton collisions. 

PARV(42) : (D=0.0 GeV) cut-off for perturbative QCD cross-sections that are divergent when the mo- 
mentum transfer Q — > 0. Here the 2-body collsion scale Q is taken equal to p±, the transverse 
momentum exchanged in a parton scattering, or equal to s for annihilation or fusion processes. For 
MSTV(42)=0, the value po = a with a = PARV(42) defines hard {Q > po) and soft {Q < po) colli- 
sions. For MSTV(42)=1, see PARV(43). Note: If PARV(42) < 0, then it is set automatically upon 
initialization to = 2.0. 

PARV(43) : (D=0.0) parameter for effective po value, if MSTV(42)=1, instead of pq = const for 
MSTV(42)=0. The effective po assumes a dependence on the total collision energy and the mass 
of the (nuclear) collision system. It is parametrized as pQ{^/s, A, B) = | • (Eh/GeV)^ , where a = 
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PARV(42), b = PARV(43), and Eh = 2^/s/{A + B) with A (B) the mass number of beam (target) 
particle, if it is a hadron or nucleus. Note: If PARV(43) < 0, then it is set automatically upon 
initialization to = 0.27 - 0.29. 

PARV(44) : (D=0.1) if MSTV(44) = 1 or MSTV(44)=2, the proportionality between mean collision time 
T of a 2-parton collision process, and the associated transverse momentum generated (for s-channel 
processes, the invariant mass of the pair). It is parametrized as r = c{l/Q), with c = PARV(45) and 
Q =p± OT s. 

PARV(45) : (D=5. GeV~^) maximum separation of two partons to be able to interact through a 2-body 
colhsion. For MSTV(45)=0, it is fixed at = PARV(45), whereas for MSTV(45)=1, it is taken 
s-dependent as = PARV(45)/(\/I/GeF). 

PARV(46) : (D=1.0) parameter in the collision-probability distribution W{x) = 9{a — x) (for 

MSTV(46)=0) and W{x) = exp(-a;/a) for MSTV(46}=1, with x = rja'^V?-!? > ^12'^^ the relative 

transverse parton separation of closest approach and the transverse radius of interaction as given 
by the parton-parton cross-section a. 

PARV(47) : (D=1.0) for 'flavor excitation' processes where a primary hadronic or nuclear (anti)quark is 
struck out by a collision, an effective minimum resolution scale is required, depending on the flavor, 
such that the momentum transfer > a ■ nif, with a = PARV(47) and nif the flavor-dependent 
quark mass. This parameter acts in addition to the constraint set by PARV(42) above. 

PARV(48) : (D=0.5 GeV) parameter in the phenomenological cross-section da/dQ^ oc al{pQ)/{Q^ + iJ?) 
for soft parton collisions, where jj, = PARV(48), pq is given by PARV(42), and = p]_ for t, u channel 
and = s for s channel. 

PARV(49) : (D=10.) parametric strength of 2 ^ 1 parton fusion processes a + b ^ c competing with 
2 — > 2 scattering processes a + b ^ c + d. It is given by the ratio of the two processes which is 
proportional to l/iR^s), where l/i?^ =PARV(49). 



4. Space-like parton branchings: 



MSTV(50) : (D=l) master switch for space-like branchings. 
= : switched off 
— 1 : on. 

MSTV(51) : (D=3) level of coherence imposed on the space-like parton shower evolution. 
= 1 : none, i.e. neither k'^ values nor angles need be ordered. 

= 2 : values at branches are strictly ordered, increasing towards the hard interaction. 

= 3 : A;^ values and opening angles of emitted (on-mass-shell or time-like) partons are both strictly 

ordered, increasing towards the hard interaction. 

MSTV(52) : (D=2) structure of associated time-like parton evolution, i.e. showers initiated by emission 
off the incoming space-like partons. 

= : no associated showers are allowed, i.e. emitted partons are put on mass-shell. 

= 1 : a shower may evolve, with maximum allowed time-like virtuality set by the phase space only. 

= 2 : a shower may evolve, with maximum allowed time-like virtuality set by phase space or by the 

scale Q^, the virtuality of the space-like parton created at the same vertex, whichever is the stronger 

constraint. 



PARV(51) : (D=0.25 GeV) A-value used in ag and in structure functions for space-like parton evolution. 
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PARV(52) : (D=1.0 GeV) invariant mass cutoff qo of space-like parton showers, below which par- 
ton branchings rc not evolved. For consistency it is taken as qo = max(Qo)PA.RV(52)) where 
Qo =PARV(23) is the initial resolution scale in the hadron structure functions. 

PARV(53) : (D=4.) Q'^-scale of the hard scattering is multiplied by PARV(53) to define the maximum 
parton virtuality allowed in space-like branchings associated with the hard interaction. This does not 
apply to s-channel resonances, where the maximum virtuality is set by = s. 



5. Time-like parton branchings: 



MSTV(60) : (D=2) master switch for time-like branchings. 
= : no branchings at all, i.e. switched off. 

— 1 : QCD type branchings of quarks and gluons. 

— 2 : also emission of photons off quarks; the photons are assumed on mass-shell. 

MSTV(61) : (D=2) branching mode for time-like parton evolution. 
= 1 : conventional branching, i.e. without angular ordering. 
= 2 : coherent branching, i.e. with angular ordering. 

MSTV(62) : not used. 

MSTV(63) : (D=l) choice of formation-time scheme for parton emission. 
= : no formation time, i.e. instantanous emission. 

— 1 : finite formation time sampled from an exponential distribution exp(— x), where x = t/r, t the 
time in the lab, and r the mean life-time of the parton in the lab, see PARV(63). 

= 2 : finite formation time sampled from Gyulassy-Wang distribution (2/t) x/{1+x'^)'^, where 
X = t/T, t the time in the lab, and r the mean life-time of the parton in the lab, see PARV(63). 

MSTV(64) : not used. 

MSTV(65) : (D=0) selection of kinematics reconstruction for branchings initiated by a single off-shell 
parton. 

= : conservation of energy and jet direction, but longitudinal momentum is not conserved. 

= 1 : with full energy-momentum conservation at each branching, but the jet direction is not conserved 

due to recoil. 



PARV(61) : (D=0.29 GeV) A-value used in running for time-like parton showers. 

PARV(62) : (D=1.0 GeV) invariant mass cutoff of time-like parton evolution, below which partons are 
not assumed to radiate. Exception: for the case MSTV(84)=3, for reasons of consistency, the value is 
automatically scaled to 1.5 x PARV(62). 

PARV(63) : (D=0.01) if MSTV(63)=1 or MSTV(63)=2, the proportionality between mean-life time r of 
an off-shell parton and its virtuality A;^, is parametrized as r = c(l/\/fc^), with c = PARV(63). (Is 

inversely related to PARV(65)). 

PARV(64) : (D=4.) Q^-scale of the hard scattering is multiplied by PARV(64) to define the maximum 
parton virtuality allowed in time-like branchings associated with the hard interaction. This does not 
apply to l-channel resonances, where the maximum virtuality is set by = s. 

PARV(65) : (D=4.) invariant virtuality of off-shell partons is multiplied PARV(65) to define the 
maximum parton virtuality allowed in time-like branchings of single parton decays. (Is inversely 
related to PARV(63)). 
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6. Parton-cluster formation and cluster-hadron decay: 



MSTV(80) : (D=l) master switch for hadronization of parton clusters. 
= : off. 

— 1 : on. 

MSTV(81) : (D=l) choice of definition of 2-parton spatial separation ri2 = ri — r2 between two partons 
1, 2, for measuring probability of pre-hadronic cluster formation of the two partons. 
= : distance ri2 measured in the global frame of the collision event. 
= 1 : distance ri2 measured in the center-of-mass frame of 1 and 2. 

— 2 : distance ri2 measured in the rest-frame of parton 1. 

MSTV(82) : (D=4) choice of definition of 2-parton momentum separation pu of two partons 1, 2, for 

prc-hadronic cluster search. 

= 1 : three-momentum, i.e. pi2 = ^pl 12 + pI 12 + pI 12- 
= 2 : transverse momentum, i.e. pi2 = \JPx 12 Py 12- 
= 3 : 'pseudo-mass', i.e. pi2 = E1E2 (1 — P1P2) jimxiri-}). 
= 4 : invariant mass', i.e. p\2 = \J -£^12 ~ 12 + Py 12 + 12- 

MSTV(83) : (D=l) choice of probability distribution from which cluster formation of 2 partons is sam- 
pled. Simulates the conversion process as a 'tunneling' of the partons through a potential barrier set 
by PARV(81) and PARV(82), separating perturbative and non-perturbative regimes. 
= : fiat distribution, i.e. ^-function. 
= 1 : exponential distribution. 

MSTV(84) : (D=3) level of including color correlations among pairs of partons in the process of cluster 
formation. 

— : none, except that clustering of two partons originating from the same mother is vetoed. 

= 1 : probability of two cluster candidates being in a color singlet state is sampled uniformly with 
factors Cgg = 1/9, Cqq = 1/3, Cgq = 1/3. 

— 2 : exact matching of colors and anticolors of two partons is required to give color singlet (s). 

= 3 : general case, where arbitrary color configuration of two partons is allowed, and color singlet 
formation is balanced by additional gluon emission(s) to conserve color locally at each vertex. 

MSTV(85) : (D=0) option allowing for diquark clustering, i.e. sequential clustering of two (anti)quarks 
plus a third (anti)quark (not yet working properly!). 
= : off. 
= 1 : on. 

MSTV(87) : (D=0) choice of effective invariant masses of partons to ensure that kinematical thresholds 
for cluster-hadron decay can be satisfied. 

= : 'curent' masses are used, i.e. m/ = (0.05, 0.01, 0.01, 0.2, 1.5, 5.0, 150) GeV for 

f = g,d, u, s, c, b, t, respectively. 

= 1 : 'constituent' masses are used, i.e. m/ = (0.5, 0.32, 0.32, 0.5, 1.8, 5.2, 150) GeV for 
f = g,d, u, s, c, b, t, respectively. 

MSTV(88) : (D=0) choice of measure for rough distinction between 'slow' and 'fast' clusters in statistics 
accumulation of exogamous cluster production with numerical value of boundary given by PARV(88) . 
= : energy Ed- 
= 1 : 3-momentum Pd 

= 2 : cluster velocity (3 = Pd/Ed- 
= 3 : Lorentz factor 7 = Ed/Md- 

MSTV(89) : (D=l) smearing of cluster rapidities within an interval \y\ < yc, where yc is determined by 
PARV(89). 
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MSTV(90) : (D=l) master switch for hadronization of beam/target remnants. 



= : off. 
= 1 : on. 



MSTV(91) : (D=l) switch for the decay of unstable hadrons, emerging either directly from cluster- 
decays, or from previous unstable hadron decays, 
= : off. 

> 1 : the decay chain is iterated MSTV(91) times, but at most 5 times. 

MSTV(92) : (D=l) switch for simulating stopping effect in nucleus-nucleus collisions of initial-state nu- 
cleons by redistributing them in transverse momentum p± according to a 1/(1 -|- distribution 
(approximately an exponential for p± < 2 GeV), and in longitudinal momentum p^ (or rapidity) ac- 
cording to a Gaussian distribution exp[— Pj^g^^/(2c^)]. The parameters b and c are given by PARV(97) 
and PARV(98), respectively. Option MSTV(92) is inactive for collisions without nuclei. 
= : off, i.e. beam nucleon distributions shaped analogous to proton-proton collisions. 
— 1 : beam protons and neutrons are redistributed. 

= 2 : as = 1, plus A's associated with beam nucleon interactions are included. 
= 3 : as = 2, plus tt's and i^'s from beam nucleon interactions are accounted. 

MSTV(93) : (D=0) option to distingiush - in beam fragmentation of colisions with nuclei - into 'un- 
wounded' and 'wounded' beam nucleons. Former are those whose initial parton content has not 
altered throughout a collisions, latter are those who have been struck at least once in a hard parton 
collision. Option MSTV(93) is inactive for collisions without nuclei. 

= : off, i.e. spectator partons are not reassembled to their initial mother nucleons, but all those 

non-interacted spectator partons that belong to an 'unwoundsed' nucleon, are lumped together with 
the spectator partons from 'wounded' nucleons and hadronize as usual via the beam-fragmentation 
scheme. 

= 1 : on, i.e. for all those nucleons that have not lost a daughter parton through parton collisions, 
their initial parton content is recollected, and these 'unwounded' nucleons appear in the final state 
with their initial energy-momentum, modulo some aquired Fermi motion. In this case, only the par- 
tons from 'wounded' nucleons participate in the beam-fragmentation, specifically, only a fraction of 
all initial-state nucleons which is set by PARV(99) (see below). 

= 2 : as = 1, but now the reassembled nucleons undergo subsequent fragmentation and are broken 
up. 

MSTV(94) : (D=l) Switch to allow to cluster together excess pions from low-mass parton recombinations, 

if their density gets above a critical value, and to convert them either into energy transferred to the 
environment, or to convert then into particle production (see MSTV(95) and MSTV(96)). Specifically, 
close-by pions lumped together in clusters of mass > 1 — 5 GeV and the energy-momentum of these 
clusters is then transferred to the remaining particle system, or, alternatively, the energy-momentum 
of TT clusters with mass > 1 — 5 GeV is converted into pair production of 7r^,p, a;. 
= : off, i.e. no action taken. 



MSTV(95) : (D=0) For MSTV(94) > 0, the critical density above which pions may percolate to clusters 
is set by MSTV(95). IF MSTV(95) < 0, its value is set automatically upon initialization. 

MSTV(96) : (D=0) For MSTV(94) > 0, the maximum mass of clusters of pions that can be percolated 
is set by MSTV(96). IF MSTV(96) < 0, its value is set automatically upon initialization. 

MSTV(97) : (D=0) Switch to select shape of p± spectrum of final-state hadrons emerging from fragmen- 
tation of beam/target clusters in collisions involving nuclei. 



MSTV(98) : (D=0) Switch to select shape of y spectrum of final-state hadrons emerging from fragmen- 
tation of beam/target clusters in collisions involving nuclei. (Notice that in case MSTV(40) > 0, the 



— 1 : on. 



= : exponential distribution oc exp(— 6mj^], m± 
= 1 : power-law distribution oc l/(m^)^, m± = a 
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camel-like option is selected by default, in order to compensate the spectra of hadrons from partonic 

clusters and those from soft beam/target clusters.) 

= : flat-top distribution, i.e. a central plateau with Gaussian tails. 

= 1 : camel-like distribution, i.e. a central dip with Gaussian on either side. 

MSTV(99) : (D=0) Switch to adjust left-right imbalance of beam/target hadron fragments, if there is 
such occurrence. 
= : off. 
= 1 : on. 



PARV(81) : (D=3.6 GeV ^) minimum separation of 2 partons, below which a cluster cannot be formed 
(denoted in Ref. [^). 

PARV(82) : (D=4. GeV~^) critical separation of 2 partons that sets the space-time scale for cluster 
formation (denoted Rcrt in Ref. Is related to the average cluster size and typically about 10 % 

above R^. 

PARV(83) : (D=1.0) 'slope' of the exponential probability distribution W{R) = e{R - 
R^) exp{-FARY {83) R/ Rcrt) for MSTV(83)> 0, describing the probability of 2 partons to form a 
pre-hadronic cluster by tunneling through the 'potential barrier' marked by R^ and Rcrt, with average 
cluster size Rd = 2i?crt/PARV(83). 

PARV(84) : (D=0.3 GeV) minimum invariant mass of 2 partons required to form a cluster being neces- 
sarily above threshold for hadron decay. 

PARV(85) : (D=1000. GeV) maximum allowed invariant mass of 2 partons to form a cluster with mass 
Mci , in order to provide option to suppress production of heavy clusters - in addition to the requirement 
of PARV(81),PARV(82). Note: default corresponds to no suppression. 

PARV(86) : (D=5. GeV~^) maximum allowed separation of 2 partons that are unambigous cluster 
candidates. Cluster formation is enforced above this value, in order to avoid unphysical separation. 
In this case the requirement Md < PARV(85) is overridden. 

PARV(87) : (D=l.) 'slope' of exponential Hagedorn-type distribution for sampling the masses of clusters 
in their 2-body decays, given by exp(— a • Mci/Mq), where a = PARV(87), Md is the cluster mass and 
Mq = 0.2 GeV is the 'Hagedorn temperature' (see Ref. [Q). 

PARV(88) : (D=l. GeV) boundary for distinguishing 'slow' and 'fast' clusters in statistics accumulation 
of exogamous cluster production (c.f. MSTV(88)). 

PARV(89) : (D=0.2) determines value of yc = PARV(89) 

Vbeam Smearing of parton cluster rapidities 

within an interval \y\ < if MSTV(89) > 0. 
PARV(90) : (D=0.0) determines value of yc = PARV(90) 

Ubeam for Smearing of beam cluster rapidities 
within an interval \y\ < yc- If PARV(90) < 0, it is ste automatically to pre-fitted value, otherwise, 
PARV(90) > as user input takes this as fixed value. 

PARV(91) : (D=0.0) Parametric factor for charged multiplicity from soft fragmentation of beam/target 
remnants that governs the magnitude of Nchi^/s) of hadrons emerging from spectator partons. If 
PARV(91) < 0, then the value of PARV(91) is set automatically to pre-fitted value, otherwise, 
PARV(91) > as user input takes this as fixed value. (Note: applies only if if MSTV(90) = 1). 

PARV(92) : (D=0.0) Factor multiplying the energy of the remnant beam/target clusters that are formed 
from the spectator partons, and that are fragmented with the total energy E=PARV(92)- (E'beam + 
Etarget)- If PARV(92) < 0, then it is set automatically upon initialization. 
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PARV(93) : (D=1.0) Nuclear dependence in eA, hA or AB collisions of the charged multiplicity N^h, 
resulting from soft fragmentation of nuclear beam or/and target remnants of the original nuclei with 
mass number A (B) (case of hadron h corresponds to A = 1 or i? = 1). It is parametrized as: 

N^.t'^HV^^ B) = [{A + B)/2f Nf\^s), where (3 = PARV(93). 

PARV(94) : (D=15.GeV~^) Time in the global collision frame when soft fragmentation starts to become 
gradually active (default corresponds to 35 fm.) 

PARV(95) : (D=1000.GeV~^) The positions of particles produced in fragmentation of both parton clus- 
ters and beam/target clusters, are assigned such they become active only after an effective formation 
time Ai = PARV{^h) x min(£;/M2,PARV(95)) (default value corresponds At < 200 fm). 

PARV(96) : (D=40.) proportionality factor in the uncertainty relation At = PARV{9Q) x E/M'^, which 
determines the mean formation time of particles produced in soft beam/target fragmentation (c.f. 
PARV(95)). 

PARV(97) : (D=0.35 GeV) parameter 6 in 1/(1 + hp\)* distribution for simulating stopping effect in 
nucleus-nucleus collisions of protons and neutrons (and optionally associated A's) by redistributing 
them in transverse momentum p_\_ (c.f. MSTV(92)). 

PARV(98) : (D=0.14 GeV) parameter c in Gaussian distribution exp[— Pj^g^^/(2c^)] for simulating stop- 
ping effect in nucleus- nucleus collisions of protons and neutrons (and associated A's) by redistributing 
them in longitudinal momentum pz (c.f. MSTV(92)). 

PARV(99) : (D=0.) Of all initial-state nucleons which have not lost a daughter parton through parton 
collisions (those that remained 'unwounded' nucleons) and appear in the final state with their initial 

energy-momentum, a fraction <PARV(93)< 1 is fragmented via the usual beam/target fragmentation 
scheme, whereas the fraction 1— PARV(99) appears as original nucleons in the final state. Is inactive 
for MSTV(93) = 0. 



7. 'Afterburner' Hadron cascade: 



MSTV(IOO) : (D=3) switch to select various levels of final-state hadron cascading in terms of scattering 
and resonance production with decay of direction of primary hadrons emerging from both, the parton 
cluster formation and the beam/target fragmentation. (See additional control switches MSTV(lOl)- 
MSTV(103)). 

= : switched off completely, i.e., no final-state hadron cascade. 

= 1 : switched on, but including only hadron cascading of hadrons from parton cluster decays and 
excluding re-interactions among hadrons from the beam/target fragmentation. 

= 2 : switched on, but in contrast to = 1, now only hadron cascading of hadrons from beam/target 
fragmentation are allowed ans re-interactions among hadrons from the parton cluster decays are ex- 
cluded. 

= 3 : switched on in full, i.e., any hadron can re- interact irrespective to its origin. 

MSTV(lOl) : (D=l) Additional selection option in addition to MSTV(IOO) that allows to include/exclude 
hadrons that have not become 'real', i.e., that have not passed through their formation time. 
= : only on-shell hadrons can interact, that is, those which have passed their formation time. 
= 1 : both on-shell hadrons and off-shell hadrons which have not passed their formation time yet, 
can take part in the hadron cascading. 

MSTV(102) : (D=3) switch to include/exclude interactions in which pre-hadronic clusters (not yet 
formed hadrons) are involved (note: applies only if MSTV(lOl) = 1). 
= : no pre-hadronic clusters can interact. 
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= 1 : only pre-hadronic clusters from parton cluster formation may participate in the hadron cacade. 
= 2 : in contrast to = 1, now only clusters from beam/target fragmentation are included in the 
hadron cacade. 

= 3 : any pre-hadronic cluster can interact with other hadrons or clusters, irrespective of the origin. 

MSTV(103) : (D=2) Sets (for MSTV(102) > 1) the maximum number of preformed hadrons for which 
is searched for in the event record in order to reconstruct the original cluster from whose decay the 
hadrons emerged. 

MSTV(105) : - not used - 



PARV(IOO) - PARV(104) : - not used - 

PARV(105) : (D=5. GeV^^) maximum separation of two hadronic particles to be able to interact through 
a 2-body collision. For MSTV(105)=0, it is fixed at = PARV(105), whereas for MSTV(105)=1, it 
is taken s-dependent as ru = PARV(105)/(-/s/GeF). 

PARV(106) : (D=1.0) parameter in the collision-probability distribution W{x) = 9{a — x) (for 
MSTV(106)=0) and W{x) = exp(-a;/a) for MSTV(106)=1, with x = rif^V^i?, ^'if^^ the rela- 
tive transverse hadron separation of closest approach and the transverse radius of interaction as 
given by the hadron-hadron cross-section a. 



8. Other settings: 



• Handling of input/output: 

MSTV(llO) : (D=l) direction of program output. 

= : all output is directed to standard output unit 6. 
= 1 : output is directed to units MSTV(111)-MSTV(113). 

MSTV(lll) : (D=10) unit number to which general output is directed. 

MSTV(112) : (D=6) unit number for writing 'on-line' information. 

MSTV(113) : (D=20) unit number for listing warnings and error messages. 

MSTV(114) : (D=0) switch for general information on unit MSTV(lll), which is directed automatically 
to file VNIRUN.DAT. Allows to select the amount of output provided on the global performance of a 
simulation consisting of of a sample of collision events: 

= : only minimal output concerning selected process and main parameters is written out, as well 

as final summary of results. 

= 1 : as =0, plus a listing the 1st event in momentum and position space. 

— 2 : as =1, plus listing of number and properties of elementary subprocesses that occurred during the 
simulation, written to files VNICOL.DAT, VNISBR.DAT and VNITBR.DAT for collisions, space-like 
and time-like branchings, respectively. 

= 3 : as =2, plus detailed event listing of particle record at the end of each single event, written 
to VNIRUN.DAT. NOTE: for large particle systems, this slows down the program substantially, and 
moreover requires a lot of disc space if a large number of events are printed). 

MSTV(115) : (D=2) switch for 'on-line' information on unit MSTV(112), i.e. writing of initialization 
and termination info, and initial and final entries for each event as the simulation goes on. 
= : no 'onhne output'. 

= 1 : only initialization and finalization info is printed. 

= 2 : 1-line info at beginning and at the end of each event to keep control of the program performance. 
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= 3 : as = 2, plus listing of entries in particle record after each single time-step in an event (useful 
for AB collisions). 

= 4 : as = 3, plus writing 'standardized' output to file VNIOSC.DAT, according to the OSCAR 
(Open Standard Codes At RHIC) conventions established at RIKEN-BNL, June 1997. 

• Miscellanous control switches for subroutines: 

MSTV(120) : (D=0) flag to indicate initialization of certain HERWIG common blocks and default values 
that are necessary for using parts of the HERWIG program for the hadronization of clusters. Is set 
equal to 1 after first intitialization call. 

MSTV(121) : (D=0) if MSTV(121) is set to 1 before a VNIROBO caU, the V vectors are reset (in the 
particle range to be rotated/boosted), i.e., they are set to before the rotation/boost. If MSTV(121) 
is equal to zero, the V vectors are not reset. MSTV(121) is inactive during a VNIROBO call and is 
set back to upon return. 

MSTV(122) : (D=0) specifies in a VNIROBO call the type of rotation of the P, R, and V vectors. 
The rotation is clockwise (active rotation) for MSTV(122)=0 and anticlockwise (passive rotation) 
otherwise. The value of MSTV(122) is set back to upon return. 

MSTV(123) : (D=0) specifies in a VNIROBO call for performing a boost, whether the vectors R are 
boosted (MSTV(123)=0) or not (MSTV(123)=1). The vectors P are always boosted. The value of 
MSTV(123) is set back to upon return. 

MSTV(124) : (D=l) pointer to lowest entry in particle record VNIREC1/VNIREC2 to be included in 
data analysis routines VNIANA1-VNIANA4. 

MSTV(125) : (D=100000) pointer to highest entry in particle record VNIREC1/VNIREC2 to be included 
in data analysis routines VNIANA1-VNIANA4. 

MSTV(126) : (D=20) number of lines in the beginning of the particle record that are reserved for internal 
event history information. The value should not be reduced, but can be increased if necessary. 

MSTV(127) : (D=l) in hsting of the current state of the particle record by calling the subroutine 
VNILIST, it can be chosen between listing either color /anticolor C A and the space time origin of a 
particle / according to the information contained V{I,5) explained above, ( = ), or, alternatively 
the mother KMO and color flow information contained in K{I, 3) — K{I, 5), ( = 1 ). All other listed 
quantities are the same for either option. 

MSTV(128) : (D=0) specifies in calls to VNIEDIT and VNIROBO the classes of particles to be included. 

The default = eax^ludes all inactive or decayed particles with K{1, 1) < or K{1, 1) > 10, whereas 
the setting = 1 iricludes any entry listed in the particle record. MSTV(128) is set back to upon 
return. 

MSTV(129) : (D=0) specifies in calls to the 'color-Hnking' routine VNIJCOL the way of encoding the 
color connections between partons in the components J = 4, 5 of the array K(I,J). For very large 
number of entries (> 35000) in the record VNIREC, MSTV(129) is set automatically to = 1 in the 
hadronization process, in order to save memory space. It is reset to = afterwards. 

• Handling of errors and warnings: 

MSTV(130) : (D=2) check on possible errors during program execution. 

= : errors do not cause any immediate action. 

= 1 : possible errors are checked and in case of problem, it is excited from the subprogram, but the 
simulation is continued from there on. For the first MSTV(lOl) errors a message is printed; after that 

no messages appear. 

= 2 : possible errors are checked and in case of problem, the the simulation is forcibly terminated. 
For the first MSTV(lOl) errors a message is printed; after that no messages appear. 
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MSTV(131) : (D=10) max number of errors that are printed. 

MSTV(132) : (D=l) printing of warning messages. 
= : no warnings are written. 

— 1 : first MSTV(103) warnings are printed, thereafter no warnings appear. 
MSTV(133) : (D=10) max number of warnings that are printed. 
• Version, date of last change: 



MSTV(140) 
MSTV(141) 
MSTV(142) 
MSTV(143) 
MSTV(144) 
MSTV(145) 



(D=l) Print-out of VNI logo on first occasion; MSTV(140) is reset to afterwards. 

VNI version number. 

VNI subversion number. 

year of last change of VNI version 

month of last change of VNI version. 

day of last change of VNI version. 



9. Statistics, event study and data analysis: 



• Particle recorc? VNIREC1/VNIREC2: Choosing event history information. 

MSTV(150) : (D=0) specifics whether the complete particle history with all particles produced at any 
time during an event, is kept throughout the simulation (rather inefficient for large systems) or whether 
it is compressed at the end of each time step (thereby loosing information on parent/child genealogy 
etc.). 

= : particle record is greatly compressed after each time step, all 'dead' particles are removed and 
genealogical history is lost. 

= 1: particle record is partly compressed, keeping some of the genelogical history. 
= 2: full particle history is kept throughout. 



• Subroutine VNIANAl 


MSTV(151) 


(D= 


=0) 


MSTV(152) 


(D= 


=0) 


MSTV(153) 


(D= 


=0) 


MSTV(154) 


(D= 


=0) 


MSTV(155) 


(D= 


=0) 


• Subroutine VNIANA2 


MSTV(161) 


(D= 


=0) 


MSTV(162) 


(D= 


=0) 


MSTV(163) 


(D= 


=0) 


• Subroutine YmAN A3 


MSTV(171) 


(D= 


=0) 
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MSTV(172) 
MSTV(173) 
MSTV(174) 
MSTV(175) 
MSTV(176) 



(D=0) Longitudinal space-time spectra dN/dz. 
(D=0) Transverse momentum spectra (l/p±) dN/dp± 
(D=0) Transverse space-time spectra {l/r±) dN/dr±. 
(D=0) Distributions of cluster sizes and masses. 
(D=0) Polarization profile of cluster density. 



• Subroutine VNIANA4: Statistics on partons and produced hadrons. 



MSTV(181) 


(D= 


=0) 


MSTV(182) 


(D= 


=0) 


MSTV(183) 


(D= 


=0) 


MSTV(184) 


(D= 


=0) 



MSTW(200), PARW(200): generated quantities and statistics 



1. General: 



MSTW(l) 
MSTW(2) 
MSTW(3) 



total number of collision events to be generated, 
total number of time steps per collision event. 



physics process A + B that is simulated. The current available beam {A) and target (B) 
particles and the collision processes are the ones listed before in Table 1, with MSTW(3)=IPR0 (see 
also selection switch MSTV(30)). 

MSTW(4) : Overall Lorentz frame of reference for event siumulation. 
= 1 : center-of-momentum frame (CMS) of beam and target particle. 
= 2 : fixed-target frame (FIXT) with target particle at rest. 

= 3 : user-defined frame (USER) with given 3-momentum of beam and target. Particles are assumed 
on the mass shell. 

= 4 : user-defined frame (FOUR) with given J^-momentum ( i.e., 3-momentum and energy) of beam 
and target particle. The particles need not to be on the mass shell. 

= 5 : user-defined frame (FIVE) with given 5-momentum ( i.e., 3-momentum, energy and mass) of 
beam and target particle. The particles need not to be on the mass shell, but 4-momentum and mass 
information must(!) match. 



MSTW(5) 
MSTW(6) 
MSTW(7) 
MSTW(8) 



KF flavour code for beam particle A. 
KF flavour code for target particle B. 

type of incoming beam particle A: 1 for lepton, 2 for hadron, and 3 for nucleus, 
type of incoming target particle B: 1 for lepton, 2 for hadron, and 3 for nucleus. 
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MSTW(9) : combination of incoming beam and target particles. 
= 1 : lepton on lepton 
= 2 : lepton on liadron 
= 3 : lepton on nucleus 
= 4 : hadron on lepton 
= 5 : liadron on liadron 
= 6 : hadron on nucleus 
= 7 : nucleus on lepton 
= 8 : nucleus on hadron 
= 9 : nucleus on nucleus 

MSTW(IO) : performance flag for current event. Is =0 at beginning, and set =1 if during the evolution 
the event turns out to be rejectable, either due to unphysical kinematics, particle combinations, etc., 
or due to numerical errors. 

MSTW(ll) : current collision event. 

MSTW(12) : current time step in this collision event. 

MSTW(13) : current number of successful collision events, i.e. those that completed gracefully with 
MSTW(10)=0. 

MSTW(14) : counter for " non-diffractive" collisions events, i.e. those that involved at least one parton 
collision (in hadron or nucleus collisions). 

MSTW(15) : Status of Lorentz transformation between different global Lorentz frames in which the 
collision event is simulated (which may be different from the initially specified frame, c.f. MSTW(4)). 
The value of MSTW(15) saves the last performed transformation: 
= 1 : from global cm-frame to fixed-target or user-specified frame; 

= 2 : from global cm- frame to hadronic center-of-mass in DIS (photon-hadron cm- frame); 
= -1 : from fixed-target or user-specified frame to global cm- frame; 

= -2 : from hadronic center-of-mass in DIS (photon-hadron cm-frame) to global cm-frame. 



PARW(l) : time increment TINC{I) in current time step /. 
PARW(2) : time TIME{I) in current time step /. 

PARW(3) : local machine time (in CPU) passed when simulation started. 
PARW(4) : local machine time (in CPU) when simulation ended. 

PARW(5) : conversion factor CPU seconds (machine-dependent: for most processors, it is equal to 
0.01, i.e. 1 sec = 100 CPU. 



PARW(ll) 
PARW(12) 
PARW(13) 



Total invariant ^/s = Eqm, i-e. the total CM energy of the collision system. 
Invariant s = E^j^^ mass-square of complete system. 



CM energy per nucleon Ecm/{^ + B) (where A (B) is the mass number of beam (target), 
and A + B is the total number of nucleons in the system) in collisions involving nuclei. Is equal to 
EcM for elementary particle- and hadronic collisions. 

PARW(14) : Invariant Sgff = E}.j^/{A + B) mass-square per nucleon in collisions involving nuclei. Is 
equal to s for elementary particle- and hadronic collisions. 

PARW(15) : mass Ma of beam particle. 

PARW(16) : mass Mb of target particle. 
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PARW(17) : longitudinal momentum P^a of beam particle in specified global frame. 

PARW(18) : longitudinal momentum Pzb of target particle in specified global frame. 

PARW(19) : angle 9 of rotation from CM frame to user-defined frame (also, in e^e~ via W^W~ , the 
rotation of the T^-pair along the z-axis). 

PARW(20) : azimuthal angle of rotation from CM frame to user-defined frame, corresponding to 
PARW(18). 



2. Initial state of collision system: 



MSTW(21) 
MSTW(22) 
MSTW(23) 
MSTW(24) 
MSTW(25) 
MSTW(26) 
MSTW(27) 
MSTW(28) 
MSTW(29) 
MSTW(30) 
MSTW(31) 
MSTW(32) 
MSTW(33) 
MSTW(34) 
MSTW(35) 
MSTW(36) 
MSTW(37) 
MSTW(38) 



number of neutrons of beam particle A. 
number of neutrons of target particle B. 
number of protons of beam particle A. 
number of protons of target particle B. 
number of d-valence quarks of beam particle A. 
number of d- valence quarks of target particle B. 
number of rt-valence quarks of beam particle A. 
number of tt- valence quarks of target particle B. 
number of s-valence quarks of beam particle A. 
number of s-valence quarks of target particle B. 

total number of initial state partons in collisions involving one or more hadron or nucleus, 
number of initial gluons. 

10^ X number of d quarks -|- number of d antiquarks. 
10^ X number of u quarks + number of u antiquarks. 
10^ X number of s quarks -|- number of s antiquarks. 
10^ X number of c quarks -|- number of c antiquarks. 
10^ X number of b quarks -|- number of b antiquarks. 
10^ X number of t quarks -|- number of i antiquarks. 



PARW(21) : absolute value of velocity Pa of beam particle A. 

PARW(22) : absolute value of velocity Pb of target particle^. 

PARW(23) : Lorentz factor ja of beam particle A. 

PARW(24) : Lorentz factor of target particle B. 

PARW(25) : rapidity Ya of beam particle A. 
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PARW(26) : rapidity Yb of target particle B. 

PARW(27) : radius Ra (in 1/GcV) of beam particle in its restframe. 

PARW(28) : radius Rb (in 1/GeV) of target particle in its restframe. 

PARW(29) : radius i?„ (in 1/GeV) of neutron within in nucleus. 

PARW(30) : radius Rp (in 1/GeV) of proton within in nucleus. 

PARW(31) : accumulated average Qq of initial hadron (nucleus) parton distributions (equals the average 
momentum scale of primary parton collisions {Qprim))- 

PARW(32) : accumulated average momentum fraction xa of beam-side initial hadron (nucleus) parton 
distributions. 

PARW(33) : accumulated average momentum fraction of target-side initial hadron (nucleus) parton 
distributions. 

PARW(34) : 3 times the total charge of beam-side particle. 
PARW(35) : 3 times the total charge of target-side particle. 
PARW(36) : ratio h/h 

max of actual impact parameter ftmjn ^ b ^ b^nax to maximum allowed impact 
parameter of beam and target particles, if variable impact parameter is chosen (c.f. MSTV(32) and 
PARV(32), PARV(33)). 

PARW(37) : azimuthal angle of colliding beam and target particles around the beam axis, if variable 
impact parameter is chosen. 

PARW(38) : sign of r^j-coordinate of beam-particle A in global frame for variable impact parameter 
selection (= — sign of target-particle B). 

PARW(39) : sign of r^-coordinate of beam-particle A in global frame for variable impact parameter 

selection (= — sign of target-particle B). 



3. Parton scatterings: 



MSTW(40) 


: Total number of 2-parton collisions a + b N {N = 1,2). 


MSTW(41) 


: Number of hard 2^2 collisions of q + q,q + q,q + q. 


MSTW(42) 


: Number of hard 2 — 2 collisions of q + g,q + g. 


MSTW(43) 


: Number of hard 2—^2 collisions of g + g. 


MSTW(44) 


: Number of soft 2^2 collisions of q + q,q + q,q + q- 


MSTW(45) 


: Number of soft 2^2 collisions of q + g,q + g. 


MSTW(46) 


: Number of soft 2 — > 2 collisions of g + g. 


MSTW(47) 


: Number of 2 — >^ 1 fusions of q + q. 


MSTW(48) 


: Number of 2 — 1 fusions of q + g,q + g. 


MSTW(49) 


: Number of 2 — 1 fusions of g + g. 
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PARW(40) : total number of 'primary' (i.e., first) collisions among primary partons. 

PARW(41) : accumulated average of the 'primary' collisions. 

PARW(42) : accumulated average ^fi of hard 2^2 collisions. 

PARW(43) : accumulated average rapidity y* = yl — J/2 of hard 2 — s- 2 collisions. 

PARW(44) : accumulated average p± of hard 2^2 collisions. 

PARW(45) : accumulated average VS of soft 2 —>■ 2 collisions. 

PARW(46) : accumulated average rapidity y* = yl — ^2 of soft 2 — 2 collisions. 

PARW(47) : accumulated average p± of soft 2^2 collisions. 

PARW(48) : accumulated average s/I of hard 2 — s- 1 collisions. 

PARW(49) : accumulated average rapidity y* = y* — yl of hard 2^1 collisions. 

PARW(50) : integrated 2-parton cross-section a{s) = J dp\{da{s,p\)/dp\). 



4. Space-like parton branchings: 



MSTW(50) : Total number of space-like branchings a^h + c. 

MSTW(51) : Number of space-like processes CMshower b + c. 

MSTW(52) : Number of space-like processes q{q) q{q) + g. 

MSTW(53) : Number of space-like processes g ^ g + g- 

MSTW(54) : Number of space-like processes g ^ q + q- 

MSTW(58) : Number of space-like processes q/q — q{q) + 7. 



PARW(50) 
PARW(51) 
PARW(52) 
PARW(53) 
PARW(54) 
PARW(55) 
PARW(56) 
PARW(57) 



average of starting scale for space-like branching evolution. 



average \J —q^ of virtuality of branching particle. 

average value oi x \n. x ^ x'x" branching. 

average fraction z = x' jx oi space-like branchings. 

average relative transverse momentum q^ of space-like branchings. 

average longitudinal momentum q^ of branching particle. 

average energy of branching particle. 

average invariant mass V^j^^ of branching particle. 



5. Time-like parton branchings: 
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MSTW(60) : Total number of time-like branchings a ^ b + c. 

MSTW(61) : Number of time-like processes C M shower h + c. 

MSTW(62) : Number of time-like processes q{q) — > q{q) + g. 

MSTW(63) : Number of time-like processes g ^ g + g- 

MSTW(64) : Number of time-like processes g ^ q + q- 

MSTW(68) : Number of time-like processes q{q) — > q{q) + 7. 



PARW(60) : average \/Q^ of starting scale for time-like branching evolution. 

PARW(61) : average ^/k^ of virtuality of branching particle. 

PARW(62) : average value of x in a; ^ x'x" branching. 

PARW(63) : average energy fraction z = x' /x of time-like branchings. 

PARW(64) : average relative transverse momentum k± of time-like branchings. 

PARW(65) : average longitudinal momentum kz of branching particle. 

PARW(66) : average energy of branching particle. 

PARW(67) : average invariant mass Vm^ of branching particle. 



6. Parton-cluster formation and cluster-hadron decay: 



MSTW(80) : Total number of 2-parton cluster formations. 

MSTW(81) : Total number of gg clusterings, gg CC + X, where X is either 'nothing', or g, or gg. 

MSTW(82) : Number of processes, gg ^ CC + g. 

MSTW(83) : Number of processes, gg ^ CC + gg. 

MSTW(84) : Total number of qq clusterings, qq — > CC + X, where X is either 'nothing', or g, or gg. 

MSTW(85) : Number of processes, qq ^ CC + g. 

MSTW(86) : Number of processes, qq CC + gg. 

MSTW(87) : Total number of gq {gq) clusterings, gq ^ C + X , where X is either q or gq. 

MSTW(88) : Number of processes, gq ^ CC + gq. 

MSTW(89) : Total number of {qq)q CC, {qq)q CC clusterings. 

MSTW(90) - MSTW(99) : Same as MSTW(80) - MSTW(89), but now for the corresponding numbers 
of 'exogamously' produced clusters, with 'exogamy index' ed = (ej -|- ej)/2 unequal to or 1 (see Ref. 
0)- 



MSTW(IOO) : Total number of 'fast' clusters, classified according to the choices of MSTV(88)] and 
PARV(88). 
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MSTW(lOl) : Total number of 'slow' clusters, the fraction complimentary to the 'fast' one. 

MSTW(llO) : Total number of primary "neutral" particles per event, produced directly by cluster-hadron 
decays. 

MSTW(lll) : Number of primary leptons and gauge bosons. 

MSTW(112) : Number of primary light mesons. 

MSTW(113) : Number of primary strange mesons. 

MSTW(114) : Number of primary charm and bottom mesons. 

MSTW(115) : Number of primary tensor mesons. 

MSTW(116) : Number of primary light baryons. 

MSTW(117) : Number of primary strange baryons. 

MSTW(118) : Number of primary charm and bottom baryons. 

MSTW(119) : Number of other particles. 

MSTW(120) - MSTW(129) : Same as MSTW(llO) - MSTW(119), but now for charged particles only. 

MSTW(130) - MSTW(139) : Same as MSTW(llO) - MSTW(119), but now for the numbers of sec- 
ondary neutral particles per event (i.e. those produced by decays of unstable particles). 

MSTW(140) - MSTW(149) : Same as MSTW(130) - MSTW(139), but now for charged particles only. 



PARW(81) : accumulated average space-time separation of clustered partons. 

PARW(82) : accumulated average energy-momentum separation of clustered partons. 

PARW(83) : minimum encountered separation of any two clustered partons with respect to chosen space- 
time measure (c.f. MSTV(81)). 

PARW(84) : maximum encountered separation of any two clustered partons with respect to chosen 
space-time measure (c.f. MSTV(81)). 

PARW(85) : minimum encountered separation of any two clustered partons with respect to chosen energy- 
momentum measure (c.f. MSTV(82)). 

PARW(86) : maximum encountered separation of any two clustered partons with respect to chosen 
energy-momentum measure (c.f. MSTV(82)). 

PARW(91) - PARW(99) : Gives the 'exo(endo)-gamy distribution' (see Ref. |jlO|) of clusters formed 
from partons with 'exogamy' index ej,ej, such that = {ei + ej)/2. Inital beam/target particles have 
e = 0(1), so that lies in the interval [0,1], which is binned as 0, 0.05, 0.2, 0.3, 0.45, 0.55, 0.7, 0.8, 0.95, 1, 
with PARW(90)-PARW(99) giving the numbers of clusters in the bins. Note: only those clusters are 
counted that are classified as 'fast' (c.f. MSTW(88)). 

PARW(lOl) - PARW(109) : Same as PARW(91)-PARW(99), but for 'slow' clusters only. The sums 
PARW(91)+PARW(101), etc., hence give the total numbers. 

PARW(llO) - PARW(149) : not used. 
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7. 'Afterburner' Hadron cascade: 



MSTW(150) : Total number of hadronic collisions involving clusters and/or hadrons. 

MSTW(151) : Number of collisions between two clusters C + C, where only light quark flavors are 
involved. 



MSTW(152) 
MSTW(153) 
MSTW(154) 
MSTW(155) 



Number of C + C collisions where strange quark flavor is involved. 
Number of C + C collisions where charm quark flavor is involved. 
Number of C + C collisions where bottom quark flavors is involved. 

- not used - 



MSTW(156) - MSTW(160) : Same as MSTW(151)-MSTW(155), but for cluster-meson collisions C + 
Ml, where Ml are light mesons only. 

MSTW(161) - MSTW(165) : Same as MSTW(156)-MSTW(160), but for cluster-meson collisions C + 
M2, where M2 are heavy or excited mesons. 

MSTW(166) - MSTW(170) : Same as MSTW(161)-MSTW(165), but for cluster-baryon collisions 
C + B, where B is any type of baryon. 

MSTW(171) - MSTW(175) : Same as MSTW(166)-MSTW(170), but for meson-meson collisions M1 + 
Ml, where Ml and M2 are light mesons. 

MSTW(176) - MSTW(180) : Same as MSTW(171)-MSTW(175), but for meson-meson collisions M1 + 
M2, where Ml is a light meson and M2 an excited or heavy meson. 

MSTW(181) - MSTW(185) : Same as MSTW(176)-MSTW(180), but for meson-baryon collisions M1 + 

B, where Ml is a light meson and B any type of baryon. 

MSTW(186) - MSTW(190) : Same as MSTW(181)-MSTW(185), but for baryon-baryon collisions 

B + B, where both B can be any type of baryon. 

MSTW(191) - MSTW(195) : Same as MSTW(186)-MSTW(190), but for baryon-meson collisions 5 + 
M2, where B is any type of baryon and M2 an excited or heavy meson. 

MSTW(196) - MSTW(200) : Same as MSTW(191)-MSTW(195), but for meson-meson collisions M2-f 
M2, where both M2 are excited or heavy mesons. 



PARW(150) 
PARW(151) 
PARW(152) 
PARW(153) 
PARW(154) 
PARW(155) 
PARW(160) 
PARW(165) 



: ( ^/s) of all cluster /hadron collisions. 

: ( ^/s) of cluster /hadron collisions involving only light flavors. 

: (\/s) of cluster /hadron collisions involving strange flavor. 

: ( of cluster/hadron collisions involving charm flavor. 

: ( \/— t ) of cluster/hadron collisions involving bottom flavor. 

- PARW(159) : - not used - 

- PARW(164) : of cluster/ hadron collisions in paralell to 

- PARW(169) : - not used - 
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PARW(170) - PARW(174) : (V^) of cluster/ hadron collisions in paralell to PARW(150)- 
PARW(154). 

PARW(175) - PARW(179) : - not used - 

PARW(180) - PARW(184) : {uprod) of particles produced in cluster/ hadron collisions in paralell to 
PARW(150)-PARW(154). 

PARW(185) - PARW(190) : - not used - 

PARW(191) : Total number of 'unwounded' nucleons at the end of current event of collision with nuclei. 

PARW(192) : Number of 'unwounded' nucleons at the end of current ewenf which are within the transverse 
area of nuclear overlap. 

PARW(193) : Number of 'unwounded' nucleons at the end of current event which are outside the trans- 
verse area of nuclear overlap. 

PARW(194) - PARW(196 : Numbers of 'unwounded' nucleons accumulated over all events in collisions 
with nuclei, in paralell to PARV(191) - PARV(193). 

PARW(197) - PARW(200) : - not used - 



3.8 Kinematics cuts and selection of subprocesses 



The commonblock VNIKIN contains the arrays VKIN and MSUB, which may be used to impose specific 
kinematics cuts and to switch on/off certain subprocesses, respectively. The notation used in the following 
is that 'hat' denotes internal variables in a partonic interaction, while '*' is for variables in the global CM 
frame of the system as-a- whole. All dimensionful quantities are in units of GeV or GeV~-^. 

Kinematics cuts can be set by the user before the VNIXRUN call. Most of the cuts are directly related 
to the kinematical variables used in the evolution of the system in each event, and in the event data analysis. 
Except for the entries VKIN(7) - VKIN(12), the assignments follow the convention of PYTHIA, in order to 
allow an easy interfacing. 

Note, that in the current version of VNI, only the values VKIN(7)-VKIN(12) are actively used. All 
other quantities are related to parton collision processes involving one or more partons, which are only of 
relevance in reactions with hadrons (nuclei) in the initial state. As stated before, the latter are presently 
worked out, but not yet included in the package. 

C0MM0N/VNIKIN/VKIN(200) , MSUB (200) , ISET (200) 

Purpose: To impose kinematics cuts on the particles' interactions and their space-time evolution. Also, to 
allow the user to run the program with a desired subset of processes. 

VKIN(l), VKIN(2) : (D=2.,-l. GeV) range of allowed m = yfs values. If VKIN(2) < 0, the upper 
limit is inactive. 

VKIN(3), VKIN(4) : (D=2.,-l. GeV) range of allowed (D=0.,-1.) range of allowed pj^ values for parton 
collisions, with transverse momentum p±_ is defined in the rest-frame of the 2-body interaction. If 
VKIN(4) < 0, the upper limit is inactive. For processes which are singular in the limit p_l — > (see 
VKIN(6)), VKIN(5) provides an additional constraint. The VKIN(3) and VKIN(4) limits can also be 
used in 2 ^ 1 ^ 2 processes. Here, however, the product masses are not known and hence assumed 
vanishing in the event selection. The actual p± range for massive products is thus shifted downwards 
with respect to the nominal one. 
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VKIN(5) : (D=0.5 GeV) lower cutoff on p± values, in addition to the VKIN(3) cut above, for processes 
which are singular in the limit p± ^ (see VKIN(6)). 

VKIN(6) : (D=l. GeV) 2-body parton collision processes, which do not proceed only via an intermediate 
resonance (i.e. arc 2 ^ 1 ^ 2 processes), arc classified as singular in the limit p± — > 0, if either or 
both of the two final state products has a mass m < VKIN(6). 

VKIN(7), VKIN(8) : (D=-100.,100.) range of allowed particle rapidities y* in the CM frame of the 
overall system. 

VKIN(9), VKIN(IO) : (D=0.,1000.) range (in GeV) of allowed particle transverse momenta k]_ per- 
pendicular to the z-axis in the global CM frame of the event. 

VKIN(ll), VKIN(12) : (D=-1000.,1000) range (in GeV-^) of allowed longitudinal positions z* of 
particles from the CM of the global collision system. 

VKIN(13), VKIN(14) : (D=0.,1000.) range (in GeV~^) of allowed transverse distances perpendic- 
ular to the z-axis in the CM frame of the global collision system. 

VKIN(21), VKIN(22) : (D=0.,1.) range of allowed xa (Bjorken-x) values for the parton of beam-side 
A that enters a parton collision. 

VKIN(23), VKIN(24) : (D=0.,1.) range of allowed xb (Bjorken-x) values for the parton of target-side 
B that enters a parton collision. 

VKIN(25), VKIN(26) : (D=-l.,l.) range of allowed Feynman xp values, where xp = xa — xb- 

VKIN(27), VKIN(28) : (D=-l.,l.) range of allowed cos((9) values in a 2 ^ 2 parton collision, where 6 
is the scattering angle in the rest-frame of the 2-body collision. 

VKIN(31), VKIN(32) : (D=2.,-l.) range of allowed m' values, where ml is the mass of the complete 
three- or four-body final state in 2 ^ 3 or 2 — > 4 processes (while rh (without a prime), constrained 
in VKIN(l) and VKIN(2), here corresponds to the one- or two-body central system). If VKIN(32) < 
0, the upper limit is inactive. 

VKIN(35), VKIN(36) : (D=0.,-1.) range of allowed abs(t) = -i values in 2 ^ 2 processes. Note that 
for deep inelastic scattering this is nothing but the scale of the photon vertex, in the limit that 
initial and final state radiation is neglected. If VKIN(36) < 0, the upper limit is inactive. 

VKIN(37), VKIN(38) : (D=0.,-1.) range of allowed abs(«) = -u values in 2 ^ 2 processes. If 
VKIN(38) < 0, the upper limit is inactive. 

VKIN(39) - VKIN(IOO) : currently not used. 

VKIN(lOl) - VKIN(200) : reserved for internal use of storing kinematical variables of an event. 



MSUB(ISUB) : array to be filled when MSTV(14) = 1 (for MSTV(14) = 0, the array MSUB is initialized 
in VNIXRIN automatically) to choose which subset of subprocesses to include in the generation. The 
ordering follows the ISUB code given before in Table 3. 
= : the subprocess ISUB is excluded. 
= 1 : the subprocess ISUB is included. 
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ISET(ISUB) : gives the type of kinematical variable selection scheme used for subprocess ISUB. 
= 1 : 2 — > 1 processes (parton fusion processes). 
= 2 : 2 — > 2 processes (hard and soft parton scattering processes). 
= 3 : 1 — > 2 processes (parton branching processes). 
= -1 : process is not yet implemented, but space is reserved. 
= -2 : process is not defined. 



3.9 Instructions on how to use the program 

There are a variety of example programs available from the website |http: / / rhic.phys.columbia.edu / rhic /vni . 
Below described is a very simple example, vnixple.f, disguised as a generic skeleton program that calls the 
steering routines VNIXRIN, VNIXRUN and VNIXFIN (c.f. Sec. 3.4) on a black-box-level. This example 
generates 1000 pp events at y/s = 546 GeV. Each event is evolved for a time range of - 30 fm in the pp 
center-of-mass frame (the global CM frame). 

PROGRAM VNIXPLE 

C... Purpose: example for p+p~ collisions at CERN collider with 546 GeV. 
INCLUDE 'vnil.inc' 
CHARACTER* ( 8 ) FRAME , BEAM , TARGET 



C . . . Required user input : 

C number of events and final time (in fm) . 

NEVT = 1000 
TFIN = 30. 

C colliding particles, global Lorentz frame, and total energy. 

BEAM = 'p+ 
TARGET= 'p- 
FRAME = 'cms ' 
WIN = 546. 

C Specify your change of default parameters here: 



C 



C. . . 



Initialize simulation procedure. 

CALL VNIXRIN (NEVT , TFIN , FRAME , BEAM , TARGET , WIN) 



C. . .Begin loop over of events. 
DO 500 IEVT=1,NEVT 



C 



Generate collision event. 



CALL VNIXRUN (lEVT, TFIN) 



C 
C 



Analysis of event as-a whole. 



C. . .End loop over events. 



500 CONTINUE 
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C. . .Finish up simulation procedure. 
CALL VNIXFINO 



END 

This example is rather self-explanatory, and will become further clear with the following subsections. 
The " . . ." indicate where one can interface on a black-box- level with the program as-a- whole without needing 
to dig deeper into its subroutine structure. In particular, one can insert here histograming (either using the 
included portable package "vnibook.f, or any other package such as the HBOOK or PAW of the CERN 
library). Note also the pre-programmed histogram options discussed in Sec. 3.7 (switches MSTV(150)- 
MSTV(200)) and Appendix B (subroutines VNIANA1-VNIANA4) . 

In the above example only the final state at the end of each event is collected. If the user is interested 
in the actual space-time development of the collision events, he/she can access this by inserting in the loop 
over events the following extension: 

C. . .Generate collision event; analyze evolution every 5 fm. 
TSTEP=5 . 
NTIM=20 

C Begin loop over time intervals. 

DO 400 ITIM=1,NTIM 
TRUN=ITIM*TSTEP 

IF (TRUN . GT . TFIN . OR . ITIM . EQ . NTIM) TRUN=TFIN 
CALL VNIXRUN(IEVT,TRUN) 

C Analysis of time development. 

C 

C End loop over time intervals. 

400 IF (TRUN. EQ. TFIN) GOTO 450 
450 CONTINUE 

Remarks: 

(i) In case one desires to interface with the PYTHIA program, one should write a little program with 

common block 

C0MM0N/LUJETS/N,K(4000,5) ,P(4000,5) ,V(4000,5) 

and which copies the variables N,K,P, and V of the particle record VNIREC1/VNIREC2 directly. 
Be aware however, that VNIREC1/VNIREC2 is dimensioned for 100000 entries, whereas LUJETS has 
space for only 4000 entries. 

(ii) In case one desires to interface with the HERWIG program, one should write a little program with 

common block 

PARAMETER (NMXHEP=2000) 
COMMON/HEPEVT/NEVHEP , NHEP , ISTHEP (NMXHEP) , IDHEP (NMXHEP) , 
& JMOHEP (2 , NMXHEP) , JDAHEP (2 , NMXHEP) , PHEP (5 , NMXHEP) , VHEP (4 , NMXHEP) 

and which calls the routine RWREF C{MCONV, IP, NP) with MCONV = 1,IP= l,and NP = N. 
Again, watch the dimensions, 100000 versus 2000. 
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3.10 Example for a typical collision event 



As a result of running the above example program vnixple.f, a number of files are created in the directory, 
all of which are called VNI???.DAT, where "???" arc cither 3 letters, or 3 digits. The latter arc data files 
that contain event statistics and other information on the performance of a run. (see Sec. 3.7, switches 
MSTV(150)-MSTV(200) and Appendix B for further options). By default only a few files are created: 

VNIRUN.DAT containing the main summary of a simulation with print-out of results, 

VNICOL.DAT containing a detailed summary of all parton collisions that occurred, 

VNITIM.DAT containing the generated time-step grid and its real-time translation, 

VNIRLU.DAT containing the status of the random number generator, 

VNIERR.DAT containing error messages in case problems occur. 

Most important is the file VNIRUN.DAT, which gives a summary of what occurred during a simulation 

on the parton cascade level, the cluster formation level and the hadron decay level. One can obtain a listing 
of an event by simply calling VNILIST(l) at the end or during an event, which gives a first impression (for 
a complete listing see Appendix D): 

(i) The event record begins with the beam/target particles, here p and p, which are then transformed to 

the center-of-mass frame (here it coincides with the global frame of reference). 

(ii) Then follow the three valence quarks of p, and after that its intrinsic gluons and seaquark-pairs. After 

that the initial-sate valence quarks and gluons and seaquarks of p are listed. 

(iii) Then comes the history of the cascade evolution, including hard scatterings with shower 
(bremsstrahlung) , subsequent cluster-formation of the materialized partons and hadron decay of these 
parton clusters. 

(iv) Finally, all non-materialized (non-interacted) partons are recombined to a beam- and a target cluster 
that undergo a soft fragmentation, again via cluster-hadron decay. 

(v) At the very end the sum of charge, three-momentum, energy, and the total invariant mass of the collision 

system is printed. 



Example of the particle record for a pp collision event at -s/s = 546 GeV. 
Particle listing (summary: momentum space P in GeV) 







event 


no . 


1 


time step 


200 


at 


time 


34. 


,75 fm 




I 


particle 


KS 


KF 


C A 


P_x 


P. 


-y 


P. 


.z 


E 


M 


1 


(p+) 


15 


2212 





.000 




.000 


272, 


,998 


273.000 


.938 


2 


(p~-) 


15 


-2212 





.000 




.000 


-272. 


,998 


273.000 


.938 


3 


(p+) 


15 


2212 





.000 




.000 


272. 


,998 


273.000 


.938 


4 


(p~-) 


15 


-2212 





.000 




.000 


-272. 


,998 


273.000 


.938 


5 


(u) 


13 


2 


1 


.650 




.048 


1. 


,786 


2.319 


.006 


6 


(d) 


13 


1 


1 


-.593 




.165 


2. 


,092 


2.587 


.010 


7 


(d~) 


13 


-1 


3 


-.135 




.452 


6. 


,724 


6.886 


-.490 


8 


(d) 


13 


1 


3 


.413 




.211 




,408 


.935 


-.484 



64 



9 (g) 



14 21 1 1 1.751 .556 2.232 3.398 -.682 



48 (d~) 

49 (u~) 

50 (u) 

51 (g) 



13 
13 
13 
13 



-10 2 

-2 3 

2 3 

21 3 2 



.526 
.121 
.252 
.110 



.150 
.114 
.044 
.229 



-16.659 
-4.941 
-1.576 

-12.733 



15.157 
4.106 
.935 
11.449 



.010 
-.174 
-.249 
-.234 



83 (g) 


14 


21 


2 


1 


2 


.052 


2, 


.858 


1. 


,999 


1. 


,045 


-3, 


,909 


84 (d) 


14 


1 


1 







.000 




.000 


-8. 


.008 


8, 


,008 




,000 


85 (d) 


14 


1 


2 







.000 




.000 


8, 


,008 


8, 


,008 




,000 


86 (CMshower) 


11 


94 








2 


.661 


3 


.728 


-3. 


,710 


8, 


,110 


5, 


,571 


87 (d) 


14 


1 


2 





1, 


.128 




.073 


-1, 


,364 


3, 


,504 


3, 


.023 


88 (cluster) 


11 


91 








1, 


.116 


1, 


.615 


-2, 


,687 


3, 


,412 




,750 


89 (cluster) 


11 


91 










.450 




.605 




,228 


1, 


,213 




,922 


90 (cluster) 


11 


91 








-1 


.528 




.662 


1. 


,381 


2. 


,181 




,280 


91 (cluster) 


11 


91 








1, 


.252 




.689 


-1. 


,272 


3. 


,136 


2, 


.484 


92 (cluster) 


11 


91 










.110 




.650 


5, 


,219 


5, 


,796 


2, 


.433 


93 piO 


6 


111 










.660 




.523 




,880 


1. 


,171 




,135 


94 piO 


6 


111 










.461 


1, 


.090 


-2, 


,094 


2. 


,377 




,135 


95 rho- 


6 


-213 










.334 




.562 




,332 


1, 


,068 




,769 


96 piO 


6 


111 










.110 




.046 




,182 




,186 




,135 


97 pi- 


6 


-211 










.678 




.312 




,548 


1, 


,066 




,140 


98 pH- 


6 


211 










.845 




.352 




,547 


1, 


,203 




,140 



109 (Beam-REM) 


11 


92 










.315 




.737 


274. 


,428 


281, 


.871 


70, 


, 164 


110 (Targ-REM) 


11 


93 










.907 


-1 


.587 


-276. 


,508 


247, 


.044 


-116, 


.181 


111 (CMF) 


15 


100 








-1 


.222 




.850 


-2. 


,080 


528, 


.915 


528, 


.909 


112 (cluster) 


11 


91 










.315 




.737 


274. 


,428 


281, 


.871 


70, 


.164 


113 (cluster) 


11 


91 










.907 


-1 


.587 


-276. 


,508 


247, 


.044 


-116, 


.181 


114 (Beam-REM) 


11 


92 










.411 




.648 


127. 


,485 


127, 


.543 


3, 


,762 


115 (a_l-) 


17 


-20213 










.493 




.228 


19. 


,200 


19, 


.250 


1, 


.275 


116 (Delta^■) 


17 


2214 










.081 




.420 


108. 


,286 


108, 


.293 


1, 


.232 


117 piO 


7 


111 










.490 




.124 


4. 


,473 


4, 


.831 




,135 


118 (rho-) 


17 


-213 










.005 




.102 


14. 


,584 


14, 


.605 




,769 


119 piO 


7 


111 










.160 




.015 


35. 


,475 


37, 


.042 




,135 


120 p+ 


7 


2212 










.085 




.402 


72. 


,525 


75, 


,579 




,938 


121 piO 


7 


111 










.122 




.268 


8, 


,519 


9, 


.014 




,135 


122 pi- 


7 


-211 










.133 




.367 


5. 


,779 


6, 


.174 




,140 


123 (Targ-REM) 


11 


93 








-1 


.252 




.510 


-63, 


,766 


63, 


.825 


2, 


.384 


124 (Delta~-) 


17 


-2214 










.779 




.303 


-21. 


,201 


21, 


,253 


1, 


,232 


125 (eta) 


17 


221 










.473 




.208 


-42. 


,565 


42, 


,572 




,549 


126 piO 


7 


111 










.072 




.254 


-3. 


,205 


3, 


, 199 




,135 


127 p~- 


7 


-2212 










.701 




.052 


-18. 


,281 


18, 


,903 




,938 


128 piO 


7 


111 










.179 




.069 


-12. 


,277 


12, 


.622 




,135 


129 piO 


7 


111 










.060 




.060 


-12. 


,901 


13. 


.269 




,135 


130 piO 


7 


111 










.225 




.220 


-17. 


,816 


18. 


.383 




,135 



65 



336 (rho+) 


17 


213 








-.263 


.404 


-5. 


.474 


5, 


.549 


.769 


337 (omega) 


17 


223 








.067 


-.237 


-31 


.485 


31, 


.496 


.783 


338 piO 


7 


111 








-.342 


.522 


-5. 


.018 


5, 


.113 


.135 


339 pi+ 


7 


211 








.085 


-.120 




.742 




.657 


.140 


340 piO 


7 


111 








-.309 


.029 


-10 


.513 


10, 


.790 


.135 




svun: 


.00 






.000 


.000 




.000 


546, 


.000 


546 . 000 



More information on pre-programmed print-out of results concerning for particle distributions, time 
evolution, observables, etc., can be switched on with switches MSTW(151) - MSTW(200), as explained in 
the previous Sec. 3.7. 
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A Further physics routines 



Aside from the three main routines VNIXRIN, VNIXRUN, and VNIXFIN, the physics routines described 
below form the heart of the program. In general, for each of the included collision processes of Table 1, there 
is an initialization routine VNI??IN that sets up the initial state, and an evolution routine VNIEV?? that 
carries out the time evolution of the particle distributions in phase-space, starting from the initial state. The 
routines VNICOLL, VNISBRA, VNITBRA, etc., therein perform the perturbative parton evolution in terms 
of parton collisions, space-like branchings, and timc-likc branchings, respectively. The routines VNICHAD, 
VNIBHAD, etc., deal with the cluster hadronization. These routines are universal and independent of the 
process under consideration. 

SUBROUTINE VNIC0L2(NC0L) 

Purpose: to test for all possible pairings of 'alive' partons, whether a collision is allowed, and, if it is allowed, 
to evaluate the quantities characterizing the 2-body collision at the parton level according to the cross-section 
and to choose one of the possible subprocesses (channels) . 

NCOL : total number of found collisions. 



SUBROUTINE VNICOLL (ISUB , Q2 , IPl , IP2 , IPS , IP4) 

Purpose: to perform a 2-body parton collision, find outgoing flavors and colors and calculate kinematics 
and parton-parton collision. 

ISUB : Subprocess of 2-body collision (c.f. Table 3). 

Q2 : Momentum transfer scale of subprocess (c.f. MSTV(43), Sec. 3.7). 

IPl, IP2 : Line numbers of paerons entering the interaction vertex. 

IPS, IP4 : Line numbers of particles emerging from the interaction vertex. For 2—^1 processes I PA is 
equal to zero. 



SUBROUTINE VNISBRA ( IP 1,IP2,QMAX) 

Purpose: to generate space-like parton evolution, by backward tracing of branching processes down to some 
initial resolution scale associated with a non-pcrturbative cut-off. The performance is regulated by the 
switches MSTV(50)-MSTV(52) and parameters PARV(51)-PARV(53) (c.f. Sec. 3.7). 

IPl, IP2 : partons in hues IPl and IP2 in the commonblocks VNIREC1/VNIREC2, are evolved by 
space-like parton branchings from initial mass scale QMAX. If one of the particles is a non-radiating 
one, then only one branch is evolved, and the spectator particle role is then only to take up recoil due 
to emitted particles by the radiating one. 

QMAX : initial mass scale, from which the two particles are evolved from Q = QMAX backwards to 
lower Q < QMAX until the cut-off Q = go =PARV(52) is reached. 



SUBROUTINE VNITBRA ( IP 1,IP2, QMAX) 

Purpose: to generate time-like parton evolution, by coherent branching processes within a finite time interval, 
Q ^ QQi Q ~^ 1^ 9 ^ 99i 9 ^ QQ- The performance is regulated by the switches MSTV(61)-MSTV(65) and 
parameters PARV(61)-PARV(65) (c.f. Sec. 3.7). 
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IPl > IP2 — : generate a time-like parton branching for the parton in line IPl in commonblocks 
VNIREC1/VNIREC2, with maximum allowed mass QMAX. 

IPl > IP2 > : generate time-like parton branchings for the two partons in lines IPl and IP2 in the 
commonblocks VNIREC1/VNIREC2, with maximum allowed mass for each parton QMAX. 

QMAX : the maximum allowed mass of a radiating parton, i.e. the starting value for the subsequent 
evolution. (In addition, the mass of a single parton may not exceed its energy, the mass of a parton 
in a system may not exceed the invariant mass of the system.) 



SUBROUTINE VNICHAD(I1 , 12) 

Purpose: to search at a given point of time the particle record for possible clustering of materialized partons, 
which are those which have been struck out from the initial state or have been produced as secondaries by 
parton collisions or radiative processes. It consists of checking for each possible pair of two partons, the 
invariant space-time distance between them, their flavor and origin, the matching of color and anticolor, 
their common invariant mass, etc.. If a pair satisfies all necessary conditions for the formation of one or two 
color singlet clusters, then they are combined and converted into a composite cluster with internal flavor 
content and invariant mass determined by the partons, which then decays either into a pair of color singlet 
clusters, or into a singlet cluster plus an additional gluon or quark. 

II, 12 : first and last entry of particle record VNIREC1/VNIREC2 in between the cluster formation is 
carried out by looping over all particles I with II < I < 12. 



SUBROUTINE VNICFRA(I1 , 12) 

Purpose: to fragment color singlet clusters with given qq substructure into final state hadrons. For each 
cluster its invariant mass and flavor content are checked, and it is determined whether it decays isotropically 
into a pair of hadrons, or converts directly into a single hadron of corresponding flavor. The hadron decay 
is described by a phase-space model, in which the probability for forming a certain hadron state is given by 
the density of hadronic states with mass below the cluster mass. Once a hadron state is chosen, it cither 
propagates on undisturbed if it is a stable particle, or it decays further according to the particle data tables, 
if it is a resonance. 

Jl > : entry of particle record VNIREC1/VNIREC2 from which analysis starts. 
J2 > 72 : entry of particle record up to which analysis extends. 



SUBROUTINE VNIBHAD(I1 , 12) 

Purpose: to perform beam cluster formation and hadronization. The remnant partons of beam/target 
particles, those that originate from the partonic coherent initial state in collisions involving hadrons or 
nuclei and have not interacted with decohereing effect. These primary partons are reassembled to clusters 
along the beam axis, which are then fragmented based on the minimum bias model of the UA5 collaboration, 
such that color, flavor and energy-momentum of the collision system as-a-whole are conserved. 

II, 12 : first and last entry of particle record VNIREC1/VNIREC2 in between the soft beam/target frag- 
mentation is carried out by looping over all particles / with II < I < 12. 
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SUBROUTINE VNIBFRA(I1 , 12) 



Purpose: to hadronize beam/target clusters consisting of remnant beam partons, using parts of the HERWIG 
5.7 program. Beam remnants are those that emerge from the coherent partonic initial state without having 
undergone any decohering interaction, and carry the remaining energy, flavor and color of the system. The 
beam remnants are lumped to beam clusters which are fragmented similarly as the clusters arising from the 
parton cascade evolution. 

/I > : entry of particle record VNIREC1/VNIREC2 from which analysis starts. 
/2 > /I : entry of particle record up to which analysis extends. 



SUBROUTINE VNIPADC(I1,I2,NUST) 

Purpose: to decay unstable particles resulting from cluster decays to hadrons. Out of these primary hadrons, 
those that are unstable resonances are decayed according to the Particle Data, which gives lower mass 
hadrons plus possibly gauge bosons and leptons. These secondary decay products may decay further. 

Jl > : entry of particle record VNIREC1/VNIREC2 from which analysis starts. 
12 > II : entry of particle record up to which analysis extends. 
NUST : returns number of decay products that are still unstable. 
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B Event analysis routines 



The following routines collect, analyze and print out information on the result of a simulation, and give the 
user immediate access to calculated particle data, spectra, and other observables. 

• The routine VNILIST provides a useful tool to obtain a quick overall view of the state of the particle 
record at any time during the simulation of one or many events, and (in parts) the space-time history 
of the system. 

• The routine VNIEDIT allows to edit the particle record, e.g. to compress, or to exclude certain particle 
species from the record. It is used frequently during the simulation procedure to clean up used but 
further on unecessary particle information. 

• Via the functions KVNI and PVNI, values of some frequently appearing variables may be obtained 
more easily. 

• The routines VNIANAL and VNIANAl - VNIANA4 carry out various data accumulation duties, 
and print out tables at the end of the simulation, which are automatically directed to files named 
VNI???.DAT where ??? is a 3-digit number, specified below. 

• Finally, the six routines VNISPHE, VNITHRU, VNICLUS, VNICELL, VNIJMAS and VNIFOWO 
allow the possibility to find some global event shape properties as discussed in detail in the PYTHIA 
documentation |Q. 

• The common block VNIANA allows to access the bin sizes, and bin ranges of histograms in the routines 
VNIANAl - VNIANA4. 



C0MM0N/VNIANA/MANA(50) ,PANA(50) 



Purpose: to provide data for variable ranges, bin sizes, etc., for histograms and tables created by 
analysis routines VNIANAl - VNIANA4. 



MANA(l) - MANA(IO) : (D=l, 21,41,61,91, 121, 151, 181, 221,301) Sequential numbers of time 
steps, at which a snapshot analysis during the space-time evolution of the collision system may 
be performed in routines VNIANA4 - VNIANA4. 

MANA(ll) : (D=40) Number of bins for energy distribution in ln(l/x) (where x is the energy 
fraction) in routine VNIANA2. 

MANA(12) : not used. 

MANA(13) : (D=40) Number of bins for pair mass distribution of Bose-Einstein correlations in 
routine VNIANA2. 

MANA(14)-MANA(20) : not used. 

MANA(21) : (D=20) Number of bins for cluster y-distribution in routine VNIANA3. 

MANA(22) : not used. 

MANA(23) : (D=20) Number of bins for cluster r^-distribution in routine VNIANA3. 

MANA(24) : not used. 

MANA(25) : (D=20) Number of bins for cluster p_L-distribution in routine VNIANA3. 

MANA(26) : not used. 

MANA(27) : (D=20) Number of bins for cluster r_L-distribution in routine VNIANA3. 

MANA(28) -MANA(30) : not used. 

MANA(31) : (D=20) Number of bins for parton/hadron y-distribution in routine VNIANA4. 

MANA(32) : not used. 
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MANA(33) : (D=20) Number of bins for parton/hadron r^-distribution in routine VNIANA4. 
MANA(34) : not used. 

MANA(35) : (D=20) Number of bins for parton/hadron p^-distribution in routine VNIANA4. 
MANA(36) : not used. 

MANA(37) : (D=20) Number of bins for parton/hadron r^-distribution in routine VNIANA4. 
MANA(38) - MANA(50) : not used. 

PANA(l) - PANA(IO) : not used. 

PANA(l) ,PANA(12) : (D=0.,6.) Lower and upper bound for energy distribution in ln(l/a;) 
(where x is the energy fraction of particles), in routine VNIANA2. 

PANA(13) ,PANA(14) : (D=0.,3. GeV) Lower and upper bound for pair mass distribution of 
Bose-Einstein correlations in routine VNIANA2. 

PANA(15) - PANA(20) : not used. 

PANA(21) ,PANA(22) : (D=-8.,8.) Lower and upper bound for cluster y-distribution in routine 
VNIANA3. 

PANA(23) ,PANA(24) : (D=-40.,40. GeV"^) Lower and upper bound for cluster r;j-distribution 
in routine VNIANA3. 

PANA(25) ,PANA(26) : (D=0.,10. GeV) Lower and upper bound for cluster p_L-distribution in 
routine VNIANA3. 

PANA(27) ,PANA(28) : (D=0.,40. GeV"^) Lower and upper bound for cluster r_L-distribution 
in routine VNIANA3. 

PANA(29) - PANA(30) : not used. 

PANA(31) ,PANA(32) : (D=-8.,8.) Lower and upper bound for parton/hadron y-distribution in 

routine VNIANA4. 

PANA(33) ,PANA(34) : (D=-100.,100. GeV"^) Lower and upper bound for parton/hadron re- 
distribution in routine VNIANA4. 

PANA(35) ,PANA(36) : (D=0.,10. GeV) Lower and upper bound for parton/hadron p^- 

distribution in routine VNIANA4. 

PANA(37) ,PANA(38) : (D=0.,100. GeV"^) Lower and upper bound for parton/hadron re- 
distribution in routine VNIANA4. 

PANA(39) - PANA(50) : not used. 
SUBROUTINE VNILIST(MSEL) 

Purpose: to list the current state of the particle record during the evolution of an event, or certain specified 
excerpts of it. Also, to list parton or particle data, or current parameter values. 

MSEL : determines what is to be listed. 

= : writes a logo with program version number and last date of change; is mostly for internal use. 
= 1 : gives a simple list of current event record, in an 80 column format suitable for viewing directly 
on the computer terminal. For each entry, the following information is given: the entry number I, 
the parton/particle name, the status code KS (K(I,1)), the flavour code KF (K(I,2)), the color and 
anticolor, C (L(I,1)), and A (L(I,2)), and the three-momentum, energy and mass (P(I,1)-P(I,5)). A 
final line contains information on total charge, momentum, energy and invariant mass. 
= 2 : gives a more extensive list of the current event record, in a 132 column format, suitable for 
printers or workstations. For each entry, the following information is given: the entry number I, the 
parton/particle name (with padding as described for MSEL = 1), the status code KS (K(I,1)), the 
flavour code KF (K(I,2)), the color and anticolor, C (L(I,1)) and A (L(I,2)), the origin of the particle 
in a 3 digit N1N2N3, where A^l is the event and N2 the time step in this event, when the particle 



71 



was first produced, and N3 is the number of the entry in the particle record that the particle occupied 
when it first appeared. This labeling allows a unique specification of its 'historical' origin. Finally, it 
follows the three-momentum, energy and mass (P(I,1) - P(I,5)). A final line contains information on 
total charge, momentum, energy and invariant mass. 

= 3 : gives the same basic listing as = 2, but with an additional line for each entry containing infor- 
mation on production vertex position and time (V(I,1) - V(I,4)). 

= -1 : gives the format of listing as = 1, but instead of momentum, energy and mass, it lists now the 
particle's position in the overall reference frame, (R(I,1) - R(I,3)), its momentum-rapidity (Y(P)) and 
its transverse momentum (FT) with respect to the z-axis, which is usually the jet-axis, or beam axis. 
= -2 : gives the format of listing as = 2, but as for = -1, instead of momentum, energy and mass, 
now the particle's positions, rapidities and transverse momenta in the overall Lorentz frame are listed. 
As for = 2, additional information on the 'historical' origin of the particle is listed. 
= -3 : gives the format of listing as = 3, again, instead of momentum, energy and mass, now the 
particle's positions, rapidities and transverse momenta are listed. 

= 11 : provides a simple list of all parton/particlc codes defined in the program, with KF code 
and corresponding particle name. The list is grouped by particle kind, and only within each group in 
ascending order. 

= 12 : provides a list of all parton/particle and decay data used in the program. Each parton/particle 

code is represented by one line containing KF flavour code, KC compressed code, particle name, an- 

tiparticle name (where appropriate), electrical and color charge, mass, resonance width and maximum 

broadening, average invariant lifetime and whether the particle is considered stable or not. 

= 13 : gives a list of current parameter values for MSTV, FARV, and for MSTW and FARW. This 

is useful to keep check of which default values were changed in a given run. 

= 14 : gives a list of time grid and size of time slizes used in the current run. 



SUBROUTINE VNIEDIT(MEDIT) 

Purpose: to exclude unstable or undetectable partons/particles from the event record. One may also use 
VNIEDIT to store spare copies of events (spcciflcally initial parton conflguration) that can be recalled to 
allow e.g. different fragmentation schemes to be run through with one and the same parton configuration. 
Finally, an event which has been analysed with VNISFHE, VNITHRU or VNICLUS (see below) may be 
rotated to align the event axis with the z-direction. 

MEDIT : tells which action is to be taken. 

= : empty {KS = 0) and documentation {KS > 20) lines arc removed. The partons/particles 
remaining are compressed in the beginning of the VNIRECl and VNIREC2 commonblocks and the 
N value is updated accordingly. The event history is lost, so that information stored in K(I,3), K(I,4) 
and K(I,5) is no longer relevant. 

= 1 : as = 0, but in addition all partons/particles that have fragmented/decayed {KS > 10) are 
removed. 

= 2 : as = 1, but also all neutrinos and unknown particles (i.e. compressed code KC = 0) are 
removed. 

= 3 : as = 2, but also all uncharged, color neutral particles are removed, leaving only charged, stable 
particles (and unfragmented partons, if fragmentation has not been performed). 

= 5 : as = 0, but also all partons which have branched or been rearranged in a parton shower and 
all particles which have decayed are removed, leaving only the fragmenting parton configuration and 
the final state particles. 

= 11 : remove lines with K(I,1) < 0. Update event history information (in K(I,3) - K(I,5)) to refer 
to remaining entries. 

= 12 : remove lines with K(I,1) = 0. Update event history information (in K(I,3) - K(I,5)) to refer 
to remaining entries. 

= 13 : remove lines with K(I,1) = 11, 12 or 15, except for any hue with K(I,2) = 94. Update event 
history information (in K(I,3) - K(I,5)) to refer to remaining entries. In particular, try to trace origin 
of daughters, for which the mother is decayed, back to entries not deleted. 

= 14 : remove lines with K(I,1) = 13 or 14, and also any line with K(I,2) = 94. Update event history 
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information (in K(I,3) - K(I,5)) to refer to remaining entries. In particular, try to trace origin of 

rearranged partons back through the parton shower history to the shower initiator. 

= 15 : remove lines with K(I,1) > 20. Update event history information (in K(I,3) - K(I,5)) to refer 

to remaining entries. 

= 21 : all partons/particles in current event record are stored (as a spare copy) in bottom of com- 
monblocks VNIREC1/VNIREC2. 

= 22 : partons/particles stored in bottom of event record with = 21 are placed in beginning of record 

again, overwriting previous information there (so that e.g. a different fragmentation scheme can be 
used on the same partons). Since the copy at bottom is unaffected, repeated calls with = 22 can be 
made. 

= 23 : primary partons/particles in the beginning of event record arc marked as not fragmented or 

decayed, and number of entries N is updated accordingly. Is simple substitute for = 21 plus = 22 

when no fragmentation/decay products precede any of the original partons/particles. 

= 31 : rotate largest axis, determined by VNISPHE, VNITHRU or VNICLUS, to sit along the z- 

direction, and the second largest axis into the xz plane. For VNICLUS it can be further specified to 

+2 axis and xz plane with x > 0, respectively. Requires that one of these routines has been called 

before. 

= 32 : mainly intended for VNISPHE and VNITHRU, this gives a further alignment of the event, in 
addition to the one implied by = 31. The "slim" jet, defined as the side (z > or z < 0) with the 
smallest summed over square root of number of particles, is rotated into the +z hemisphere. In 
the opposite hemisphere (now z < 0), the side of x > and a; < which has the largest summed 
absolute is rotated into the z < 0, x > quadrant. Requires that VNISPHE or VNITHRU has been 
called before. 



FUNCTION KVNId.J) 

Purpose: to provide various integer-valued event data. Note that many of the options available (in particular 
/ > 0, J > 14) which refer to event history will not work after a VNIEDIT call. 

I = 0, J = : properties referring to the complete event. 
= 1 : N, total number of lines in event record. 

= 2 : total number of partons/particles remaining after fragmentation and decay. 
— 6 : three times the total charge of remaining (stable) partons and particles. 

/ > 0, J = : properties referring to the entry in line no. / of the event record. 

= 1 - 5 : K(I,1) - K(I,5), i.e. parton/particle status KS, flavour code KF and origin/decay 

product/color-flow information. 

= 6 : three times parton/particle charge. 

= 7 : 1 for a remaining entry, for a decayed/fragmented/ documentation entry. 

= 8 : KF code (K(I,2)) for a remaining entry, for a decayed/fragmented/documentation entry. 

= 9 : KF code (K(I,2)) for a parton (i.e. not color-neutral entry), for a particle. 

= 10 : KF code (K(I,2)) for a particle (i.e. color- neutral entry), for a parton. 

= 11 : compressed flavour code KC. 

= 12 : color information code, i.e. for color neutral, 1 for color triplet, -1 for antitriplet and 2 for 
octet. 

= 13 : flavour of "heaviest" quark or antiquark (i.e. with largest code) in hadron or diquark (including 
sign for antiquark), else. 

= 14 : generation number. Beam particles or virtual exchange particles are generation 0, original 
partons/particles generation 1 and then 1 is added for each step in the fragmentation/decay chain. 
= 15 : line number of ancestor, i.e. predecessor in first generation (generation entries are disre- 
garded) . 

= 16 : rank of a hadron in a jet system it belongs to. Rank denotes the ordering in flavour space, with 
hadrons containing the original flavoTir of the jet having rank 1, increasing by 1 for each step away 
in flavour ordering. All decay products inherit the rank of their parent. Whereas the meaning of a 
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first-rank hadron in a quark jet is always well-defined, the definition of higher ranks is only meaningful 
for independently fragmenting quark jets. In other cases, rank refers to the ordering in the actual 
simulation, which may be of little interest. 

= 17 : generation number after a collapse of a jet system into one particle, with for an entry not 
coming from a collapse, and -1 for entry with unknown history. A particle formed in a collapse is 
generation 1, and then one is added in each decay step. 

= 18 : number of decay /fragmentation products (only defined in a collective sense for fragmentation). 

= 19 : origin of color for showering parton, else. 

= 20 : origin of anticolor for showering parton, else. 

= 21 : position of color daughter for showering parton, else. 

= 22 : position of anticolor daughter for showering parton, else. 

= 23 : number of 2-body collisions of parton/particle. 

= 24 : number of space-like decays of parton/particle. 

= 25 : number of time- like decays of parton/particle. 

= 27 : entry of direct mother of parton/particle. 

= 28 : timestep of production of parton/particle. 

— 29 : index of generation in the cascade tree of parton/particle. 



FUNCTION PVNI(I,J) 

Purpose: to provide varioTis real- valued event data. Note that some of the options available (/ > 0, J = 
20 — 25) are primarily intended for studies of systems in their respective CM frame. 

/ = 0, J = : properties referring to the complete event. 

= 1 - 4 : sum of Px,Py,Pz and E, respectively, for the stable remaining entries. 

— 5 : invariant mass of the stable remaining entries. 

= 6 : sum of electric charge of the stable remaining entries. 

J > 0, J = : properties referring to the entry in line no. / of the event record. 

= 1 - 5 : P(I,1) - P(I,5), i.e. noimally Px,Py,Pz, E and m for parton/particle. 
= 6 : electric charge q. 

— 7 : momentum squared p^ = Px + Py + p1- 



= 10 : transverse momentum p^. 

= 11 : transverse mass squared m']_ = vr? +p'^+Py. 

— 12 : transverse mass m±. 

13 - 14 : polar angle 9 in radians (between and tt) or degrees, respectively. 
= 15 - 16 : azimuthal angle (f) in radians (between — tt and tt) or degrees, respectively. 
= 17 : true rapidity y = ^ln{{E + pz)/{E —Pz))- 

= 18 : rapidity y.^ obtained by assuming that the particle is a pion when calculating the energy E, 
to be used in the formula above, from the (assumed known) momentum p. 
= 19 : pseudorapidity rj = ^ ln((p -|- pz) /{p — Pz))- 

— 20 : momentum fraction Xp = 2p/W, where W is the total energy of initial parton/particle 
configuration. 

= 21 : Xp = 2pz/W (Feynman-x if system is studied in CM frame). 



= 22 : x_L = 2p±/W. 
^23 : XE = 2E/W. 
^24: z+ = {E + pz)/W. 
= 25 : z_ = {E-p^)/W. 



SUBROUTINE VNIANAL(MSEL) 
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Purpose: to administrate in each event, which of the available analysis duties (described below) are carried 
out, and when during the time evolution of the particle system data samples are taken and accumula- ted 
to perform a statistical analysis at the end. It calls the routines VNIANAl - VNIANA4 to do this job, for 
the options that are switched on by default or by the user. The switches and their effect are described in 
Sec. 3.7. After the last event of a run, the results are automatically printed out in tabular form to files 
called VNI???.DAT, with numbers ???, corresponding to the selected analysis options. 

MSEL : gives the type of action to be taken. 

= : all arrays for data analysis in VNIANAl - VNIANA4 are reset. 

= 1 : data samples are accumulated, stored and analyzed in VNIANAl - VNIANA4. 

= 10 : final results are written to output files called VNI.???.DAT, where '???' is a 3-digit number 

specified in the descriptions below. 



SUBROUTINE VNIANAl (MSEL) 

Purpose: to provide a number of global analyses referring to time- independent quantities of events as- 
a- whole, as e.g. total charged multiplicities, etc.. Is adopted from PYTHIA and only slightly altered. 
Accumulated statistics is written out at the end of a run. When errors are quoted, these refer to the 
uncertainty in the average value for the event sample as a whole, rather than to the spread of the individual 
events, i.e. errors decrease like one over the square root of the number of events analyzed. 

MSEL : determines which action is to be taken. Generally, a last digit equal to indicates that the 
statistics counters for this option is to be reset. Last digit 1 leads to an analysis of current event with 
respect to the desired properties. The statistics accumulated is output in tabular form with last digit 
2, and written to disc. 

= 10 : statistics on parton multiplicity is reset. 

= 11 : the parton content of the current event is analyzed, classified according to the flavour content 

of the hard interaction and the total number of partons. 

= 12 : gives a table on parton multiplicity distribution, written to file VNI110.DAT. 
= 20 : statistics on particle content is reset. 

= 21 : the particle/parton content of the current event is analyzed, also for particles which have 
subsequently decayed and partons which have fragmented. Particles are subdivided into primary and 
secondary ones, the main principle being that primary particles are those produced directly, while 
secondary come from decay of other particles. 

= 22 : gives a table of particle content in events, written to file VNI120.DAT. 
= 30 : statistics on factorial moments is reset. 

= 31 : analyzes the factorial moments of the multiplicity distribution in different bins of rapidity and 
azimuth. 

= 32 : prints a table of the first four factorial moments for various bins of pseudorapidity and azimuth, 
written to file VNI130.DAT. The moments are properly normalized so that they would be unity (up 

to statistical fluctuations) for uniform and uncorrelated particle production according to Poissonian 
statistics, but increasing for decreasing bin size in case of "intermittent" behaviour. 

= 40 : statistics on energy-energy correlation is reset. 

= 41 : the energy-energy correlation (EEC) of the current event is analyzed. Events are assumed 
given in their CM frame. The weight assigned to a pair i and j is 2EiEj/W'^, where W is the sum of 
energies of all analyzed particles in the event. Energies are determined from the momenta of particles 
and masses. Statistics is accumulated for the relative angle 9ij, ranging between and 180 degrees, 
subdivided into 50 bins. 

= 42 : prints a table of the energy-energy correlation (EEC) and its asymmetry (EEC A), with er- 
rors, written to file VNI140.DAT. The definition of errors is not unique. Each event is viewed as one 
observation, i.e. an EEC and EECA distribution is obtained by summing over all particle pairs of an 
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event, and then the average and spread of this cv cul-distrib ution is calculated in the standard fashion, 
with the quoted error containing a factor of 1/ y/no. events. 

= 50 : statistics on complete final states is reset. 

= 51 : analyzes the particle content of the final state of the current event record. During the course 
of the run, statistics is thus accumulated on how often different final states appear. Only final states 
with up to 8 particles are analyzed, and there is only reserved space for up to 200 different final states. 
Most high energy events have multiplicities far above 8, so the main use for this tool is to study the 
effective branching ratios obtained with a given decay model for e.g. charm or bottom hadrons. 
— 52 : gives a list of the (at most 200) channels with up to 8 particles in the final state, with their 
relative branching ratio, written to file VNI150.DAT. The ordering is according to multiplicity, and 
within each multiplicity according to an ascending order of KF codes. The KF codes of the particles 
belonging to a given channel are given in descending order. 



SUBROUTINE VNIANA2 (MSEL) 

Purpose: to analyze the time evolution of particles yields, energy ad momentum composition. Also, to 
analyze final state hadron spectra in the energy variable ln(l/x) and to model Bose-Einstein correlations. 

MSEL : determines which action is to be taken, with the convention noted before. 

= 10 : statistics on time evolution of particle numbers, energy partitions, longitudinal momentum 
fractions and average transverse momentum is reset. 

= 11 : the particle content is analyzed in each single time step during an event, for the different 
particle species individually. 

— 12 : gives a table on the obtained time evolution profiles, written to file VNI210.DAT. 

= 20 : statistics on y = ln{l/x) distribution of final hadrons is reset, where x = 2E/yfs is the 
energy fraction of a particle. 

= 21 : hadron content at the end of each event is binned in y = ln(l/,x) individually for all, and for 
charged, hadrons. Covered range in y is given by PARW(161) and PARW(162). 

— 22 : gives a table of the obtained y = ln{l/x) spectra, written to file VNI220.DAT. 

= 30 : statistics on Bose-Einstein correlations of final state pious and kaons is reset. 
= 31 : Bose-Einstein effects are simulated according to a simple parametrization, with an algorithm 
adopted from PYTHIA. The distribution of pair masses is stored for both cases, without and with 
Bose-Einstein correlations taken into account. Covered range in pair mass is given by PARW(163) 
and PARW(164). 

= 32 : gives a table of the obtained pair mass spectrum of 'like-sign' pions and kaons, without and 
with Bose-Einstein effects, written to file VNI230.DAT. 



SUBROUTINE VNIANA3(MSEL) 

Purpose: to analyze the dynamics of pre-hadronic clusters formed from partons during the parton-hadron 

conversion stage of the system. The cluster size- and mass-distributions, and momentum and spatial spectra 
are accumulated at various time steps during the time evolution of the system, so that the time change of 
the spectra gives a history of the cluster formation processes, locally in phase-space. 

MSEL : determines which action is to be taken, with the convention noted before. The statistics accu- 
mulated is output in tabular form with bins, ideally covering the range of momenta and coordinates 
of relevant particles. Depending on the kinematics, however, the limits of this range may require an 
adjustment, and/or a finer resolution of bins may be necessary. 
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= 10 : statistics on cluster sizes and their invariant mass spectrum is reset. 

= 11 : the cluster content of the system is analyzed each single time step during an an event, and 
cluster sizes and masses are stored in bins. Covered range of analysis is < < 1.5 fm, and < M < 7 
GeV. 

= 12 : gives a table on the obtained cluster size- and mass-spectra, written to file VNI310.DAT. 
= 20 : statistics on rapidity (y) spectra of clusters is reset. 

= 21 : the cluster content at certain pre-selected time steps during an an event is analyzed. Covered 
range in y is given by PARW(171) and PARW(172). 

= 22 : gives a table on the obtained y-spectra, written to file VNI320.DAT. 
= 30 : statistics on longitudinal (r^) spatial distribitions of clusters is reset. 

= 31 : the cluster positions at certain pre-selected time steps during an event are analyzed, with 
respect to the overall CM frame and to the z-axis, usually the jet-axis, or beam-axis. Covered range 
in V;, is given by PARW(173) and PARW(174). 

= 32 : gives a table on the obtained rz-spectra, written to file VNI330.DAT. 

= 40 : statistics on transverse momentum {p±) distribitions of clusters is reset. 

= 41 : the cluster momenta p± at certain pre-selected time steps during an an event are analyzed, 
transverse with respect to the z-axis. Covered range in p±_ is given by PARW(175) and PARW(176). 

— 42 : gives a table on the obtained p^-spectra, written to file VNI340.DAT. 

— 50 : statistics on transverse {r±) spatial distribitions of clusters is reset. 

= 51 : the cluster content at certain pre-selected time steps during an an event is analyzed, with 
respect to the distance transverse to the z-axis. Covered range in r± is given by PARW(177) and 
PARW(178). 

— 52 : gives a table on the obtained r_L-spectra, written to file VNI350.DAT. 



SUBROUTINE VNIANA4(MSEL) 

Purpose: to analyze momentum and spatial distributions of partons and hadrons at various time steps during 
the time evolution of the system, so that the time change of the spectra gives a history of the dynamics 
from the initial state all the way to the final state of hadrons. 

MSEL : determines which action is to be taken, with the convention noted before. The statistics accu- 
mulated is output in tabular form with bins, ideally covering the range of momenta and coordinates 
of relevant particles. Depending on the kinematics, however, the limits of this range may require an 
extension. 

= 10 : statistics on rapidity (y) spectra of partons and hadrons is reset. 

= 11 : the particle content at certain pre-selected time steps during an an event is analyzed, for gluons, 
quarks, mesons and baryons seperately. Covered range in y is given by PARW(181) and PARW(182). 
= 12 : gives a table on the obtained y-spectra, written to file VNI410.DAT. 

= 20 : statistics on longitudinal (r^) spatial distribitions of partons and hadrons is reset. 
= 21 : the particle positions at certain pre-selected time steps during an event are analyzed, with 
respect to the overall CM frame and to the z-axis, usually the jet-axis, or beam-axis. Covered range 
in is given by PARW(183) and PARW(184). 

= 22 : gives a table on the obtained r^-spectra, written to file VNI420.DAT. 

= 30 : statistics on transverse momentum {p±) distribitions of partons and hadrons is reset. 
= 31 : the particle momenta p± at certain pre-selected time steps during an an event are analyzed, 
transverse with respect to the z-axis. Covered range in p±_ is given by PARW(185) and PARW(186). 
= 32 : gives a table on the obtained pT-spectra, written to file VNI430.DAT. 
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= 40 : statistics on transverse {r±) spatial distribitions of partons and hadrons is reset. 

= 41 : the particle content at certain pre-selected time steps during an an event is analyzed, with 
respect to the distance transverse to the z-axis. Covered range in r± is given by PARW(187) and 
PARW(188). 

— 42 : gives a table on the obtained r^-spectra, written to file VNI440.DAT. 



SUBROUTINE VNISPHE(SPH, APL) 

Purpose: to diagonalize the momentum tensor, i.e. find the eigenvalues Ai > A2 > A3, with sum unity, and 

the corresponding eigenvectors. Momentum power dependence is given by PARU(41); default corresponds 
to sphericity, PARU(41)=1. gives measures linear in momenta. Which particles (or partons) are used in the 
analysis is determined by the MSTU(41) value. 

SPH : 3(A2 + A3)/2, i.e. sphericity (for PARU(41)=2.). Is returned = -1, if analysis not performed because 
event contained less than two particles (or two exactly back-to-back particles, in which case the two 
transverse directions would be undefined). 

APL : 3/2A3, i.e. aplanarity (for PARU(41)=2.). Is returned = -1, as SPH = -1 

Remark: the fines N + 1 through iV + 3 (iV - 2 through N for MSTU(43) = 2) in VNIREC1/VNIREC2 
will, after a call, contain the following information: 

K{N + i,l) = 31; 
K(iV + i,2) = 95; 

K(iV + i, 3) : i, the axis number, i = 1, 2, 3; 
K(N + i,4), K(iV + i,5) =0; 

P(A^ + i, 1) - F{N + i, 3) : the i-th eigenvector, x,y and z components; 

P(A^ + ^,4) : Aj, the i-th eigenvalue; 

P(A^ + i,5) = 0; 

\{N + i,l) - Y{N + i,5) = 0. 

Also, the number of particles used in the analysis is given in MSTU(62). 



SUBROUTINE VNITHRU (THR , OBL) 

Purpose: to find the thrust, major and minor axes and corresponding projected momentum quantities, in 
particular thrust and oblateness. The performance of the program is affected by MSTU(44), MSTU(45), 
PARU(42) and PARU(48). In particular, PARU(42) gives the momentum dependence, with the default value 
1. corresponding to linear dependence. Which particles (or partons) are used in the analysis is determined 
by the MSTU(41) value. 

THR : thrust (for PARU(42)=1.). Is returned = -1, if analysis cannot be not performed, because event 
contained less than two particles. Is returned = -2, if remaining space in VNIREC1/VNIREC2 (partly 
used as working area) is not large enough to allow analysis. 

OBL : oblateness (for PARU(42)=1.). Is returned = -1, respectively = -2, as for THR. 

Remark: the lines + 1 through A^ + 3 (A^ - 2 through N for MSTU(43) = 2) in VNIREC1/VNIREC2 
will, after a call, contain the following information: 

K{N + i,l) = 31; 
K(Ar + i,2) = 96; 

K(A^ + ^, 3) : i, the axis number, i = 1, 2, 3; 
K{N + i,4), K(Ar + i,5) =0; 

P(A' + i, 1) - P(A^ + ^, 3) : the thrust, major and minor axis, respectively, for i = 1, 2 and 3; 
P(A^ + ^,4) : corresponding thrust, major and minor value; 

P{N + i,5) =0; 

V(A^ + i,l) - Y{N + i,5) = 0. 

Also, the number of particles used in the analysis is given in MSTU(62). 
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SUBROUTINE VNICLUS (N JET) 



Purpose: to reconstruct an arbitrary number of jets using a cluster analysis method based on particle 
momenta. Two different distance measures are available, the traditional one roughly corresponding to 
relative transverse momentum, and one based on the JADE method, which roughly corresponds to invariant 
mass (in both cases with some important modifications). The choice is controlled by MSTU(46). The 
distance scale djcnn, above which two clusters may not be joined, is normally given by PARU(44). In general, 
djoin may be varied to describe different "jet resolution powers". With the mass distance, PARU(44) can 
be used to set the absolute maximum cluster mass, or PARU(45) to set the scaled one, i.e. in y = m^/W"^, 
where is the total invariant mass-squared of the particles being considered. It is possible to continue the 
cluster search from the configuration already found, with a new higher djoin scale, by selecting MSTU(48) 
properly. In MSTU(47) one can also require a minimum number of jets to be reconstructed; combined with 
an artificially large djoin this can be used to reconstruct a predetermined number of jets. Which particles 
(or partons) are used in the analysis is determined by the MSTU(41) value, whereas assumptions about 
particle masses are given by MSTU(42). The parameters PARU(43) and PARU(48) regulate more technical 
details (for events at high energies and large multiplicities, however, the choice of a larger PARU(43) may 
be necessary to obtain reasonable reconstruction times). 

NJET : the number of clusters reconstructed. Is returned = -1, if analysis cannot performed, because event 
contained less than MSTU(47) (normally 1) particles, or analysis failed to reconstruct the requested 
number of jets. Is returned = -2, if remaining space in VNIREC1/VNIREC2 (partly used as working 
area) is not large enough to allow analysis. 

Remark: if the analysis does not fail, further information is found in MSTU(61) - MSTU(63) and PARU(61) 
- PARU(63). In particular, PARU(61) contains the invariant mass for the system analyzed, i.e. the 
number used in determining the denominator of y = m? /W"^ . PARU(62) gives the generalized thrust, 
i.e. the sum of (absolute values of) cluster momenta divided by the sum of particle momenta (roughly 
the same as multicity). PARU(63) gives the minimum distance d (in or m) between two clusters 
in the final cluster configuration, in case of only one cluster. Further, the lines iV + 1 through 
N + NJET [N - NJET+l through N for MSTU(43) = 2) in VNIREC1/VNIREC2 will, after a cah, 
contain the following information: 
K(iV + i, 1) = 31; 
K(iV + i,2) = 97; 

K(iV + i, 3) : i, the jet number, with the jets arranged in falling order of absolute momentum; 
K{N + i, 4) : the number of particles assigned to jet i; 

K(iV + f,5) = 0; 

P(A^ + - P{N + i, 5) : momentum, energy and invariant mass of jet i; 
\{N + i,l) - V(iV + i,5) =0. 

Also, for a particle which was used in the analysis, K(/,4) = i, where / is the particle number and i 
the number of the jet it has ben assigned to. Undecaycd particles not used then have K(1, 4) = 0. An 
exception is made for lines with K(/, 1) = 3 (which anyhow are not normally interesting for cluster 
search), where the color flow information stored in K(I,4) is left intact. 



SUBROUTINE VNICELL(NJET) 

Purpose: to provide a simpler cluster routine more in line with what is currently used in the study of high-p± 
collider events. A detector is assumed to stretch in pscudorapidity between -PARU(51) and +PARU(51) 
and be segmented in MSTU(51) equally large rj (pscudorapidity) bins and MSTU(52) (p (azimuthal) bins. 
Transverse energy E_i for undecayed entries are summed up in each bin. For MSTU(53) nonzero, the energy 
is smeared by calorimetric resolution effects, cell by cell. This is done according to a Gaussian distribution; 
if MSTU(53) = 1 the standard deviation for the E± is PARU(55)V^, if MSTU(53) = 2 the standard 
deviation for the E is PARU(55)v^, E± and E expressed in GeV. The Gaussian is cut off at and at 
a factor PARU(56) times the correct E± or E. All bins with E± >PARU(52) are taken to be possible 
initiators of jets, and arc tried in falling E2_ sequence to check whether the total E± summed over cells no 
more distant than PARU(54) in ^ (Ary)^ + (A^)^ exceeds PARU(53). If so, these cells define one jet, and 
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axe removed from further consideration. Contrary to VNICLUS, not all particles need be assigned to jets. 
Which particles (or partons) are used in the analysis is determined by the MSTU(41) value. 

NJET : the number of jets reconstructed (may be 0). Is returned = -2, if remaining space in 
VNIREC1/VNIREC2 (partly used as working area) is not large enough to allow analysis. 

Remark: the Hncs + 1 through N + NJET {N - NJET + 1 through N for MSTU(43) = 2) in 
VNIREC1/VNIREC2 will, after a call, contain the following information: 

K{N + = 31; 
K(A^ + i,2) = 98; 

K(A^ + ^, 3) : i, the jet number, with the jets arranged in falling order in E_\_\ 

K{N + i, 4) : the number of particles assigned to jet i; 

K{N + i,5) = 0; 

Y{N + i,l) - Y{N + i,5) =0; 

Further, for MSTU(54) = 1: 

P(A^ + i, 1), P{N + i,2) = position in rj and cj) of the center of the jet initiator cell, i.e. geometrical 
center of jet; 

P(iV + i,3), P(N+i,4) = position in rj and cf) of the £^_L-weighted center of the jet, i.e. the center of 

gravity of the jet; 

P(iV + i, 5) = E -E± of the jet; 

while for MSTU(54) = 2: 

P(A'^ + - P{N + i, 5) : the jet momentum vector, constructed from the summed E± and the rj and 
4> of the £^_L -weighted center of the jet as {px,Py,Pz, E,m) = £^_l(cos^, sin^, sinh?7,cosh?7, 0); 

and for MSTU(54) = 3: 

P(A'^ + 1) - P{N + i,5) : the jet momentum vector, constructed by adding vectorially the momentum 
of each cell assigned to the jet, assuming that all the E± was deposited at the center of the cell, and 
with the jet mass in P(A'' + i, 5) calculated from the summed E and p as m? = E"^ — p^ — p^ — p^. 

Also, the number of particles used in the analysis is given in MSTU(62), and the number of cells 
hit in MSTU(63). 



SUBROUTINE VNIJMAS (PMH , PML) 

Purpose: to reconstruct high and low jet mass of an event, by dividing the collision event into two hemi- 
spheres, transversely to the sphericity axis. Then one particle at a time is reassigned to the other hemisphere 
if that reduces the sum of the two jet masses- squared rrijj + m\. The procedure is stopped when no further 
significant change (see PARU(48)) is obtained. Often, the original assignment is retained as it is. Which 
particles (or partons) are used in the analysis is determined by the MSTU(41) value, whereas assumptions 
about particle masses are given by MSTU(42). 

PMH : heavy jet mass (in GeV). Is returned = -2m, if remaining space in VNIREC1/VNIREC2 (partly 
used as working area) is not large enough to allow analysis. 

PML : light jet mass (in GcV). Is returned = -2, if remaining space in VNIREC1/VNIREC2 (partly used 
as working area) is not large enough to allow analysis. 

Remark: light jet mass (in GeV). After a successful call, MSTU(62) contains the number of particles used 
in the analysis, and PARU(61) the invariant mass of the system analyzed. The latter number is helpful 
in constructing scaled jet masses. 



SUBROUTINE VNIFOWO (HIO , H20 , H30 , H40) 
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Purpose: to do an event analysis in terms of the Fox- Wolfram moments The moments Hi are normalized to 
the lowest one, Hq. Which particles (or partons) are used in the analysis is determined by the MSTU(41) 
value. 



if 10 : H\/Hq. Is = if momentum is balanced, 
if 20 : H2/HQ. 
H30 : H3/H0. 
H40 : Hi/Ho. 

Remark: the number of particles used in the analysis is given in MSTU(62). 



81 



C Further commonblocks 



The subroutines and commonblocks that one normally comes in direct contact with have already been 
described. For the sake of completeness, some further commonblocks are here described, which might be of 
interest to a user that desires a more detailed insight. 

The common block VNISYS contains information concerning the type and status of the collision system 
A + B defined by the beam and target particles, c.f Sec. 3.4. 

COMMON/VNISYS/NSYS (2,5), KSYS (2,5), PSYS (2,5), RSYS (2 , 5) , AB (500 , : 10) 

Purpose: to store information and data concerning beam and target particles. It is in particular of interest 
if one desires to use the user-specified frames USER, FOUR or FIVE, in which case the user has to give 
the initial momenta PSYS and coordinates RSYS of the colliding particles by assigning values before the 
initialization. Furthermore, for collisions involving nuclei, the array AB provides information on how many 
initial-state nucleons have been 'wounded' and how often they have been struck. For collisions without 
nuclei, the array AB is empty! 

NSYS(I,1) : number of neutrons in beam/target particle / = 1,2. 
NSYS(I,2) : number of protons. 
NSYS(I,3) : number of d valence quarks. 
NSYS(I,4) : number of u valence quarks. 
NSYS(I,5) : number of s valence quarks. 

KSYS(I,1)- KSYS(I,5) : status code, particle id, origin of beam/target particle / = 1,2, in analogy to 
the array K in the particle record VNIREC1/VNIREC2. 

PSYS(I,1)- PSYS(I,5) : initial momentum, energy, mass of beam/target particle, in analogy to the 
array P in the particle record VNIREC1/VNIREC2. 

RSYS(I,1)- RSYS(I,5) : initial x, y, z position before collision of of beam/target particle, in analogy to 
the array R in the particle record VNIREC1/VNIREC2. (RSYS(I,4) and RSYS(I,5) are currently not 
used). 

AB(IH,0) : number of collisions of nucleon IH, that is, the number of times an initial-state parton 
belonging to that nucleon gets struck and liberated from its mother nucleon. Here 1 < IH < NnuciA) + 
Nnuc{B), with A and/or B being the initial beam and/or target nucleus. AB(/iJ,0) = initially for 
all IH , and at the end allows to discriminate 'wounded' (> l)and 'unwounded'(= 0) nucleons. 

AB(IH,1) : initial number of partons in nucleon IH, before collision. 

AB(IH,2) : current number of partons in nucleon IH, during time evolution. For 'unwounded' nucleons 
AB(/i^,2)=AB(/i7,l). 

AB(IH,3) : KF particle id of nucleon IH (either proton or neutron). 

AB(IH,4) : KF particle id of mother nucleus to which nucleon IH belongs. 

AB(IH,5) : initial P± of nucleon IH (before collision) with respect to collision axis of beam/target. 

AB(IH,6) : initial Pz of nucleon IH with respect to collision axis. 

AB(IH,7) : initial E of nucleon IH with respect to chosen Lorentz frame. 
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AB(IH,8) : current P± of nucleon IH (during time evolution of collision) with respect to collision axis of 
beam/target. 

AB(IH,9) : current Pz of nucleon IH with respect to collision axis. 
AB(IH,10): current E of nucleon IH with respect to chosen Lorentz frame. 

The common block VNINIA contains detailed information on the number of parton interactions for each of 
the included sub-process channels. 

C0MM0N/VNINIA/NIA(5,0:10) 

Purpose: to store accumulated statistics on type and number of parton interactions. 

NIA(1,0) : total number of hard parton collisions a + h ^ c + d. 

NIA(1,1) : number of hard q + q ^ q + q. 
NIA(1,2) : number of hard q + q ^ q + q. 
NIA(1,3) : number of hard q + q ^ g + g- 
NIA(1,4) : number of hard q + q ^ g + ^■ 
NIA(1,5) : number of hard g + g — 7 + 7. 
NIA(1,6) : number of hard q + g ^ q + g- 
NIA(1,7) : number of hard q + g ^ q + ^■ 
NIA(1,8) : number of hard g + g q + q- 
NIA(1,9) : number of hard g + g ^ g + g- 

NIA(2,0) : total number of soft parton collisions a + b —>■ c + d. 

NIA(2,1) : number of soft q + q ^ q + q- 

NIA(2,2) : number of soft q + q ^ q + q- 

NIA(2,3) : number of soft q + q ^ g + g. 

NIA(2,4) : number of soft q + q ^ 5 + 7- 

NIA(2,5) : number of soft + g ^ 7 + 7. 

NIA(2,6) : number of soft q + g ^ q + g. 

NIA(2,7) : number of soft q + g ^ q + ^■ 

NIA(2,8) : number of soft g + g ^ q + q- 

NIA(2,9) : number of soft g + g ^ g + g- 

NIA(3,0) : total number of parton fusions a + 6 — c*. 

NIA(3,1) : number of g + q — > 5*. 
NIA(3,2) : number oiq + g^q*. 
NIA(3,3) : number oi g + g ^ g* . 

NIA(4,0) : total number of space- like branchings a ^ b + c. 

NIA(4,1) : number oi q ^ q + g. 
NIA(4,2) : number oi g ^ g + g. 
NIA(4,3) : number oi g ^ q + q. 
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NIA(4,4) : number of g — > g + 7. 

NIA(5,0) : total number of time-like branchings a — > 6 + c. 

NIA(5,1) : number of q ^ q + g. 

NIA(5,2) : number of g ^ g + g. 

NIA(5,3) : number of g ^ q + q. 

NIA(5,4) : number of g — > g + 7. 



The common block VNITIM contains the (upon initialization) pretabulated time increments and time slizes 
used for following the time-evolution of an collision event. 

CDMMDN/VNITIM/TINC (0 : 2500) , TIME (0 : 2500) 

Purpose: to store time increments TINC and time steps TIME (in 1/GeV). 

TINC(I) : i-th time increment between time steps i — 1 and i, where i = MSTW(12) and 1 < i < 
MSTW(2). PARW(l) is assigned that increment value at a each time step during a collision event. 

TIME(I) : i-th time step TIME(i) = TIME(i - 1) + TINC(i), where i = MSTW(12) and 1 < i < 
MSTW(2). PARW(2) is assigned that time value at a each time step during a collision event. 



The common block VNIINT serves as a temporary storage of quantities during the simulation of a 
specific parton collision plus its associated branching evolution. It is very similar to a corresponding array 
in the PYTHIA program. 

C0MM0N/VNIINT/MINT(400) , VINT (400) 

Purpose: to provide an intermediate storage for integer and real valued variables as work area, used in the 
program to carry out the evolution of a single partonic subsystem at a time, consisting of a parton collision 
and its associated space-like (initial state) and time-like (final state) branchings. These variables must not 
be changed by the user. 

MINT(l) : specifies the general type of partonic subprocess that has occured, according to the ISUB code 
given in Table 3. 

MINT(2) : whenever MINT(l) (together with MINT(15) and MINT(16)) are not enough to specify the 
type of a 2 ^ 2 parton collision uniquely, MINT(2) provides an ordering of the different possibilities. 



related to different color-flow topologies. The different possibilities are discussed in Refs. |57]. With i 
-- MINT(15), j = MINT(16) and k = MINT(2), they are: 

ISUB = 1, i= j, qiqi q^qi; 
k = 1 : color configuration A. 
k = 2 : color configuration B. 

ISUB = 1, i = j, qiqj^ q^qj; 
k = 1 : only possibility. 

ISUB = 2, qiqi qiqr, 
k = 1 : only possibility. 

ISUB = 3, qiqi gg; 

k = 1 : color configuration A. 
k = 2 : color configuration B. 
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ISUB = 6, Qig Qig; 

k = 1 : color configuration A. 
k = 2 : color configuration B. 

ISUB = 8, gg ^ qiqi; 

k = 1 : color configuration A. 
k = 2 : color configuration B. 

ISUB = 9,gg^ gg; 



k = l 

k = 2 
k = 3 



color configuration A. 
colour configuration B. 
color configuration C. 



MINT (3) : number of partons produced in the hard interactions, i.e. the number n of the 2 — ^ n matrix 
elements used, where currently only n = 1 or n = 2 processes are included. 

MINT (4) : number of documentation lines in commonblocks VNIREC1/VNIREC2 that are reserved and 
used as work space for simulating a parton collision. 

MINT(5) : not used. 

MINT(6) : current frame of event as initialized, cf. MSTW(4). 

MINT (7) : line number for documentation of outgoing particle 1 from a parton collision 2 — > 2 or 2 — ^ 1. 

MINT (8) : line number for documentation of outgoing particle 2 from a parton collision 2—^2. Is for 
2^1. 

MINT(9) : not used. 

MINT(IO) : not used. 

MINT(ll) : KF flavour code of initial beam (side 1) particle (cf. MSTW(5)). 
MINT(12) : KF flavour code of initial target (side 1) particle (cf. MSTW(6)). 

MINT (13) : KF flavour code of side 1 parton that initiates a space- like cascade (only for primary partons, 
is else). 

MINT(14) : KF flavour code of side 2 parton that initiates a space-like cascade (only for primary partons, 
is else). 

MINT(15) : KF flavour code for side 1 incoming parton to collision subprocess. 
MINT(16) : KF flavour code for side 2 incoming parton to collision subprocess. 

MINT(17) - MINT(20) : not used. 

MINT(21) : KF flavour code for parton 3 outgoing from a parton collision. 

MINT(22) : KF flavour code for parton 4 outgoing from a parton collision (is for 2 1 processes). 

MINT(23) : KF flavour code for parton 5 outgoing from a parton collision - not used currently (relevant 
only for 2 — 3, 2 — 4 processes). 

MINT(24) : KF flavour code for parton 6 outgoing from a parton collision - not used currently (relevant 
only for 2 — 4 processes). 

MINT(25) - MINT(30) : not used. 
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MINT(31) - MINT(40) : not used. 



MINT(41) : type of initial beam particle; 1 for lepton and 2 for hadron and 3 for nucleus (c.f. MSTW(7)). 

MINT(42) : type of initial target particle; 1 for lepton and 2 for hadron and 3 for nucleus (c.f. MSTW(8)). 

MINT(43) : combination of incoming beam and target particles (c.f. MSTW(9)). 

MINT (44) : as MINT (43), but a photon counts as a lepton. 

MINT(45), MINT(46) : structure of incoming beam and target particles. 
= 1 : no internal structure, i.e. lepton carries full beam energy. 
— 2 : defined with structure functions, i.e. hadrons and nuclei. 

MINT (47) : combination of incoming beam and target particle structure function types. 
= 1 : no structure function either for beam or target. 
= 2 : structure functions for target but not for beam. 
= 3 : structure functions for beam but not for target. 
= 4 : structure functions both for beam and target. 

MINT (48) : total number of subprocesses switched on. 

MINT(49), MINT(50) : not used. 

MINT(51) : internal flag that the simulation of a 2-parton collision system, including space-like initial 
state and time-like final state cascades, has failed cuts. 
= : no problem. 
= 1 : failed and to be skipped. 

MINT(52) - MINT(60) : not used. 
MINT(61) - MINT(70) : not used. 
MINT(81), MINT(82) : not used. 

MINT(83) : number of lines in the event record already filled by previously considered pileup events. 

MINT (84) : MINT (83) + MSTP(126), i.e. number of lines already filled by previously considered events 
plus number of lines to be kept free for event documentation. 



VINT(l) : EcM, CM energy. 

VINT(2) : s = E^j^, mass-square of complete system. 

VINT (3) : mass of beam particle. 

VINT (4) : mass of target particle. 

VINT (5) : momentum of beam (and target) particle in CM frame. 

VINT (6) : angle for rotation from CM frame to fixed target frame. 

VINT(7) : angle $ for rotation from CM frame to fixed target frame. 

VINT(8) : Px velocity for boost from CM frame to fixed target frame. 

VINT (9) : Py velocity for boost from CM frame to fixed target frame. 
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VINT(IO) : f3z velocity for boost from CM frame to fixed target frame. 

VINT(ll) : Px velocity for boost from overall CM frame to 'microscopic' cm frame of a parton-parton 
collision subsystem. 

VINT(12) : (3y velocity for boost from overall CM frame to parton-parton cm frame. 

VINT(13) : $z velocity for boost from overall CM frame to parton-parton cm frame. 

VINT (14) : angle for rotation from overall CM frame to 'microscopic' cm frame of a parton-parton 
collision subsystem. 

VINT(15) : angle (/) for rotation from overall CM frame to parton-parton cm frame. 

VINT(16) : Lorentz boost factor 7 to parton-parton collision cm frame. 

VINT (17) : Rapidity jjcm of 2-parton cm sytem with respect to global CM frame. 

VINT(18) : Transverse momentum p± cm of 2-parton cm sytem in global CM frame. 

VINT(19) : Longitudinal separation Afz of 2 partons in their cm. frame of collision, at the point of 
'closest approach'. 

VINT (20) : Transverse separation Af± of 2 partons in their cm. frame of collision, at the point of 'closest 
approach'. 

VINT(21) : r = X1X2S with xi^2 the longitudinal momentum fractions (Bjorken-x) of 2 colliding partons. 
VINT(22) : y* = 1/2 ln(xi/a;2), the 'relative' rapidity of 2 colliding partons. 
VINT(23) : cos(^)) (scattering angle) in 2 — 2 parton collision. 
VINT (24) : (f) (azimuthal angle) in in 2 2 parton collision. 

VINT(25) - VINT(30) : not used. 

VINT (31) : ti, lower value of kinematic limit on t in parton collision. 
VINT (32) : tu, upper value of kinematic limit on t in parton collision. 

VINT(33) : upper absolute value of limit of integration over dt of 2-parton cross-section. 
VINT(34) : \t\2, lower absolute value of limit of integration over dt of 2-parton cross-section. 
VINT(35) : cos{§)max for ^ < 0. 
VINT(36) : cos(^)„a^ for ^ > 0. 

VINT(37) : integrated 2-parton cross-section a{s) = J dp\{da{s,p\) / dp]_. 
VINT(38) : differential 2-parton cross-section da{s,p\)/dp'j_. 

VINT(39) : characteristic proper collision time Tcol{p\) oc l/p± in the 2-body c.m.-frame. 
VINT (40), not used. 

VINT (41) : the momentum fraction xi taken by parton 1 at the 2-parton vertex. 
VINT (42) : the momentum fraction X2 taken by parton 2 at the 2-parton vertex. 
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VINT(43) : m = \/?, mass of 2-parton collision subsystem. 
VINT (44) : s of the collision subprocess (2 ^ 2, 2 1). 
VINT(45) : t of the collision subprocess (2 ^ 2, 2 ^ 1). 
VINT (46) : u of the collision subprocess (2 ^ 2, 2 ^ 1). 

VINT (47) : p± of the collision subprocess, i.e., transverse momentum evaluated in the rest frame of the 
collision. 

VINT(48) : p]_ of the collision subprocess; see VINT(47). 

VINT(49) : p±min of process, i.e. PARV(42), or max(VKIN(3),VKIN(5)), depending on which is larger, 
and whether the process is singular in — > 0, or not. 

VINT (50) : not used. 

VINT (51) : Q of the hard subprocess. The precise definition is process-dependent, see MSTV(43). 
VINT(52) : ^^e hard subprocess; see VINT(51). 

VINT(53) : Q of the outer hard-scattering subprocess. Agrees with VINT(51) for a 2 ^ 1 or 2 ^ 2 
process. 

VINT(54) : of the outer hard-scattering subprocess; see VINT(53). 

VINT(55) : Q scale used as maximum virtuality in parton showers. Is equal to VINT(53), except for 
deep-inelastic-scattering processes . 



VINT(56) 
VINT (57) 
VINT (58) 

VINT(59) 
ing. 



scale in parton showers; see VINT(55). 
cnem value of hard process. 
as value of hard process. 

sin^ (cf. VINT(23)); used for improved numerical precision in elastic and diffractive scatter- 



VINT(61), VINT(62) : nominal values, i.e. without initial-state radiation effects, for the two 
partons entering the hard interaction. 

VINT(63), VINT(64) : nominal values, i.e. without final-state radiation effects, for the two (or 
one) partons/particles leaving the hard interaction. 

VINT(65) : Pinit, i-e. common nominal absolute momentum of the two partons entering the hard inter- 
action, in their rest frame. 

VINT(66) : Pf m, i c. common nominal absolute momentum of the two partons leaving the hard interac- 
tion, in their rest frame. 

VINT(67) - VINT(IOO) : not used. 

VINT(lOl) - VINT(105) : 6, (p and (3 for rotation and boost from cm. frame to hadronic cm. frame 
of a lepton-hadron or lepton- nucleus event. 

VINT(106) - VINT(llO) : 6, (p and /3 for rotation and boost from cm. frame to frame of equal 
nucleon-nucleon velocity in the case of a hadron-nucleus or an asymmetric nucleus-nucleus event with 
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The common block VNICOL contains accumulated data on parton collisions during a beam/target event. 

COMMON/VNICOL/MCOL (1000,0:10), VCOL (1000,0:10) 

Purpose: to provide an intermediate storage for a minimal set of integer and real valued variables for all 
detected parton collisions that occur during a single timestep during an event. The first index I in the 
arrays MCOL and VCOL labels the sequence number of a collision and the second index J is a pointer to 
interface with the above MINT and VINT arrays wherein the calculated quantities are written to. 

MCOL(I,l) : line number in VNIREC1/VNIREC2 of parton 1 entering the collision subprocess. 

MCOL(I,2) : line number in VNIREC1/VNIREC2 of parton 2 entering the collision subprocess. 

MCOL(I,3) : = MINT(l) 

MCOL(I,4) : = MINT(2) 

MCOL(I,5) : = MINT(3) 

MCOL(I,6) : = MINT(15) 

MCOL(I,7) : = MINT(16) 

MCOL(I,8) : = MINT(21) 

MCOL(I,9) : = MINT(22) 

MCOL(I,10): not used. 



VCOL(I,l) : 
VCOL(I,2) : 
VCOL(I,3) : 
VCOL(I,4) : 
VCOL(I,5) : 
VCOL(I,6) : 
VCOL(I,7) : 
VCOL(I,8) : 
VCOL(I,9) : 
VCOL(I,10): 



= VINT (23 
= VINT(24 
= VINT(26 
= VINT(41 
= VINT(42 
= VINT (44 
= VINT (45 
= VINT (46 
= VINT(48 
= VINT(52 
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D Sample output 



1. General information on the simulation as-a- whole 

The following output stems simulating 500 pp collisions at CERN collider energy ^/s = 900 GeV. Unless 
switched of by setting MSTV{114:) = 0, this type of output information is provided automatically and 
written to the file VNIRUN.DAT at the end of a simulation. The output includes a summary of the 
initialization overall collision system, the gluon and quark content of the initial state, the type and number 
of occured interaction processes during the parton-cascade development, the number and type of produced 
pre-hadronic clusters from partons, and the numbers of resulting final-state hadrons. 
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Warning: this program is a Monte 
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Carlo program! So when you play 
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roulette, dont do it the Russian 
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way - avoid fatal results! 
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~ Copyright KayKayGee (1998) ~ 
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Depository of program 


versions and documentation: 


** 



** \protect\vrule widthOpt\protect\href {http: //rhic.phys . Columbia. edu/rhic/vni}{htt 

** ** 

** Official reference: ** 

** K. Geiger, Computer Physics Communications 104 (1997) 70. ** 

** (Further relevant references may be found therein) ** 

** ** 

** The Crew: ** 

** ** 

** KLAUS KINDER-GEIGER RON LONGACRE DINESH K. SRIVASTAVA ** 

** ** 

** Mailing address: ** 
** Physics Dept., Brookhaven National Laboratory, Upton, N.Y. 11973, U.S.A. ** 

** ** 

** Contact (e-mail / phone) : ** 

** klaus@bnl.gov / USA + (516) 344-3791 ** 

** longacreObnl.gov / USA + (516) 344-3713 ** 

** dks@vecdec.veccal.ernet.in / India + (33) 337-1231 ** 

** ** 
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******************* VNIXRIN : initialization of VNI routines ****************** 



1 1 

1 VNI will be initialized for a I 

I I 

I p+ on p- collider I 

1 1 

1 at 900.00 GeV center-of-mass energy I 

1 1 

1 corresponding to 900.00 GeV sqrt(s)/hadron(nucleon) I 

1 1 

1 h + h -> jets + X -> hadrons I 

1 1 

1 Number of collision events: 500 I 

I Number of time steps per event: 200 I 

I I 

I Initial time: .00 1/GeV = .00 fm = O.OOD+00 sec I 

I Final time: 176.10 1/GeV = 34.75 fm = 0.23D-22 sec 1 

1 1 

1 1 

1 Initial parton distributions from VNI : I 

I nucleon structure functions GRV set 1 I 

I pion structure functions DO set 1 I 

I I 

I Minimum initial resolution scale Q_0: 3.42 GeV I 

I Minimum sqrt(s) in parton collisions: 2.00 GeV I 

I Minimum pT_0 in parton collisions: 3.60 GeV I 

I Minimum q_0 in space-like branchings: 3.42 GeV I 

I Minimum mu_0 in time-like branchings: 1.48 GeV I 

I I 

I I 

I Initial Collision System (in global CM frame) : I 

I I 

I Beam Target I 

I Particle: p+ p- I 

I I 

I Number of neutrons : 1 

I Number of protons: 1 -1 I 

I Number of d_v quarks: 1 -1 I 

I Number of u_v quarks: 2-2 1 

I Number of s_v quarks: 1 

I Radius R_ms (f m) : .810 .810 I 

I I 

I Velocity: .999998 .999998 I 

I Gamma factor: 479.591 479.591 I 

I Rapidity: 6.866 -6.866 I 

I I 

I p_x (GeV) : .000 .000 I 

I p_y (GeV) : .000 .000 I 

I p_z (GeV): 449.999 -449.999 I 
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I E (GeV): 450.000 450.000 I 

I M (GeV): .938 .938 I 

I I 

I charge px_tot py_tot pz_tot E_tot M_tot I 

I .0 0. 0. 0. 900. 900. I 

I I 

I Initial separation parallel to beam: Delta Z = .25 fm I 

I Min/max impact transverse to beam: .00 < B < 1.19 fm I 

I I 

I I 



********************** VNIXRIN : initialization completed ********************* 
Welcome on board! Hope you survive this flight with VNI 



*********************** VNIXFIN : run finished normally *********************** 



I 
I 
I 


Total 


CPU time needed: 


hours 


21 minutes 


11 seconds 


I 
I 
I 






INITIAL 


STATE 











Internal parton structure functions selected: 

set 1(1) for nucleons (pions) at initial QO = 5.939 GeV 



Minimum resolvable x-value resulting from spread of confined 
partons with 1/xP < R_conf = .20 fm: 

Beam p+ : x_min = 0.222D-02 

Target p~- : x_min = 0.222D-02 

Average sum of momentum fractions per hadron (nucleon) : 

Beam p+ : <sum x> = .9853 

Target p~- : <sum x> = .9908 



Average number of initial state partons : 
gluons : 51.8160 
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u+u~+d+d~ : 
u-u~+d-d~ : 
s+s~ : 
c+c~ : 
b+b~ : 
t+t~ : 

total : 



49.4000 
.0000 
6.7800 
1.0520 
.0320 
.0000 

109.0800 



PARTON CASCADE EVOLUTION 



Total number of events N_evt : 500 

Total number of time steps N_tim: 200 
Initial/final time: .0 fm / 34.7 fm 
Analysis based on 500 event samples 



All niombers per non-dif f . event N_nd — N_nd/N_evt = 0.8880D+00 
Actual minimvim momentvim transfer for hard collisions: 3.6000 GeV 



hard 2 -> 2 collisions: 
soft 2 -> 2 collisions: 
fusions 2 -> 1 : 
space-like branchings 1 
time-like branchings 1 



switched on 
switched off 
switched on 

> 2: switched on 

> 2: switched on 



Average number of hard parton collisions a + b -> c + d : 



q+q 


-> 


q+q 


.0000 


q+q~ 


-> 


q+q~ 


.0000 


q+q~ 


-> 


g+g 


.0000 


q+q~ 


-> 


g+gamma 


.0000 


q+q~ 


-> 


2 gamma 


.0000 


q+g 


-> 


q+g 


.0135 


q+g 


-> 


q+gamma 


.0000 


g+g 


-> 


q+q" 


.0045 


g+g 


-> 


g+g 


.2095 


hard 


a+b 


-> c+d 


.2275 



- Average number of soft parton collisions a + b -> c + d : 



q+q 


-> 


q+q 


.0000 


q+q- 


-> 


q+q" 


.0000 


q+q~ 


-> 


g+g 


.0000 


q+q" 


-> 


g+gamma 


.0000 


q+q~ 


-> 


2 gamma 


.0000 


q+g 


-> 


q+g 


.0000 


q+g 


-> 


q+gamma 


.0000 


g+g 


-> 


q+q" 


.0000 
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g+g -> g+g : .0000 
total soft a+b -> c+d : .0000 



- Average number of parton fusions a + b -> c : 

q+q~ -> g* : .0000 

q+g -> q* : .0000 

g+g -> g* : .0135 

total fusion a+b -> c : .0135 



- Average number of space-like branchings a -> b + c: 



( 


CMshower -> b+c 


.1464 ) 




q -> q+g 


.0045 




g -> g+g 


.0045 




g -> q+q~ 


.0000 




q -> q+gamma 


.0000 


total 


sp-like a -> b+c 


.0090 



- Average number of time-like branchings a -> b + c: 



( 


CMshower -> b+c 


.8829 ) 




q -> q+g 


.0608 




g -> g+g 


.7050 




g -> q+q~ 


.0586 




q ~^ q+gamma 


.0000 


total 


tm-like a -> b+c 


.8243 



PARTON CLUSTER FORMATION 



NOTE: numbers refer to clusters from interacted partons only; 

beam/target remnant clusters (if present) are not included. 

Total number of events N_evt : 500 
Total number of time steps N_tim: 200 
Initial/final time: .0 fm / 34.7 fm 
Analysis based on 500 event samples 

All numbers per non-dif f . event N_nd — N_nd/N_tot = 0.8880D+00 



Total numbers of produced clusters: 
all: 1.6171 

fast: 1.3131 ( > PARV(88)= 1.000 MSTV(88)= ) 

slow: .3041 ( < PARV(88)= 1.000 MSTV(88)= ) 

exogamous: .0068 

Nvimber of elementary cluster formation processes 
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type 












total 


exogamous 


g + g 


-> 


C 


+ 


c 






.1712 


.0000 


g + g 


-> 


C 


+ 


g 






.4234 


.0000 


g + g 


-> 


c 


+ 


g 


+ 


g 


.7275 


.0000 


q + q' 


-> 


c 


+ 


c 






.0270 


.0000 


q + q- 


-> 


c 


+ 


g 






.0428 


.0023 


q + q- 


-> 


c 


+ 


g 


+ 


g 


.0000 


.0000 


g + q 


-> 


c 


+ 


q 






.0811 


.0045 


g + q 


-> 


c 


+ 


q 


+ 


g 


.1441 


.0000 


qq + q -> 


c 


+ 


c 






.0000 


.0000 


STom: 














1.6171 


.0068 



- Distribution of exogamous cluster processes: 



origin 




all 


fast 


slow 


0.00 - 0. 


.05 


.2658 


.2072 


.0586 


0.05 - 0, 


,20 


.0000 


.0000 


.0000 


0.20 - 0, 


.30 


.0000 


.0000 


.0000 


0.30 - 0, 


.45 


.0000 


.0000 


.0000 


0.45 - 0. 


.55 


.0023 


.0000 


.0023 


0.55 - 0, 


.70 


.0000 


.0000 


.0000 


0.70 - 0, 


.80 


.0023 


.0000 


.0023 


0.80 - 0. 


.95 


.0023 


.0023 


.0000 


0.95 - 1. 


.00 


1.3446 


1.1036 


.2410 



CLUSTER HADRONIZATION AND UNSTABLE PARTICLE DECAYS 



NOTE: numbers refer to both, clusters from interacted partons 
as well as beam/target remnant clusters (if present)DO 

Total number of events: 500 
Total number of time steps: 200 
Initial/final time: .0 fm / 34.7 fm 
Analysis based on 500 event samples 

parton-cluster fragmentation: switched on 

beam/target-cluster fragmentation: switched on 



- Average number of primary particles from cluster decays : 

type charged + neutral charged only 



leptons , gauge bosons : 




,0000 




.0000 


light mesons : 


18, 


,2700 


8 


.8700 


strange mesons : 


3, 


,8100 


1 


.7880 


charm, bottom mesons: 




,0440 




.0140 


tensor mesons : 


2, 


,2340 


1, 


.1220 


light baryons : 


3, 


,0220 


2 


.0860 


strange baryons : 


1, 


,4140 




.6960 


charm, bottom baryons: 




.0000 




.0000 


other: 




.0040 




.0040 
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total : 



49.2540 14.5800 



- Average number of secondary particles from unstable particle decays: 

type charged + neutral charged only 



leptons, gauge bosons: 


1 


9660 




0060 


licrht; mesons* 


51 


1120 


32 


4320 


strange mesons: 


5 


5240 


2 


8480 


charm, bottom mesons: 




0200 




0080 


tensor mesons: 




0060 




0040 


light baryons: 


5 


1600 


2 


5200 


strange baryons: 


1 


2200 




0060 


charm, bottom baryons: 




0000 




0000 


other : 




0000 




0000 


total : 


65 


0040 


37 


8220 



HADRON CASCADE EVOLUTION 



Total number of events : 500 
Total number of time steps : 200 
Initial/final time: .0 fm / 34.7 fm 
Analysis based on 500 event samples 

All numbers per non-dif f . event N_nd — N_nd/N_tot = 0.8880D+00 

Average number of cluster/hadron collisions a + b -> n (n >= 1) : 
C = prehadronic cluster 

Ml = low lying pseudoscalar & vector mesons 
M2 = excited & tensor mesons 
B = baryons 

channel all light strange charm bottom 



c + c 


-> 


n 





0000 





0000 





0000 





0000 





0000 


C + Ml 


-> 


n 





0000 





0000 





0000 





0000 





0000 


C + M2 


-> 


n 





0000 





0000 





0000 





0000 





0000 


C + B 


-> 


n 





0000 





0000 





0000 





0000 





0000 


Ml + Ml 


-> 


n 





0000 





0000 





0000 





0000 





0000 


Ml + M2 


-> 


n 





0000 





0000 





0000 





0000 





0000 


Ml + B 


-> 


n 





0000 





0000 





0000 





0000 





0000 


B + B 


-> 


n 





0000 





0000 





0000 





0000 





0000 


B + M2 


-> 


n 





0000 





0000 





0000 





0000 





0000 


M2 + M2 


-> 


n 





0000 





0000 





0000 





0000 





0000 


other unidentified 





0000 


















sum 









0000 





0000 





0000 





0000 





0000 
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Average kinematical quantities and average number of particles 
per cluster / hadron collision: 



sqrt(s) (GeV) sqrt(-t) (GeV) sqrt(-u) (GeV) # prod. part. 



light : 





0000 





0000 





0000 





0000 


strange : 





0000 





0000 





0000 





0000 


charm : 





0000 





0000 





0000 





0000 


bottom : 





0000 





0000 





0000 





0000 


all : 





0000 





0000 





0000 





0000 



Thanxxxxxxxxx 4 flying with VNI 
('We do know that you have the choice') 
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2. A typical event listing 



The event listing below is part of the standard output obtained for 500 pp collisions at CERN collider 
energy -y/s = 900 GeV. The listing, which is is written to the file VNIRUN.DAT, embodies the complete 
history of the first event from the initial beam particles p and p, their original gluon and quark content, 
the newly produced partons by hard scattering and radiation during the parton-cascade development, the 
formed pre-hadronic clusters and resulting final state particles. The latter emerge from both the parton 
clusters that result from parton-cascade products, and the beam/target clusters that consist of the remnant 
initial partons which have not participated in the parton cascades. As explained in Sec. 3.6, each entry 
contains name, status, flavor, and origin of a particle, as well as it's three-momentum, energy and mass. 
Particles in brackets are inactive, decayed or fragmented particles. The final line of the listing shows the 
total charge, three-momentum, anargy and invariant mass of the particle system, given by the sums over all 
active, stable particles. 

Particle listing (summary: momentum space P in GeV) 





event 


no. 


1 time step 


200 


at 


time 


34. 


75 fm 




I particle 


KS 


KF 


KMO 


P_x 


P. 


-y 


P. 


.z 


E 


M 


1 (p+) 


15 


2212 





.000 




000 


449 


999 


450 


000 


.938 


2 (p~-) 


15 


-2212 





.000 




000 


-449 


999 


450 


000 


.938 


3 (p+) 


15 


2212 


1 


.000 




000 


449 


999 


450 


000 


.938 


4 (p~-) 


15 


-2212 


2 


.000 




000 


-449 


999 


450 


000 


.938 


5 (u) 


13 


2 


3 


6.468 


- 


283 


67 


618 


67 


928 


.006 


6 (d) 


13 


1 


3 


-2.033 


7 


351 


87 


053 


87 


387 


.010 


7 (g) 


13 


21 


3 


. 135 


- 


110 


6 


530 


6 


530 


-.175 


8 (g) 


13 


21 


3 


.410 




205 


2 


324 


2 


324 


-.459 


9 (g) 


13 


21 


3 


-.009 




066 


3 


136 


3 


136 


-.067 


10 (g) 


13 


21 


3 


-.098 




153 


14 


160 


14 


160 


-.184 


11 (g) 


13 


21 


3 


-.094 




012 


2 


851 


2 


851 


-.095 


12 (g) 


13 


21 


3 


-.159 


- 


130 


4 


587 


4 


587 


-.205 


13 (g) 


13 


21 


3 


.745 


- 


058 


5 


911 


5 


910 


-.747 


14 (g) 


13 


21 


3 


.015 




208 


1 


636 


1 


636 


-.208 


15 (g) 


13 


21 


3 


-.002 




082 


1 


964 


1 


964 


-.082 


16 (g) 


13 


21 


3 


-.133 


- 


088 


1 


742 


1 


742 


-.159 


17 (g) 


13 


21 


3 


-.067 


- 


580 


4 


587 


4 


587 


-.584 


18 (g) 


13 


21 


3 


-.343 




022 


1 


636 


1 


636 


-.344 


19 (u~) 


13 


-2 


3 


-.355 




198 




988 




988 


-.406 


20 (u) 


13 


2 


3 


.132 




153 




988 




988 


-.202 


21 (g) 


13 


21 


3 


.007 




003 


2 


992 


2 


992 


-.010 


22 (g) 


13 


21 


3 


-.154 




440 


3 


746 


3 


746 


-.466 


23 (g) 


13 


21 


3 


-.038 




415 


17 


137 


17 


137 


-.418 


24 (g) 


13 


21 


3 


.054 




036 


2 


714 


2 


714 


-.065 


25 (g) 


13 


21 


3 


.054 




037 


1 


338 


1 


338 


-.065 


26 (u~) 


13 


-2 


3 


-.406 




103 


49 


084 


49 


083 


-.431 


27 (u) 


13 


2 


3 


-.268 




098 


1 


245 


1 


245 


-.285 


28 (g) 


13 


21 


3 


-.009 




030 


1 


742 


1 


742 


-.031 


29 (s~) 


13 


-3 


3 


.082 




062 


9 


842 


9 


842 


-.225 


30 (s) 


13 


3 


3 


.314 




241 


1 


070 


1 


070 


-.443 


31 (s~) 


13 


-3 


3 


-.008 




001 


7 


634 


7 


634 


-.200 


32 (s) 


13 


3 


3 


.016 




105 




988 




988 


-.226 


33 (d) 


13 


1 


3 


-.029 




227 


1 


851 


1 


851 


-.229 


34 (d~) 


13 


-1 


3 


-.030 




025 




988 




988 


-.040 


35 (g) 


13 


21 


3 


-.015 




008 


1 


636 


1 


636 


-.017 


36 (g) 


13 


21 


3 


.218 




061 


6 


744 


6 


744 


-.227 


37 (g) 


13 


21 


3 


.147 




042 


2 


581 


2 


581 


-.153 


38 (g) 


13 


21 


3 


.100 




031 


2 


450 


2 


450 


-.105 


39 (s~) 


13 


-3 


3 


-.080 




236 


1 


245 


1 


245 


-.319 


40 (s) 


13 


3 


3 


-.057 




207 




988 




988 


-.292 
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41 




13 


21 


3 




044 


-.095 


2 


992 


2 


992 




105 


42 


Cu; 


13 


2 


3 




343 


- . Oil 


1 


433 


1 


433 




343 


43 


Cu ; 


13 


-2 


3 




058 


. 067 




988 




988 




089 


44 




13 


21 


3 




084 


- . 195 


1 


245 


1 


245 




212 


45 




13 


21 


3 




048 


. 017 


1 


533 


1 


533 




051 


46 


(g) 


13 


21 


3 




008 


. 022 


2 


200 


2 


200 




024 


47 


(g) 


13 


21 


3 




019 


. 610 


4 


765 


4 


765 




611 


48 


(■g) 


13 


21 


3 




032 


- . 050 


2 


081 


2 


081 




060 


49 


Cu ; 


13 


-2 


3 




119 


- . 036 


1 


636 


1 


636 




125 


IT A 


Cu; 


lo 


Z 


o 
O 




-1/1/1 

144 


- . 129 




noo 

9oo 




noo 

9oo 




H no 

lyo 


51 


Ca; 


lo 


1 


o 
O 






- . 119 


o 
2 


CO H 

5ol 


o 
2 


CO H 

5ol 




oc c 

255 


52 


Cd ; 


13 


-1 


3 




209 


- . 049 




988 




988 




215 


53 


(g) 


13 


21 


3 




017 


. 334 


50 


278 


50 


278 




351 


54 


Cd ; 


13 


-1 


3 




127 


. 001 


1 


156 


1 


156 




127 


55 


(d) 


13 


1 


3 




182 


. 007 




988 




988 




182 


56 


(g) 


13 


21 


3 




026 


-.076 


18 


553 


18 


553 




089 


57 


(g) 


13 


21 


3 




095 


.035 


3 


746 


3 


746 




101 


58 


(g) 


13 


21 


3 




248 


- . 167 


4 


240 


4 


240 




299 


59 


(.g) 


13 


21 


3 




035 


. 128 


1 


338 


1 


338 




133 


60 


Cg) 


13 


21 


3 




170 


. 029 


6 


320 


6 


320 




173 


ol 


cgj 


lo 


21 


o 
O 




ceo 

552 


A ^ o 

.019 


A A 

14 


A TT 

'±77 


A A 

14 


A TT 

477 




ceo 

55o 


b2 


cg) 


lo 


^1 


o 
O 




■1 TO 

17o 


ooo 
- . 2o2 


o 
2 


00/1 

o24 


o 
2 


00/1 

o24 




O O /I 

oo4 


CD 

DO 


Cd ; 


lo 


-1 


o 
O 




210 


. 044 


o 
2 


OA A 

200 


o 
2 


OA A 

200 




o ^ /I 
214 


64 


Cd; 


13 


1 


3 




012 


- . 049 




988 




988 




051 


65 


Cu ; 


13 


-2 


3 




056 


- . 557 


2 


714 


2 


714 




560 


66 


Cu; 


13 


2 


3 




137 


-.062 




988 




988 




151 


67 


(g) 


13 


21 


3 




067 


- . 147 


2 


851 


2 


851 




161 


68 


(g) 


13 


21 


3 




409 


- . 114 


10 


915 


10 


915 




425 


69 


(g) 


13 


21 


3 




157 


- . 266 


6 


961 


6 


961 




310 


70 


Cg) 


13 


21 


3 




363 


- . 114 


5 


133 


5 


133 




380 


71 


(g) 


13 


21 


3 




110 


.285 


8 


338 


8 


338 




306 


72 


Cu; 


13 


2 


3 




093 


. 272 


6 


320 


6 


320 




288 


To 


Cu ; 


lo 




o 
O 




5d1 


. 05o 




ooo 

9oo 




noo 

9oo 




CCA 

5d4 


74 


Cs; 


13 


3 


3 




252 


- . 121 


7 


406 


7 


406 




344 


10 


Cs ; 


lo 


o 

-o 


o 
O 




ATO 

07o 


. OlD 




noo 

9oo 




noo 

9oo 




o ^ /I 
214 


76 


Cd; 


13 


1 


3 




095 


. Oil 


4 


072 


4 


072 




097 


77 


Cd ; 


13 


-1 


3 




013 


.005 




988 




988 




017 


78 


(g) 


13 


21 


3 




315 


. 023 


5 


711 


5 


711 




317 


79 


Cu ; 


13 


-2 


3 




376 


- . 047 


1 


245 


1 


245 




379 


80 


Cu) 


13 


2 


3 




339 


. 142 




988 




988 




368 


81 


(g) 


13 


21 


3 




009 


. 110 


1 


964 


1 


964 




111 


82 


eg) 


13 


21 


3 




143 


- . 143 


12 


337 


12 


337 




203 


83 


Cd ; 


13 


-1 


3 




279 


.015 


2 


200 


2 


200 




280 


84 


Cd; 


13 


1 


3 




093 


- . 020 




988 




988 




096 


O d 

ob 


cg) 


lo 


^1 


o 
O 




llo 


. lo9 


4 


d ov 

bo f 


4 


d ov 

bo 1 




^ nn 

19a 


oc 
OD 


Cu ; 


lo 




o 
O 




005 


- . 196 


■1 

1 


o /I c; 
24b 


■1 

1 


O /I d 

24b 




19d 




Cu; 


lo 


Z 


o 
O 




012 


- . 449 




AOO 

9oo 




noo 

9oo 




A A n 
449 


88 


Cu; 


13 


2 


3 




203 


. 190 


4 


947 


4 


947 




278 


89 


Cu ; 


13 


-2 


3 




086 


- . 351 




988 




988 




361 


90 


Cg) 


13 


21 


3 




041 


- . 083 


4 


240 


4 


240 




093 


91 


(g) 


13 


21 


3 


_ 


038 


.260 


3 


434 


3 


434 




262 


92 


(g) 


14 


21 


206 


- 


149 


.025 


5 


467 


5 


467 


- 


147 


93 


(g) 


13 


21 


3 




001 


- . 168 


7 


406 


7 


406 




169 


94 


Cd ; 


13 


-1 


3 




095 


. 010 


1 


636 


1 


636 




096 


95 


Cd; 


13 


1 


3 




029 


. 014 




988 




988 




033 


96 


(u ) 


13 


-2 


4 




538 


6.772 


-468 


861 


468 


910 




006 


Q7 




1 ^ 

lO 


_-i 
± 


n 




zoy 


Rn9 

. Ouz 


— 1 A"^ 

i*±o 


yoo 




QRQ 
yoy 






98 


(u-) 


13 


-2 


4 




050 


.087 


-1 


742 


1 


742 




101 


99 


(u) 


13 


2 


4 




149 


-.244 




988 




988 




286 


100 


(g) 


13 


21 


4 




204 


.198 


-2 


324 


2 


324 




284 


101 


(g) 


13 


21 


4 




021 


.084 


-1 


433 


1 


433 




087 


102 


(g) 


13 


21 


4 




078 


-.026 


-2 


324 


2 


324 




082 



99 



103 


(g) 


13 


21 


4 




001 




009 


-4 


412 


4 


412 




013 


104 


(u; 


13 


2 


4 




080 




216 


-4 


240 


4 


240 




230 


105 


Cu ; 


13 


-2 


4 




126 




018 
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o 

J , 


V d c 

, /5b 




H /I A 
. 140 


2/2 KO 


■7 

7 


311 


A 




C A -t 

. 041 




OCA 

. 3b4 


-b , 


/I A H 

. 401 


b , 


/ICO 




/I AO 

. 498 


0*70 «^ 

2^3 pi+ 


■7 

7 


211 






AOO 

. 032 




O H A 

. 219 


-2 , 


dec 

. bbb 


o 
2 , 


d O A 

, 580 




H /I A 

. 140 


274 K- 


7 


-321 







. 135 




. 027 


-5 , 


. 088 


5 , 


, 113 




. 494 


275 (.Targ-REM; 


11 


93 


226 




. 342 




. 955 


-20 , 


. 126 


20 , 


, 346 


2 


.810 


276 pi- 


7 


-211 







.063 




. 297 


-4 , 


. 395 


4, 


,407 




. 140 


277 Cf_2; 


17 


225 


275 




. 518 




. 135 


-18 , 


. 057 


18 , 


, 110 


1 , 


. 273 


278 pi- 


7 


-211 







.897 




.386 


-16 , 


. 770 


16 , 


, 799 




. 140 


279 pi+ 


7 


211 







. 078 




.064 


-6 , 


. 170 


6 , 


, 172 




. 140 


280 (Targ-REM; 


11 


93 


226 




. Oil 


2 


. 337 


-277 , 


. 032 


277 , 


, 048 


1 , 


. 826 


281 (omega) 


17 


223 


280 


_ 


.282 


1 


.678 


-154, 


.896 


154, 


,908 




.783 


282 (rhoO) 


17 


113 


280 




.293 




.659 


-122, 


.136 


122, 


,141 




.769 


283 ((gamma)) 


-7 


22 


285 




.076 




.029 


-16, 


.715 


16, 


,715 




.000 


284 piO 


7 


111 







.272 




.281 


-27, 


.169 


27, 


,172 




.135 


285 pi- 


7 


-211 







.065 




.210 


-17, 


,101 


17, 


,103 




.140 


286 pi+ 


7 


211 







.015 




.300 


-48, 


,046 


48, 


,047 




.140 




sum: 


.00 






.000 




.000 


0, 


,000 


900, 


.000 


900 


.000 
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E Summary list of subroutines 



c * 

c * 

C List of subprograms grouped according to general program structure * 

C * 

C (S = subroutine, F = function, B = block data) * 

C * 

C * 

C PART 1: Main steering routines and general utilities * 

C * 

C B VNIDATA to give default values to particle and decay data * 

C S VNIXRIN to initialize the overall run and simulation procedure * 

C S VNIXRUN to steer the simulation of a single collision event * 

C S VNIXFIN to finish up the overall run and simulation procedure * 

C S VNICBIN to interface parameters sets between commonblocks * 

C S VNILOGO to write program header * 

C S VNILIST to list event record, particle, or parameter values * 

C S VNIEDIT to perform global manipulations on the event record. * 

C S VNICOPY to save/restore copy of commonblocks & particle record * 

C F KVNI to provide integer-valued event information * 

C F PVNI to provide real-valued event information * 

C F RLU to provide a random number generator * 

C S RLUGET to save the state of the random number generator * 

C S RLUSET to set the state of the random number generator * 

C S VNIERRM to write error messages and abort faulty run * 

C S VNIFILE to open and close output files * 

C S VNIDUMP to dump (read) state of all commonblocks to (from) file * 

C * 

C * 

C PART 2: Initialization routines (global and process dependent) * 

C * 

C S VNICSIN to initialize collision system and kinematics * 

C S VNIABXS to give parametrized cross-sections for beam/target * 

C S VNITMIN to initialize time slizes and form of time increment * 

C S VNI2JIN to set up an initial 2- jet in 1+1- -> gamma/ZO -> qq~ * 

C S VNI2JFL to select flavour for initially produced qq~ 2-jet * 

C S VNI2JPY to get energy-momentum of initial 2-jet using PYTHIA * 

C S VNI4JIN to set up an initial 4-jet in 1+1 > W+W > qq~ q'q'~* 

C S VNI4JFL to initialize W+ cind W- decay vertices and the 4 jets * 

C S VNI4JPY to get energy-momentum of W+W- and 4 jets using PYTHIA * 

C S VNILHIN to set up initial system for 1- + h -> 1- + jet(s) + X * 

C S VNILAIN to set up initial system for 1- + A -> 1- + jet(s) + X * 

C S VNIXDIS to generate the variables of DIS photon-hadron (nucleus)* 

C F VNIGAQX to give elementary photon-quark cross-sections * 

C S VNIGAQQ to give qq~ fluctuation of photon from VDM * 

C S VNIHHIN to set up initial system for h + h'-> jets + X * 

C S VNIHAIN to set up initial system for h + A -> jets + X * 

C S VNIABIN to set up initial system for A + A'-> jets + X * 

C S VNIBTIN to set up initial state of hadronic or nuclear beams * 

C S VNINPDT to tabulate nuclear neutron/proton distribution * 

C S VNIRPIN to initialize single hadrons, or nucleons in nucleus * 

C S VNIQGDT to tabulate distributions in x and flavor of partons * 

C F VNIEFFQ to give parsunetrization of ef f . parton resolution Q(x) * 

C F VNIQGSH to give simple parametrization of parton shadowing * 

C S VNIVQKF to extract the valence quark content of a hadron * 

C S VNIQGXF to initialize x and flavor of partons in parent hadron * 

C S VNIQGCS to arrange a set of partons as color singlet hadron * 

C S VNIQGIN to distribute initialized partons inside parent hadron * 

C S VNIPRKT to give primordial kT to initial partons * 
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c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
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c 
c 
c 
c 
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c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 



S VNIXTIN to find spatial extension of initial particle system 

PART 3: Particle structure function package 
- portable, self-contained program unit with interface to PDFLIB - 

p , hyperon 



VNISTFU calls structure functions: photon, pi+, 

VNISTGA gives photon structure function 

VNISTAG gives anomalous part of photon structure function 

VNISTGS parametrized parton distributions in a virtual photon 

VNISTPI gives pi+ structure fimction 

VNISTFL gives proton structure function at small x and/or Q"2 

VNISTPR gives proton structure function 

VNIFCTQ CTEQ2'94 parametrization of parton distributions 

VNIFGRV GRV'94 parametrizations of parton distributions 

FV Function for GRV parametrizations 

FW Function for GRV parametrizations 

FWS Function for GRV parametrizations 

FHS Function for GRV pareimetrizations 

VNIHFTH gives threshold factor for heavy flavour production 

VNIDILN gives Dilogarithm function 

VNIGAMM gives ordinary Gamma function 



PART 4: Main routines for cascade evolution eind hadronization 



VNIEVND 

VNIEVDK 
VNIEV2J 
VNIEV4J 
VNIEVLH 
VNIEVLA 
VNIEVAB 
VNICOLP 
VNICOLX 
VNICOLI 
VNICOLL 
VNISBRA 
VNITBRA 
VNICHAD 
VNICFOR 
VNICREC 
VNICFRA 
VNIBHAD 
VNIBFRA 
VNIBMPS 
VNIBMMO 
VNIPADC 
VNIPAKF 
VNINSPE 
VNISTPG 
VNIPICL 
VNICDLH 
VNIHACL 



to administer general duties at the end of an event 
to decide whether event is o.k. (accept or regenerate) 



to perform evolution for 1+ + 1- -> 
to perform evolution for 1+ + 1- -> 
to perform evolution for 1- + h -> 
to perform evolution for 1- + A -> 

to perform evolution for A + B -> 



2 jets -> hadrons 
4 jets -> hadrons 
jets -> hadrons 
jets -> hadrons 
jets -> hadrons 



VNIALEM 
VNIALPS 



to find 2-parton collisions aind their kinematics 
to give cross-sections for 2 -> 2 and 2 -> 1 collisions* 
to set up a found 2-parton collision subsystem * 
to find outgoing flavors and type of parton collision * 
to generate space-like parton branchings * 
to generate time-like parton branchings * 
to orgcinize parton cluster formation and hadronization * 
to form one or two cluster (s) from a pair of partons * 
to reconnect cluster-partons according to color flow * 
to fragment found parton clusters into hadrons * 
to do beam remnant cluster formation and hadronization * 
to fragment beam/target remnant clusters into hadrons * 
to give phase-space dist. of hadrons from cluster decay* 
to balaince beam/target side momentum of hadrons 
to decay unstable particles according to particle data 
to give a hadronic particle from a flavor pair 
to recollect/fragment partons of spectator nucleons 
to simulate baryon stopping in nuclear A+B collisions 
to cluster pions together if above critical density 
to find & perform 2-body collisions of clusters/hadrons* 
to reconstruct cluster from its hadronic decay products* 

* 
* 

PART 5: Auxiliary routines I: * 

* 

- parton cascades & cluster formation - * 

* 

to give the alpha_electromagnetic value * 
to give the alpha_strong value * 
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C F VNIANGL to give the angle from known x and y components * 

C F INVCHGE to give three times the electric charge * 

C F INVCOMP to compress standard KF flavour code to internal KG * 

C S VNICOLB to boost/rotate 2-body collision subsystem along p_z * 

G S VNIDEC2 to generate 2-body decay * 

G S VNIFLAV to give flavor composition & beiryon number of system * 

G S VNIFRAM to perform boosts between different frames * 

G S VNIKLIM to find kinematic limits in 2-parton collsions * 

G S VNICCOL to connect color indices in 2-parton collisions * 

G S VNIJCOL to connect partons with color flow indices * 

G S VNILCOL to label color indices according to color flow tracing * 

G F VNIMASS to give the mass of a particle or parton * 

G F VNIMTHR to give 'current' or 'constituent' masses for partons * 

G S VNINAME to give the name of a particle or parton * 

G F VNIGORI to give the origin of a particle in space-time history * 

C F VNISORI to sey the origin of a particle in space-time history * 

G S VNIONSH to put final-state off -shell hadrons on mass-shell * 

G S VNIPROP to propagate a given number of particles * 

G S VNIPSEA to generate a hadron pair from the Dirac 'sea' * 

G S VNIQSEA to generate a q or q~, or a qq~-pair from the 'sea' * 

C S VNIROBO to rotate and/or boost a given number of particles * 

C S VNIRSCL to rescale particle momenta for numerical precision * 

G F VNIRCAM to generate raindom number from a 'ccimel' distribution * 

G F VNIRMTP to generate random number from a power-law distribution* 

G F VNIRMTX to generate rcindom number from 2-dim exp. distribution * 

C F VNIRPLT to generate random number from flat-top distribution * 

C F VNIRGAU to generate random number from Gaussicin distribution * 

G S VNISMRP to smear out classical particle positions and momenta * 

G S VNISDRT to count and sort particles * 

G S VNIWIDT to calculate width of particle resonances * 

G * 

G * 

G PART 6: Auxiliary routines II: * 

G * 

G - cluster-hadron decay and beam/target fragmentation - * 

G * 

G S HWHEPC to convert VNI event record to or from standard HEP * 

G S HWLIST to list event record in HEP format * 

G S HWIGIN to set parameters of HERWIG * 

G S HWUINC to compute constants and lookup tables * 

G S HWUINE to perform re-initialization of a new event * 

G S HWURES to load common blocks with particle data * 

G S HWUIDT to treinslate particle identifiers * 

G S HWCCUT to cut up heavy cluster into two clusters * 

G S HWCFLA to to set up flavors for cluster decay * 

G S HWCHAD to hadronize prepared beaim/target clusters * 

G S HWDHAD to decay unstable hadrons from beam/target remnants * 

G S HWDTWO to generate two-body decay 0->l+2 * 

G S HWDTHR to generate three-body decay 0->l+2+3 * 

G F HWDPWT to give matrix element squared for phase-space decay * 

G F HWDWWT to give matrix element squared for V_A weeik decay * 

G F HWMNBI to compute negative binomial probability * 

G S HWMLPS to generate cylindrical phase-space * 

G S HWMULT to pick charged multiplicity of beami/target diffraction * 

G F HWUPCM to compute cm. momentum for masses in 2-body decay * 

G S HWULDB to transform momenta from rest frame into lab fraime * 

G S HWULOF to transform momenta from lab frame into rest frame * 

G S HWUMAS to put mass in 5th component of momentum vector * 

G S HWUROB to rotate vectors by inverse of rotation matrix R * 

G S HWURDT to perform rotation with rotation matrix R * 

G S HWUSOR to sort n-tupel A(N) into ascending order * 

G F HWUSQR to give square root with sign retention * 

G S HWUSTA to make a particle type stable * 
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C S HWVDIF to give vector difference * 

C F HWVDOT to give vector dot product * 

C S HWVEQU to give vector equality * 

C S HWVSCA to give product of vector times scalar * 

C S HWVSUM to give vector sum * 

C S HWVZRO to give zero vector * 

C F HWRGEN to provide random number generator * 

C S HWRAZM to randomly rotate 2-vector pT * 

C F HWREXP random number from 1-dim. exponential distribution * 

C F HWREXT reindom number from 2-dim. exponential distribution * 

C F HWRGAU random number from 2-dim Gaussiein distribution * 

C F HWRINT random integer in finite interval * 

C F HWRLOG random logical * 

C S HWRPOW random number from power-law distribution * 

C F HWRPOT random number from 2-dim power-law distribution * 

C F HWRCAM random number from 'camel' distribution * 

C F HWRUNG random number from flat-top distribution * 

C F HWRUNI uniform random in finite interval * 

C S HWUTIM dummy routine * 

C S HWWARN dummy routine * 

C B HWUDAT to provide HERWIG common blocks with data * 

C * 

C * 

C PART 7: Integrated package WIZ * 

C * 

C - hadron / cluster cascade - * 

C * 

C All components of WIZ were developed and put together by the * 

C "Wizard of Poquott" Ron Longacre * 

C * 

C * 

C S WIZDATA to give default values to switches & parameters of WIZ * 

C S WIZCOLL to perform cluster / hadron collisions in a time step * 

C S WIZADQM to get cross section using the additive quark model * 

C S WIZSIGM to get cross section using the parametrized exp. data * 

C S WIZMMXS to evaluate total meson-meson interaction cross-section* 

C S WIZNNXS to evaluate total nucleon-nucleon cross-section * 

C S WIZCNQU to determine baryon, strange, charm aind bottom number * 

C S WIZCHCK to check energy -momentum, charge, flavor conservation * 

C S WIZXSEC to set up collsion of two (pre)hadronic particles * 

C S WIZYSEC to perform collision of two (pre-)hadronic particles * 

C * 

C * 

C S ANNIL to annihilate a baryon and antibaryon to a meson pair * 

C S MAKBBAR to make (create) a baryon- antibaryon from a meson pair * 

C S TNICES to generate outgoing particles of (pre-)hadronic coll. * 

C S ROTIT to rotate coordinate system for boosts * 

C S RBIAS to generate minimum bias event or beam jets * 

C S FXQUDD to determine strangeness number in HEP id code * 

C S FXCSBY to determine charm number in HEP id code * 

C S CHNGID to perform isospin rotation for minimum bias events * 

C * 

C S SCATS to do low-energy coll.'s - meson-nucleon, Delta-N * 

C S SCATKN to do low-energy coll.'s - K-N * 

C S SKATH to do low-energy coll.'s - K-Lambda, K-Sigma * 

C S SWAVE to rotate colliding particle momenta creating a S-wave * 

C F BARPRO to get probability for Delta-Pi or Rho-N channels * 

C F BCHPRO to get probability for 3-body decay out of collision * 

C S DELKST to determine K* or Delta modes in K-N collision * 

C F ETAPRO to get probability for eta production * 

C F KNPRO to get probability for K-N production * 

C F KSTPRO to get probability for K* production * 

C F LAMPRO to get probability for Lambda production * 
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c 


F 


LP I PRO 


c 


F 


PKNPRD 


c 


F 


RHOSC 


c 


F 


SCPPRO 


c 


F 


SCPRO 


c 


F 


DELPR 


c 






c 


F 


SIGBR 


c 


F 


SIGLAM 


c 


F 


SIGPRO 


c 


S 


SIGRHO 


c 


F 


SPIPRO 


c 


F 


SRHPRO 


c 


F 


SWAPRO 


c 


S 


RECPRO 


c 


F 


BCHKH 


c 


F 


LAMKH 


c 


F 


SIGKH 


c 






c 


S 


MMSCAT 


c 


S 


SCAM 


c 


F 


SCPPI 


c 


F 


SCPI 


c 


F 


XEHC 


c 


F 


CHEX 


c 


F 


ETAETA 


c 


F 


KKBAR 


c 


F 


PIPI 


c 


F 


ETPI 


c 


F 


PIK 


c 


F 


BCHTB 


c 


S 


RECP 


c 


S 


ISOSCAT 


c 






c 


F 


SCPAl 


c 


S 


SCAIT 


c 


F 


PIAl 


c 


F 


RAIPI 


c 


F 


CHEXPI 


c 


F 


KKB4PI 


c 


S 


MAKISO 


c 


S 


AlPRB 


c 






c 


S 


PUTPNT 


c 


s 


GETPNT 


c 


F 


CHARJ 


c 


F 


AIMASS 


c 


F 


PBARGT 


c 


S 


ANMESID 


c 


F 


PANNIL 


c 


S 


BBARID 


c 


F 


YSTPRO 


c 






c 


S 


PHIEXT 


c 


S 


PHIEXTM 


c 


F 


PHIPI 


c 


F 


PHIK 


c 


F 


PHIKST 


c 


F 


PHIETA 


c 


F 


PHIPHI 


c 


F 


PHIPIAl 


c 


S 


CHAGID 


c 







to get probability for Leunbda-Pi production in K-M 

to get probability for pi K-N production 
to determine if F.F. for rho-N is called 
determines if F.F. is called for non-MM scattering 
determines if elastic scattering for non-MM channels 
determine the probability for Delta M -> M N 



get branching ratio for Sigma-Pi charge states 
get branching to K-Lambda from Pi-N 

determine if Sigma is produced in low energy M B coll. 
to determine Sigma* or Rho modes in K-M collision 
determine Sigma pi for low energy K-N collision 
determine Sigma pi pi for low energy K-M collision 
determine three modes from low energy K-N collision 
to determine Delta or Rho modes in Pi-N collision 
determines if three body in Pi-Hyperon interaction 
determine if Lambda is produced in K-Hyperon scattering* 
determine if Sigma is produced in K-Hyperon scattering * 

* 

to collide 2 mesons and generate outgoing particles 
to do low-energy coll.'s M M scattering 
determines if F.F. is called M M scattering 
determines if elastic scattering for M M 
determine charge exchange for K* system 
determine charge exchcinge low energy Pi Pi K Pi scatt. 
determine if eta eta is produced from M M scattering 
determine if K Kbar is produced from M M scattering 
determine if pi pi is produced from M M scattering 
determine if eta pi is produced from M M scattering 
determine if K pi is produced from M M scattering 
determines if three body in M M interaction 
determine delta or rho isobars in M M collision 
collide isobar states N* Y* Al ect. with isobar or reg 



determines if F.F. is called Al meson scattering 
to do low-energy coll.'s - Pi-Al 
determines if Pi-Al is produced 
determines the charge mode for Pi-Al 
determines the correct charge for the outgoing Pi Pi 
determines if K Kbar should be produced from Pi-Al 
makes isobars out of 2 particles 
to determine probability for Al formed from 2 particles* 

* 

reuse array position that stores isobar decay products 
gets a reused position to store isobar decay products 
returns charge of particle cind works for isobars 
returns mass of particle and works for isobars 
determines energy dep. for diquarks at Ecm in F.F. 
B Bar annihilation into mesons the meson id's are made 
determines the B Bar annihilation probability, 
to decompose baryon id into flavor/spin of a baryon 
determines Y* production. 



determines phi production in non M M scattering, 
determines phi production in M M scattering, 
determines if a phi pi state is produced in M M scatter* 
determines if a phi K state is produced in M M scatter * 
determines if a phi K* state is produced in M M scatter* 
determines if a phi eta state is produced in M M scatt.* 
determines if a phi phi state is produced in M M scatt.* 
determines if a phi pi state decays from the rho prime.* 
Puts in the correct OZI suppression for phi B scatters.* 

* 
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C F IDISA to convert particle id's from HEP to WIZ codes * 

C F IDHEP to convert particle id's from WIZ to HEP codes * 

C F HEPMASS to provide particle masses according to HEP standard * 

C F RANF to provide internal random number generator for WIZ * 

C * 

C S MBSET Set parameters for generating minbias or beam jets. * 

C S MDVLEV Moves one array into an other array. * 

C S DECAY The WIZ decay of particles not used outside of WIZ. * 

C S FLAVOR The flavor decode of WIZ id's. * 

C F LABEL Returns character*8 string for WIZ id's. * 

C F AMASS Mass of particles with WIZ id's. * 

C F CHARGE Charge of particles with WIZ id's * 

C F DBLPCM Calculate com momentum for A->B+C with double precision* 

C F GETPT generates pt for quark pairs for the F.F. chain. * 

C S IDGEN Does nothing. * 

C S L0REN4 Does a general Lorentz transformation. * 

C S JETGEN Fragments jet into primary hadrons using Field Feynman.* 

C S PRTEVT Prints the event stored in the ISAJET part of WIZ. * 

C S QCDJET Final state QCD jet evolution using Fox and Wolfram. * 

C S READIN Setup parameters for the ISAJET part of WIZ. * 

C S RESCAL Rescales momentum of particles to the total energy. * 

C S RESET Resets all user defined variables for ISAJET part of WIZ 

C S SETCON Sets the constants for the ISAJET part of WIZ. * 

C S SETDKY Sets the decay table that is used inside WIZ. * 

C S SETTYP Sets jetype flags in ISAJET. * 

C S HEAVYX Generates X for heavy particle fragmentation (Peterson)* 

C S QCDT Auxiliary routine for QCDJET. * 

C S QCDZ Auxiliary routine for QCDJET. * 

C * 

C * 

C PART 8: Event study and data analysis routines * 

C * 

C S VNIANAL to steer pre-programmed event analysis * 

C S VNIANAl statistics on global particle properties & observables * 

C S VNIANA2 statistics on time evol., ln(l/x) spectra, BE effect * 

C S VNIANA3 statistics on pre-hadronic clusters * 

C S VNIANA4 statistics on partons and produced hadrons * 

C S VNISPHE to perform sphericity analysis * 

C S VNITHRU to perforin thrust analysis * 

C S VNICLUS to perform three-dimensional jet cluster cinalysis * 

C S VNICELL to perform jet cluster ciaalysis in (eta, phi, E_T) * 

C S VNIJMAS to give high and low jet mass of event * 

C S VNIFOWO to give Fox-Wolfram jet moments * 

C * 

C * 

C PART 9: Dummy routines and system dependent utilities * 

C * 

C S VNIPRIN dummy routine: remove when linking 'vniprob.f ' * 

C S VNIEVPR dummy routine: remove when linking 'vniprob.f * 

C S VNIPRND dummy routine: remove when linking 'vniprob.f * 

C S VNIPRIN dummy routine: remove when linking 'vniprob.f * 

C S PYINIT dummy routine: remove when linking ' vnipythia. f ' * 

C S PYEVNT dummy routine: remove when linking 'vnipythia. f ' * 

C S PDFSET dummy routine: remove when linking PDFLIB * 

C S STRUCTM dummy routine: remove when linking PDFLIB * 

C S VNICLOC to measure and return the CPU time used (sys. dep.) * 

C S VNIDATE to give current date and time (sys. dep.) * 

C * 

C * 
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